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03/31/80 


The Soft Wsirehouse 
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4 .. 
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11/01/79 

12 

PLES5.TRA 

11/01/79 
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FILE mmc casfjmncti 


1. Files of type COM are unprintable directly-executable machine* 
language program COMmand files. 

2. Files of type SYS are unprintable machine*language memocy^images 
which can be LCADed from within muSIMP. 

3. Files of type DOC are printable files Documenting the usage of a 
program file having the same first name. 

4. Files of type TXT are printable text files containing information 
implied by the first name. 

5. Files having a first name of the form CLESn are interactive 
Calculator-mode LESscxis to be executed from within muSIMPr in the 
order indicated by numeric moffix n. 

6. Files having a first name of the form PLESn are interactive 
Programming-mode LESsons to be executed from within muSIMP, in the 
order indicated by the numeric suffix n. 

7. All files listed above of a type other than COM, SYSr DOCr and ^CT 
are muSIMP program source files. The type name denotes the first 
three letters of the first name of the most immediate prerequisite 
program file. For example: 

a) Tt» nujSIMP file named Alj5£BRA.ARl implemaits adgebra, requiring 
the muSIMP file named ARITH.MOS as a prerequisite. 

b) File ALSEBSAJXX: is the reference documentation for usage of the 
f^ilities implemented by file ALSIBPA.ARI. 

c) Files CLES3JVLG and CLES4.ALG are interactive calculator-mode 

lessons teaching use of the facilities implemented by file 
ALGSBPA.ARI. 

% SDS (}$ 
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File READlSr.TSCT (c) 


10/30/79 


Tbe Soft Waxebouse 


General infocnation 

Cor^ratulations <xi your purchase of the iDuSIMP/muMAXB-79 Symbolic 
Mathematics System. This package is a revolutioiary and s^histicated 
software s^t&a. for 8080 a^ Z80 based microcomputers. Therefore some 
degree of study and patiaice is rec^ired to properly build a ^tem from 
muMATH source filesr and then save the result as a memory image file. 
However once built and saved, it is a simple matter to load and interact 
with the systea at any later time as described in file INI!ESACr.ascr. 

Unfortunately your first task is the rather uninteresting one of 
building the system as explained in the remainder of this file. The 
immediately following background information is provided to make the 
process seem less inysterious. 


SasL miSSSSzlS^ giggranning, language 

^ muSIMP-79 (micro-<;omputer Structured iMPlanentation language) 
systffis provides a high level programming language waitable for a wide 
variety of ai^lications. It is implemented using an efficient and 
versatile interpreter requiring 7K t^es of machine code. Ttm current 
version of muSIMF also requires a bootstrap file which is loaded 
immediately after the machine coded portion. The interpreter is 
distributed as two disk files: 

MUSlMP79.CX3i an executable OCnnand file 

HUSMQRE.MUS the muSIMP bootstrap file 


The file named FILES.!Z3a' indexes and briefly describes the imoSIMP 
documentaticsi package. The documenation is only distributed in printed 
form. FILES.TXT also lists the lesson files required to become 
proficient in the programming language. Since using muMATB in the 
calculator mode requires no programming knowledge, learning muSXMF can 
be safely postponed. 


^ muMa3a~79tm 




UsL mih sss&m 


The muMATH-‘79 System provides the facilities to perform a wide 
variety of symbolic mathematical ^rations efficiently and accurately 
on a axoputer. It is implemented as a set of muSIHP program packages. 
These source file packages are organized in a very modular fashion in 
order to accommodate both differing mathematical needs and differing 
computer memory sizes. More sophisticated mathematical packages retire 
prerequisite files as indicated by the following dependency diagram. 
Each file requires those abov% it in the diagram as a prerequisite. (The 
significance of the numbers will be explained below.) 
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MDSH-IPTS.CCM 


/ 

/ 


\ 

\ 


aSACE.MOS 


ARI1B.MUS 


H/A,551 1536,3115 


/ 

y / 

ARRAY.ARI 


\ 

\ 

ALSEBR^.ARI 


375,756 1160,1946 

/ I 

/ --- 

V y / I I 1 

HAS3tlX.ARR EQN.ALS | LG6.AI£ | 

497,892 86,208 I 186,345 I 

J / I I 

SaVE.BQN TFGPOS.ALG 7SGME1G.ALG 
452,734 327,602 372,680 


\ . 

DIF.AIJG WO-i, 
306,580 

INT.DIP 
794,1420 
I ^ 
QCSOlE.QZr 
914,1760 


System Generation groce^ilfi 

If you are proficient in the use of the microcomputer's disk 
operatin 9 system (DOS), the foliating procedure should be sufficient 
explanation to build and save a. complete muMAIH system. However if you 
are a novice or questions arise, additional information can also be 
found in files BACKUP.TXT and INTERACI.TXT, and in the documentation 
provided with the disk grating system. 

I. G&ierate a muSBQ’/muMAIH-79 badcup disk. 

Using the computer's DOS, transfer a copy of the source 
files depicted in the above diagram from the Soft Warehouse 
master diskette to a backup diskette. Since the total disk 
^>ace required to store the files is approximately 96 K bytes, 
more than one diskette may be required. Often it is 
convenient to generate a DOS on the backup disk(s}. 

II. Build a riuSlMP-79 System. ^ ' S 

Execute the MUSIMP79 COMmand file by the entering the 
following DOS command: 

I4US1MP79 

It Should then di^lay a versiof\/copyright message. Jot 
down the number following the word "SAVE” in the message for 
later use. The bootstrap file will begin to load 
automatically. After about 5 minutes load time, the system 
should respond with the ”? ” prompt character indicating that 
nu:SII4F has been successfully constructed. 


2 



III. Build a iiuMftIH-79 Systan. 

Two muSIMP commands ace retired to build a muHATE 
^tem. All commands must be tecminated by a semicolon and a 
carriage return in order to initiate tbe action. 

Ihe source file "nameVtype” on disk ’drive* is loaded 
by issuing a command of the form 

BDS (name, type, drive); 

Note that the file name and t^ are separated a comma 
rather than a period, and the drive requires only a single 
letter without the customary colon. The drive is an option^d 
argument which defaults to the drive currently logged in. 

Each muMATH source file requires a given amount of 
computer memory to store its function definitions. The unit 
of storage in muSlHF is the "node” which is described in file 
MDS>A3AJ1US. The pairs of lumbers in the dqpeiH3«icy diagram 
above are the approximate number of nodes required to store 
the corresponding muHATH package in both condensed and 
unccndensed form. See the next section for the significance 
of the smaller ccxidensed nuxsbers. 

The number of free {i.e. unemployed) nodes is given by 
the muSXMP command 

RECLAIM {); 

When building a system you must ensure that there is 
adequate storage for the program plus, at least 500 additicxial 
nodes to store mathematical expressions. This can be 
accomplished through the use of ^iCLAIM and r^ersue to the 
unccndensed number in the dependency diagram corresponding to 
the program package. Once this has been verified, ececute the 
appropriate EDS load command. This procedure is illustrated 
by the following dialogue building a muMA3H ALGEBRA systsn: 

? RECLAIM 0 ; 

q 6360 

? BDS (ARITH, MUS) ; 

@ ARUH 

? RECLAIM 0 ; 

€ 3245 va.i-v; - 

? BDS (ALGEBRA, ARI) ; 

@ ALGEBRA 

? RECLAIM 0 ; 

q 1300 ^ ■ 
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IV. Save the i!uMA 3B~79 Systea. 

In order to avoid repeating the tedious build routine in 
the futuref a memory image file of the muMAIH system should be 
saved immediately after the build using either method A or B 
described below. Whichever method is used, it is advisable to 
use a file name representative of the most sophisticated 
mathenatical capabilities loaded. For example, names such as 
ALSIBBA, miG, CALCULaS, EQUAnCN, etc. 

*** WARNING *** Before attempting a save, ensure that a 
diskette with sufficieit free ^ace to store almost an entire 
ffi«nory image is properly mounted on a disk drive. 

A. Return control to the DOS by typing CTRL-C. 

Then generate a COM type file by issuing the DOS 
SAVE commemd. Use the decimal number, N, of 256 
byte records recorded in step II above. The file 
size will be (N/4) K tytes. 

Although most desirable from a convenience standpoint 
(see file INTERACT.TXT), this method may NOT work due to 
limitations of the DOS. The easily recognizable syn^oms are 
either the inability of the DOS to save that large a file, or 
erratic behavior when the file is subseguertly executed. 
If this happeis, you will have to build the ^stem again from 
step II and use method B in the future. 

•I 

B. Type the muSIMF ccnmand of the form 

SAVE (name, drive); 

This will save a file called name.STS on the given 
drive. Again the drive is optional, defaulting to 
the current drive. The SYS file will be 
approximately (N-28)/4 K bytes, where N is the 
number recorded in step n. 




Generatioi of CCXi or SYS files is also useful for checkpointing a 
lengti^ dialogue for purposes including: 

1. coitinuation at a later time, 

2. preservatioi of an environment so that uncertain exploratory 
ccn^tations which might endanger the envircnmsrt can be 
safely pursued, 

3. preservaticn of a "program” (meaning an environment) produced 
interactively rather than using a text editor. 

O 
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Cflndenggd Systea Ggngration 


Although muMATH is much more compact than any previous general- 
purpose symbolic math system, it is still a large set of programs for 
microcomputers. Fortunately, the use of a "condensing" technique to 
load programs can ecoiomize on me&ory consumption by a factor of almost 
two. If the control variable named CCNDEllSE is assigned the value TRUE, 
then common subexpressions of function definitions are autcmiatic 2 my 
shared as the source files are read iiv infix muSIM? colon aerator 
is used to make the assignment to CONDENSE as described in the muSIHP 
lesson and documentation files. Thus, to build a condensed ALGEBRA 
system, the following commands should be issued in place of ti»»e in the 
example of stqp II above: 

CONDENSE: TRUE; 

RDS (AHUB, ICS); 

ms (ALGIBRA, ARI); 

CCNDENSE: FALSE; 

The exhaustive searching makes condensation too slow for inter¬ 
active use, so that is why it is advisable to set CCNDQiSE to FALSE just 
before the save. The above condensed load requires an hour or so, 
depending on the processor speed. Since it is possible to type the 
above commands on one line, you are free to take a long break while the 
ocxidensatiai takes place. 

Condensation can be regarded as an qptional sort of "compilatioi" 
stage, and one of the principal reascxis for generating COH or SYS files 
is to preserve this investment of time. However, for simplicity, we 
suggest not worrying about prodiibing condensed COM or SYS files until it 
becomes impossible to otherwise fit all of the desired files into memory 
smiltaneously. 

Many of the tmoMATB files contain optional sections identified by 
con^icuous comments. To save space, appropriately pruned versions of 
these files can be created using a text editor. For example, a user 
with oily 32 kilot^es may have to delete portions of the arithmetic and 
algebra packages in order to do algebra comfortably, even with 
c^deisation. As with the set of muMAIE files, each file is internally 
organised "bottom up", with the most expendable and highest-level 
features late in the file. Consequently, the files can be truncated 
between almost any two commands, without risk of invoking undefined 
functions or uninitialized cmtrol variables. 
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%File; INIERACT.ECT (c) 10/30/79 The Soft Warehouse 


miTIATINS A ibuHAZE CCH FILE 

Initiating muMAra is easiest if someone has already saved a O^tmand 
file having a memory image including all of the muMMH or other programs 
which you wish to use at that time. The name of such a file is of 
course up to wt»ever creates it, but in general the name assigned is of 
the highest-level math package loaded. Thus MJSEBBli,OOH would be the 
name of a command file containing ARTBLMUS and ALGEBRAAEI. 

For example when using a CP/M (tm) type (grating systaa with file 
AXiGE2PA.COM on the current drive, one merely enters the following 
system-level command terminated by a carriage return: 

ALGEBBA 

After a minute or so of loading from the diskette, the re^nse should 
be a message of the form: 

muSIMP-79 (Versiem roenth/day/year) SAVE: sise 
Copyright (c) 1979 by The SCFT WAREHOUSE phone 

where appropriate numbers appear as the «itries "month", "day", "year", 
"size”, "^one". You are now free to enter mathematical expressions 
as described below. .. 


INITIATING A muMATH SYS FILE 

Qifortunately for reascxis described in file READlST<T3Cr, it be 
ispossible to construct a muMAlH GOI4mand file. However, the following 
means of initiating muMATE is always possible and quite easy provided 
someone has saved a SYS-type file containing a memory image of the 
mul'lATB packages which are needed at that time. Here too the name of such 
a file is up to whoever creates it, but in all probability the same 
naming convention describe above for COM files was used. Therefore 
ALSEBRA.SYS would have had the same source files loaded in as the file 
ALGEERAXCM would have. 

For example, if HUSIMF79.COM is on the oirr^t drive and ALGEBRA.SYS 
is on disk drive B, the appropriate operating system load command would 
be as follows: 

MDS1MP79 B:ALGEBRA 

About half a minute after the muSIMP log^ message appears, the mMATH 
system should respond with the "? " prompt characters. Now you can 
begin your interactive dialogue with muMATH. 
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USE INHSACnCN OCLE 

nuSIMF prompts tiie user with a (^estion mark indicating readiness to 
accept a command entered f rom the terminal. The user then types an 
ei^ression followed by a semicolon and a carriage return. First muSDIP 
parses the expression and converts it into an internal representation. 
After printing an to herald the ”@nswer"r the expression is 
evaluate and then a ^pace is printed to indicated the evaluation phase 
is conplete. Finally the result is d^arsed and printed in mathanatical 
notation. This interaction ^cle is repeated indefinitely until a CZE^'-C 
is typed (i.e. a "C” typed while depressing the CTBL key). 

For example, here is a segment of a trivial muSXHP dialogue: 

7 5; 

€ 5 

7 2 + 2 ; 

e 4 

7 JCES » MAR^; 

@ FALSE 

7 ME-BER (APPLE, ' (GRAPE, APPLE, PLUM)); 

e TSVE 


CCRRECmC 1TPCGRAPBICAL ERRORS 

Since muSIMP uses the operating system's console 1/0 routines, all 
the line-editing features of that system are inherited by muSIMP. 
Backspacing is usu^dly accomplished by typing either a CTRL-B, or a 
RUBout, or a DELete key. Some systems echo the deleted character; 
whereas, others erase the character from the screen and backspace the 
cursor. Entire lines can be deleted or flushed by typi^ a CZRL-U or a 
CTRL-X. As a note of caution: there is no way to modify a line once a 
carriage return has been-typed. If this happais, ths ^ire e:^ressicri 
can be flushed by typing a senicolon. 


INZERBUPmC EVALUATION 

An evaluation in progress can usually be aborted by typing a CIRL-Z, 
Escape, or ALTmode. An options available message will then be 
displayed. The usual choice is typing another CTRL-Z, ESCape, or 
ALTmode which allows you to enter expressions as before. The other 
alternatives are fully explained in file ERBCRS.TXT. As a last resort, 
the computer can of course be RESET and a "cold start" performed to 
reload the operating syst^ 

% SDS (} $ 
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File: BACKQP.ISCT (c) 


11/01/79 


The Soft Warehouse 


How to Backu|> the Master Diskettes 
allied by the Soft Warehouse 


1. The following information is provided for those who are not 
throughly familar with their computer's disk grating system. Since 
there are many different ^)erating s^tems and computer configuratioiSf 
it is necessarily a gei^ral guide which should be supplemented by study 
of the documentation ^i^lied with the disk operating system. 

2. Obtain an a^ropriate number of blank new, high-quality diskettes 
suitable for your drive. 

3. Become thoroughly familiar with the terminal, computer, disk drives, 
and operating system. Most cases of accidental erasure of the master 
diskettes or oti»r irrecoverable errors are committed in tl:» first few 
moments by eager users, inexperienced with the system on which th^ aure 
installing the new software. In particular, practice initiadizing a 
diskette, then generating a disk operating system on it. Use the 
largest version of the system in terms of the space available to user 
programs. Finally transfer to the new diskette files from a spare, 
write-protected diskette. 

4. Due to wear and inadequate industry-wide manufacturing sta n da r ds, 
there are sli^t or not-so-slight'mechanical and electrical differences 
between various nominalli^ compatible drives and diskettes. 
Consequently, we suggest that if you have two or more drives, you first 
try placing the new diskette in the drive on which it will be used most 
often. The Soft Warehouse diskette can then be tried on each of the 
other drives, trading with the new diskette if ncxie of the other drives 
are successful. 

5. For two or more drives, some grating systems provide a convenient 
COFY DISK command which automates most of ti% copying protocol otherwise 
necessary. Alternatively, the PIP (Peripheral Interchange Program) or 
XFEE (File Transfer) commands of most operating systems usually permit 
copying all files item drive A to drive B by a command such as 

PIP A:*«;*.* 

or XFER A:^;*.* 

6. Copying diskettes is much more laborious on systems having only one 
drive. Gaieradly, it involves repetitively reading a portico into main 
memory from the old disk, switching disks, then writing the portioi onto 
the new disk, then switching back to the old disk. The process 
generally involves using a resident bootstrap monitor in read-only 
memory, or using a "DDT"-like "Dynamic Debugging Tool" program. 
Moreover, the process may pad out the otherwise oily psu^ially filled 
last page or record of a file with au±>itrary garbage w^ch is harmless 
in a commamd file but annoying in a text file. Thus, text files 
transferred this way may require some text editing to clean-up. 
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7. For some operating systems it may be necessary to remove the write 
protect tape from the old diskette temporarily, even though it is only 
to be read from. (Perhaps this is true only if the old diskette is in 
the "principal" drive.) Moreover, it be awkward to copy a diskette 
without first removing the write protect tape in order to copy the 
operating system onto that diskette, especially if there is only one 
drive. 

8. We use costly highest«quality diskettes, and we endeavor to record 
them on the most precisely adjusted drives available. C^iseguently, if 
you cannot read our diskettes after several attempts and after carefully 
restudying our directions and those provided for the hardware and 
operating systen, then 

a) Carefully check whether or not you correctly specified 
all of the details on the order fora and purchased the 
proper type of blank diskettes. 

' b) Use an alignnent test diskette or have your drive 
checked professionally. 

c) Get help fron an experienced professional or friend. 



2 



%FiIe: LESSCNS.OXT (c) 


10/30/79 


Tbe Soft Waurehouse 


Osiiig Ibe interactive Lessons Files 

This file explains how to use the interactive muSIMP and muMATE 
lesson files. muMATH and the lessons are designed to serve a broad 
range of math levels from arithmetic through calculus, and to serve a 
broad range of programming backgrounds from none to professional 
programmer, low is this scope possible? Read oi: 

The sequence of lesson files CLESl, CLSS2, etc. explains how to use 
muMAIS as an arithmetic or symbolic calculator, for successively more 
sophisticated mathematical (perati(»s. The sequence of lessons FLESl, 
PLES2, etc. explains how to write programs in muSIMP, in order to 
enhance the suite of built-in operations or for ai^ other purpose. 

The calculator sequence is ordered according to the most common 

sequence in which the corresponding math subjects are taught. It is 
intended that a user proceed in this sequence only as far as their math 
background, before ^ionally beginning the programming sequence. Due 
to slight variations in math curricula sequences, some users may prefer 
to skip certain calculator lessons in the middle of the sequence as well 
as at the end. 

muMATH has such a rich set of built-in capabilities that many users 
will be content to postpone study of the programming sequence 
indefinitely. However, many users evenntually will want to proceed to 
the programming sequence, perhaps for one or more of the following 
reasons: -> 

1. to eihance the built-in mul^ATH capabilities, 

2. to understand how the underlying muMATS algorithms work, 

3. to learn cooputer programming, 

4. to use muSIMP for seme other a^lication. 

In order to make the programming sequence most useful to users of 
all mathenatical backgrounds, the sequence begins with muSIMP examples 
which are non-mathematical, or arithmetic at most. Host general 
programming techniques and their realization in muSIMP are indep^ent 
of higher-level math. Ttus, only the last lessons in this sequence deal 
with muMATH specifically, explaining how to extend it, alter it, and 
even replace it with alternative symbolic math systems. 

There are three ways to experience the lessons. For most people, 
the best way is to execute them interactively, trying out examples at 
the opportunities provided; the second best way is to read the printed 
record of a dialogue produced by someone else executing the lessois; and 
the third best way is to read the files containing tl:» original lessons, 
which contain only one side of an inteivied dialogue. 


1 



As iiidicated in file F1LES.!]XT, the first lesson is file CLESIJ^RI. 
Consec^ently, to cosun^ice the lesson you must first initiate a nuNAIlE 
system co ntai ning at least file ARUHJIUS. Bow to do this is explained 
in file ornUfiCIJISCI. Then, you simply issue the muSIMP command 

PDS (CLE31, ARI, drive); 

where "drive" is the name of the drive on which the CLESl.Am is 
mounted. The lesson will tell you what to do from then on. If the 
interactive lesson for any reason becomes hopelessly confusing, you can 
always cease the lesson and simply read it. Also, it may help if you 
take the lesson along with a companion, because your possible confusions 
may be disjoint. 

Have funi 

% RDS 0 $ 







%FUe: CLES1.ARI (c) 


10/30/79 


The Soft Warehouse % 


LOtELOlTni (78}$ «ECaO: ECSOS ECBO: TESTES 

% If this lesson is b eing display too fast, it can be temporarily 
stop ped by t^ing a CIRL-S (i-e. typing the letter "S" while depressing 
the CTRL icey). Then type it again vhai you are reat^ to resume. 

If yaj have not yet read files LESSCNS.TKT and INTEHACT.T2CT/ it is 
advisable to ^rt this lesson and read those files first. To abort the 
lesson, enter an ESCape or a CTRL-Z character followed by a CITL-C. 

In shiMATB a "comment” is a percent sign followed by any number of 
other characters terminated by a matching percent sign. Thus, this 
explanation is a comment which has not yet been terminate Commaits do 
not cause computation; they are merely used to explain programs and 
examples to human readers. Here is an example of an actual computati<»% 

1/2 + 1/6 ; 

% Note how luiMAIH uses exact raticxial arithmetic, reducing fracti^is 
to lowest terms. 

In muHATH, arithmetic expressions can be formed in the usual 
manner, using parentheses together with the operators 
"/", and respectively for addition, subtraction or negation, 
multiplication, division, and raising to a power. Et>r example: % 

(3*4 - 5) " 2 ; 

% some terminals, looks like an upward-pointing arrow; on 
others it looks like a shallow upside-down letter V; and some terminals 
may employ an utterly different looking character which you may have to 
determine by e^qperimentation. 

The reason for using and * is that standard terminals do not 
provide superscripts or centered dots or special multiplication crosses 
distinct from the letter X. 

To prevent certain ambiguities, railtiplicatim cannot be implied by 
mere juxtapositicsu Gne of the most frequent mistakes of beginners is 
to omit asterisks. 

Later, in order to give you an opportunity to try some examples, 
we will "assign” the value FALSE to the variable named RLS. When you 
are ready to resume the lesson, tht "assignm^t” 

BDS; TKJE ; 

ii^udi^ the semicolon and carriage return. This revises the value of 
the variable named RDS to the value TPETE. We will explain assignment in 
more ^tail later. 

Don't forget that yai can use local editing to correct mist^ings 
on the current line. For example, on many operating systems, the key 
marked lUTBout or DELete cancels the last character typed on the line, 
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and typing a CZRL'-U cancels the current line. There is no way to modify 
a line after striking the R£Tum key, but an expression can always be 
flushed by typing a final line containing a "grammaticeLl" or "syntax" 
error such as 

Now we are going to turn control over to you by setting RDS to 
FALSE. Try some examples of your own similar to the above. Also we 
suggest that you make a few intentional errors in order to become 
familiar with how they are treated. For example, try 

5 7; 5+ /It 5/0; and 0/0; 

Havefuni: % BDS: FALSE ; 

% The value resulting from the last ii^ut expression is automatically 
saved as the value of a variable named #ANS, which can be used in the 
next eaqpression. For example: % 

3 ;«ANS * #ANS ;«ANS " IANS; 

% As this acample illustrates, muMATB can treat very large numbers 
exactly and quickly. In fact, muMAIH can accomodate numbers to about 
611 digits. To partiedly appreciate how large this is, compute the 
distance in feet or in meters to the star Alpha Centauri, which is 4 
light years away, then use IANS to compute the distance in inches or in 
centimeters without starting all over. (In case you forgot, the speed 
of light is 186,000 mile^second or 300,000,000 meters/secon^} % 

RDS: FALSE ; 

% Our answers are about 123,883,499,520,000,000 feet or 
1,486,601,994,240,000,000 inches or 37,843,200,000,000,000 meters or 
3,784,320,000,000,000,000 centimeters. Another dramatic comparison with 
10''611 is that there are thought to be aU^out 10"72 electrons in the 
entire universe. (Whoever oxjhted them must be exhausted!) 

Often one performs an intermediate computation or a trivial 
assignment for which there is no need to display the result. Whai this 
is the case, the display of the result can be suppressed by using a 
dollar sign rather than a semicolon as a terminator. For example, type 

RDS: TEDE $ 

a;nd note the difference frcm when you previously typed RDS:TRUE ; % 

RDS: FALSE $ 

t It is often conveueit to save values longer than IANS saves them, 
for use beyond the next input expressi^ The colon ASSISlMEIir operator 
provides a means of doing so. The neme on the left side of the 
assignment operator is BOUM) or SET to the value of the expression on 
its right. This value is saved as the value of the name until the name 
is bound subsequently to some other value. The name can be used as a 
variable in subsequent expressicns, as we have used #ANS, in which case 
the name coitributes its value to the expression. For example: % 

RATE: 55 $ TIME: 2 $ DIST?^: RATE * TIME ; 

% Alphabetic characters irclude the letters A through Z, both u^per 
and lower case, and the character *#”. Note that the upper and lower 
case version of a letter are entirely distinct. Names can be smy 
sequence of alphabetic characters or digits, provided the first 
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character is alphabetic. Thus X, *9, and ABC3 are valid names. Make an 
assignment of 3600 to a variable named SECPEKHCXJRr then use this 
variable to help compute the number of seconds in 1 d^ and 1 week: % 
roS; FALSE $ 

% Congratulations on completing CLES1.ARI. To execute the next 
lesson, merely enter the muMATH command 

FDS (CL£S2, ARIr drive}; 

where drive is the name of the drive on which that lesson is mounted. 
Alternatively, it may be advisable to repeat this lesson, perhaps 
another day, if this lesson was not perfectly clear. The use of any 
computer program tends to become much clearer the seccnd time. 

In order to experience the decisive learning reinforcement afforded 
by meaningful persal^U. examples that eure not arbitrarily ccntrived, we 
urge you to bring to subsequent lessons appropriate examples from 
textbooks, tables, articles, or elsewhere, i^o, you are encouraged to 
experiment further with the techniques learned in this lesson: % 

ECSO: »ECHO $ 

ms 0 $ 
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%FUe: CLES2.ARI (c) 


10/30/79 


^ Soft Warehouse % 


LINELEtGIE (78) S *EC50: ECSOS ECHO: TEDES 

% This file is the secmd of a segueice of interactive lessons ibout 
the muMJ^-79 system for computer symbolic math. This lessoi presumes 
that the muMATB files through ARITBJfUS have been loaded. 

For positive integer Nr the "postfix* factorial operator named *!” 
returns the predict of the first N successive integerSr and 0! returns 
1. For example, 3i yields 6, which is 1*2*3. Use this operator to 
determine the predict of the first 100 integers: % 

BDS: FALSE $ 

% The nus^ser base used for input and output is initially ten, but the 
RADIX function can be used to change it to any base from two through 
thirty^ix. For example, to see what thirty looks like in base two: % 

mm: 30 S RftDDC (2) ; THim ; 

% As you can see, the radix function returns the previous base, which 
is, of course, displa^^ in the new number base. This inforation helps 
to get back to a previous base. In base two, eight is written as 1000, 
so to see what thirty looks like in base eight: % 

RADIX (1000) ; THim ; 

% In base eight, sixteen is written as 20, so to see what thirty 
looks like in base sixteen: % 

RADIX (20) ; THim ; 

% As you can see, the letters A, B, ... are used to represent the 
digits ten, eleven, ... for bases exceeding ten. Now can you guess why 
we limit the base to thir^ six? 

In input ^ressiens, integers beginning with a letter as the most 
significant digit must begin with a leading zero so as not to be 
interpreted as a name. . For example, in base sixteen, tei is the letter* 
digit A, so to return to base ten: % 

RADIX (OA) ; 

% Why don't you now see what ninety-nine raised to the ninety-nine 
power looks like in base two and in base thirty-six, then return to base 
ten: % RDS: FALSE $ 

% As you may have discovered, it is ea^ to become confused and have 
a hard time retumii^ to base ten. Two is represented as 2 in any base 
exceeding 1, so a fooJ^roof wiy to get from ary base to any other is to 
first get to base two, then express the desired new base in base two. 
For example: % 

RADIX (2) ; RADIX (1010) ; 

% Now we are guaranteeably in base ten, no matter bow badly you got 
lost. 

New osnsider irrational arithmetic: Did you know that 
(5 + 2*6*(l/2))*(l/2) - 2‘(V2) - (3/2)*(l/2) 
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can be simplified to Or provided we make certain reasonable choices of 
brauK:hes for the square roots? In generalr simplification of arithmetic 
expressions containing fractional powers is quite difficult, but mul'iAlB 
makes a valiant attempt. For example: % 

4 * ( 1 / 2 ) ; 12 * ( 1 / 2 ) ; 1000 * ( 1 / 2 ) ; 

% Try siji^lifying the square roots of increasingly large integers to 
gain a feel for how the computation time increases with the complexity 
of the input and answer: % RDS: FALSE $ 

% An input cf the form (m/n) '"(p/q) is treated in the usual manner 
as (m*(l/<3))*P / (0^(1/^)‘“p . For example: % 

(V9) * (3/2) ; 

% For geometrically similar people, surface area increases as the 2/3 
power of the mass. Veronica wears a 1 square-meter bikini, and she is 
50,653 grsBS, whereas her look-alike iwsther is 132,651 grams. Use muHAIB 
to determim the area of her mother's similar bikini: % BDS: FALSE $ 

% 4*(1/2) could simplify to either -2 or +2, but muIiATH picks the 

positive real branch if one exists. Otherwise, muMATH picks the 
negative real branch if one exists, as illustrated ky the example: % 

(-8) * (V3) ; 

% What if no real branch exits? Then muHATH uses the unbound 
variable named #I to represent the IMAGINARY number (-l)''(l/2), and 
expresses the answer in terms of tl, using the branch having smallest 
positive argument. For example: % 

(-4) * (V2) ; ^ 

% Decent simplification of expressions containing imaginary ranobers, 
as described in lesson CLES41ALG, requires that file ALGEBRAJ^ be 
loaded. Meanwhile if you believe in imaginary numbers and you can't 
CCTitain your curiosity, wiy don't you experiment with them to see what 
muMATB knows about them: % EDS: FALSE S 

% As with maiwial ccm^tation, picking a branch of a multiply-branched 
functicn is hazardous, so answers therety obtained should be verified by 
substitution into the origineil problem or by physical reasoning. For 
this reason, there is a COriHOL VARIABLE named FBRCH, initially TRUE, 
which suppresses Picking a BRan(Z if FALSE. For example: % 

PBBCB: FALSE $ 4 * (1/2); 

% Users hav^ a conservative temperament might prefer to do most of 
their computation with FBRCB FALSE. 

This brings us to end of CXES2.ARI. Thou^ arithmetic, some of 
the features illustrated in this lesson may be foreign to you, because 
sometimes they are taught during algebra rather than before. Thus, if 
you have any algebra background whatsoever, we urge you to proceed to 
lesson CLES3.AIG evei if some of CLES2.ARI was intimidatii^. Naturally, 
as implied by its type, file CLES3*ALG requires a muldATH system 
containing files through ALGEBRAJtfO. 

If you decide not to proceed to algebra, but want to learn how to 
program using muSII4P, then proceed to li^sw PLESUmi. % 

BC20: #ECS0$ PBECH: TRUE? RDS () $ 
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%File: CLES3.ALG (c) 


10/30/79 


The Soft Warehouse % 


LINELEN3IH (78)?. tECaO: BCaO$ ECHO: FALSE? 

HJMNQH: DQ^VM: 6? DOCC!: 2? NUMDEN: IWHEXH): 0? FBBCE: TRUE? 

X: 'X? ECSO: TRUE? 

% *01X8 file is the third of a sequence of interactive lessons about 

inuHATB-79. This lesson presumes that the muMATH files through 
ALG£BRA«ARI have been loaded and that the user has studied the 
arithmetic lessees CLESl,iARI and a:<£S2.ARI. 

kn. VARIABLE is one to whic h no value has been assigned. 

Hath^aticians call such variables QDEITRMIRATES. You ma^ have already 
inadvertently discovered that if you use an unbound variable in an 
expressionr muMATH treats the variable as a legitimate algebraic 
unJ^wn. Moreover, muMATH attempts to simplify expressions containing 
such unbound variables by collecting similar terms and similau: factors, 
etc. For example: % 

2*X - X"2/X i 

% See if muMAIE automatically simplifies the expressions 
0+Y, Y+0, 0*Y, Y*0, 1*Y, Y*l, Y"l, 1*Y, and 2*(X+Y) - 2*X. % 

ICS: FALSE ? 

% Sometimes it is desirable to change a bound variable back to 
unbound status. This can be done by using the single-quote prefix 
operator, ', which looks like an apostrophe on many terminals. For 
example: % 

EG: X + 5; EG: 'EG; EG + 2; 

% Try assigning the value M*C''2 to E, then change E back to 
unbound status: % RDS: FALSE ? 

% You ma^ have xxsticed that some of the more draistic transformations, 
such as expanding products or integer powers of sums, are not automatic. 
The reason is that such transformations are not always advantageous. 
They may make an expression much larger and less comprehensible. 
However, they may be necessary in order to permit cancellations which 
make an e^qpressicn smaller and more con^reh&isible. Accordingly, there 
are a few control variables whose values specify whether or not such 
trsuisformations are performed. For example, the variable controlling 
expansion of integer powers of sums is called Ft'JREXID. This variable is 
conservatively initialized to zero, so that Integer powers of sums eure 
not expanded. For example: % 

EG; (X+1) *2 - (X*2-2*X-1) ; 

% Clearly this is an instance where expansion is desirable. When 
FWREXFD is a positive integer multiple of 2, then positive integer 
powers of sums are expanded, so let's try it: % 

IWREXFD: 2 ? EG; 

% Nothing happened! 

The reason is that muMATH does not automatically reevaluate 
previously evaluated expressions just because we change a control value. 
Not only would this be rather time consuming, but the ability to form 
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aqpressions from other expressions constructed under different control 
settings provides a valuable flexibility for constructing partially 
es^^anded expressions. 

Cn the other hand^ it is often desirable to reevaluate repressions 
under the influence of new control settings, and the built-in EVAL 
function enables this: % ^ 

EVAL (EG) ; 

% New that EWEEXPD is 2, see how {X+Y)*2 - {X-Y)'‘2 simplifies: % 

BDS: FALSE $ 

% In IBUNA3IS-79, dnominators are represaited internally as negative 
powers, and negative integer powers of sums are es^aanded if FWSSXFD is a 
positive integer multiple of 3. For example: % 

FWBEXPDs 3 $ 1 / (X^^l) ^2 # 

% What happens if 1 / ({X+l)*2 - X)'‘2 is evaluated under the influ¬ 
ence of FWREXPD being 3? For a little surprise, try it.% KDS: FALSE $ < 
% Even though (X+l)‘‘2 is WITHIN a negative power, it is itself a 
positive power, so how about trying again with PWHEXH) being 2*3: % 

BDS: FALSE $ 

% NOW, we would like to suggest a little experiment for you: The 
size limitation on algebraic expressions depends primarily upon the 
amount of unemployed mmory available for storing names, numbers, and 
program or algebraiic structure. Memory for the structural use is 
measured in units called NCDES, which happen to correspond to 4 bytes in 
muSIHF-79 on microcomputers. Node-space tends to be the limiting 
resource for algebraic expressions, and we can always determine the 
number of unemployed nodes t^ing the command: % 

RECLAIM {); 

% lAimbers and nodes which are no longer a part of any value that we 
can retrieve are automatically recycled intermittently, but the RECLAIM 
function forces this "g 2 uddage collection” process. The collection takes 
on the order of a seccxid, dependijig on memory size and processor ^eed; 
and these slight pauses are sometimes noticable in the middle of a 
printout or a trivial computation. On a computer with front panel 
lights, the coUectiens are also usually recognizable by the change in 
light patterns. 

Naturally, if we load an extravagant number of muMATH files into a 
single muMATH dialogue or if we save a number of relatively large 
expressions as the values of variables, then there will be relatively 
little unemployed space for our next computation. Not only does this 
limit the size of the next expression, but computation time also 
increases dramatically as ^ce becomes scarce, because relatively more 
time becomes devoted to increasingly fregumit collections. The moral of 
the story is: don't unnecessarily load too mai^ muMAini files or retain 
numerous expressions as the values of variables. 

Now, for the esqperlment: In order to gain an e^reciation for how 
computation time depends on the size of the input expression, answer, 
and unemployed storage, try timing each computation in the following 
sequence, until it appears that your space or patience is nearly 
exhausted: 
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BG:(1+X)*2; REXXAIMO; BG:rX3"2; REXIADK); BG:EX3"2; ... % 

JDS: FALSE $ 

% These polynomials are called "daise”, because there are no missing 
terms less than the maximum degree in each unbound variable. In 
contrast "sparse” polynomials are missing a large percentage of the 
possible terms less than the maximum degrees. If you are still in an 
experimeital mood, you may wish to try the following analogous sequence 
which produces extremely ^arse results: 

RECLAIMO ; (A+Br2; RECLAIMO: (A+B+C)"2j RECLAIM ();... % 

RDS: FALSE $ 

% Distributim of sums over sums is another transformation which can 
dramatically increase expression size but is sometimes necessary to 
permit cancellations. For example, this transformation is clearly 
desiradsle in the e:qpressicxi: % 

EG: X'"! - 1 - (X+1)*{X-1) ; 

% When the control variable named NUMMUM is a positive integer 
multiple of 2, then integers in MUMerators are distributed over sums in 
MQMerators. Similarly whei the variable is a positive integer multiple 
of 3, then monomials in numerators are distributed over sums in 
numerators, whereas when the variable is a positive integer multiple of 
5, th^ sjrns in numerators are distributed over mms in numerators. 

The reason for using the successive primes 2, 3, and 5, is that it 
provides a convenient way to independently control the three types of 
distributioi using one easily remembered coitrol variable name. 

The initiail value of HDMMUM is 6, because numeric and monomial 
distribution are recoverable .(as we shall see), because neither 
distributicn dramatically increases expression size, and because a lack 
of these distributions often prevents annoyingly obvious cancellations. 
For instance the expression 2*(X+1) - 2*X will not simplify uixless 
NUMMUM is a positive multiple of 2. Similarly XH • (XH) will not 
simplify to 0, since the egression is r^resented internally as 
XH -l*(X4l), which requires the -1 to be distributed over the man. 

Thus, to return to our exanple, % 

EG; 5 * MOMMUM; EVAL(EG} ; 

% To witness the great varied of possible expansions, we set % 
«M3UM: 0 $ EG: 4 * X'^ * (1+X) * (1-X) ; 

% Mow, successively EVAL EG with NUMMUM being 2, 3, 5, 6, 10, 15, and 
30: % RDS: FALSE $ 

% In interpreting these results, it is important to recall that 
negations are represented internally as a product with the integer 
coefficient >1, so NUMIXJM must be a positive multiple of 2 to distribute 
negations over sums. 

If positive values of MUMNUM cause expansion in ixmierators, how do 
we request factoring in numerators? 

Negative values of MUMNUM cause factoring of numerators. Moreover, 
tkm specific negative values csu^se factoring of the type which reverses 
the corresponding expansion. For example: % 
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X: 'X $ Y: 'Y $ NDMNUM: -2 $ BS: 10*X*2*Y + 15*X*3; 

MUMNUM: 3*NUMNUH; EVAl.(EG} ; 

% What about negative fflult4>les o£ 5? Sony folks, that's hard for 
computers as well as humans. However, we are working on it for future 
releases. Meanwhile, try out our semifactoring on the example 
3*X*Y*3/7 - 15*X*Y*V14 + 9*X*4*Y^2/7 % RDS: FALSE $ 

% As you may have guessed, there is a flag named DEMDEN which 
controls expansion and factoring among negative powers in a manner 
entirely analogous to NOMNIJM. Use it together with NUI4MCJM to expand the 
denominator then sa&ifactor the denominator of the expression 
X*2/{{X-Y)*(X+Y) + ^^2 + X*2*Y) % BDS; FALSE $ 

% You may have wondered why we chose the names IKJMMUM and DEMDEN. 
The reason is that there is another closely related control variable 
named OEMMUM, which controls the distribution of various kinds of 
denominator factors over the terms of corresponding numerator factors: 

A positive multiple of 2 causes integer denominator factors to 
be distributed; a positive multiple of 3 causes monomial 
factors to be distributed; and a positive multiple of 5 causes 
sum factors to be distributed. For example: % 

Y: 'Y $ DENDEN: NCMNUM: 0 $ EG: (5 + 3*X*2) *(Y+1)/(15*X*(44X)); 

DENNUM: 2 $ EVAL{BG) ; 

DENNDI4: 3*DENNUr4; EVAL(BG); 

DQJNUM: 5*DQm4; EWAL(B3) ; 

% Positive setting of DENNUM and MHMMJM are particularly useful for 
work with truncated series or partial fracticn expansions. For example, 
see if you can put the expression (6 + 6*X + 3*X''2 + X*3)/6 into a more 
attractive form: % EDS: FAL^ '$ 

% What about negative valuW of DE21EUM? 

A little reflection confirms that forming a common denominator 
reverses the effect of distributing a denomimtor. T!»js, expressions 
are put over a common integer denominator when DENMDM is a negative 
integer multiple of 2, expressions are put over a common monomial 
denominator when OENWUM is a negative integer multiple of 3, and 
expressions are put over a common sum denominator when DENNUM is a 
negative integer multiple of 5. For example: % 

X: 'X S DEjaOM: DQDEN: 0 $ EG: 1 + X/3 + (1+X)/X + (1-X)/(1+X); 
DENNUIi: -2 $ HS: EVAL(BG); 

DENNUM: 3*DQ4MJM; EG: £VAL(EG}; 

DENtEM: 5*DEMN:ii; EG: EVAL(BG); 

% Try fully simplifying the ecpressioa X*4/(X'^+X*2) + V(X+1) - 1 
by expanding over a commoi denominator, then fi^toriz^g: % EDS: FALSE S 
% As with NDMMUM and OEMDEH, the initial setting of DENNUM is 6, 
whidi causes distribution of numeric arui monomial de^mijiator factors 
over numerator sums. This tends to give attractive results for 
polynomials or series with rational-number coefficients, but the 
relatively costly common-denominator operation may be necessary for 
prdslems involving ratios of polynomials. 

You have now been exposed to the four most important algebraic 
control variables in muMATH. Together with EVAL, the various 
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combinations of settings of these variables give rather fine control 
over the form of algebraic expressions. muMAlE cannot read the user's 
mind, so it is important for the user to thoroughly master the use of 
these variables in order to achieve the desired e^^ects. 

Here are the most frequently useful combinations of settings for 
these three variables: 

PWRDCFD: 0; NUMMJM: DENDEN: DEiiNUH: 6} These initial values are 
usually good for gaieral-purpose work, when the user wants to view sme 
results before doing anything drastic or potentially guite time 
consuming. 

' FWFEXH): 6; NOMNDM: DE213E21: 30; DEXMJ14: -30; These settings yield 
a fully expanded numerator over a fully expand^ common denominator. 
This form gives the maximum chance for combination of similar terms. 
Moreover, a rational function equivalent to 0 is guaranteed to simplify 
to 0. However, valuable factoring information may be irrecoverably 
lost. 

PWHD3D: 0; NUMNtIM: DQ®EU: -6; DENNUM: •30; These settings yield 
a semifactored numerator over a semi-factored commcxi denominator. This 
form gives the maximum chance for cancellation of factors between a 
numerator and denominator. However, the factoring is done 
incrementally, term by term, so it may be necessary to first expand over 
a common denominator so that all cancellable terms have an opportunity 
to cancel before attempting factorization. 

FWHE3CH3: 2; NDMICTM: 30; DPOEM: -6; DEMNOM: -30; These settings 
ax& a good compromise between the advantages of expai^sion and factoring. 
Semi-factoring is done in the denominator where it is usually most 
important, but there is a maximum ofportunity for combination of similar 
terms in tt» numerator. 

PWBEXFD: 6; MUMMUM: DEMDEM: DEMMUM: 30; These settings are good 
for series expansions or partial fractions, because each term is fully 
eepanded over its own denominator. 

Again, we can't overen^hasize the importance of mastering the use of 
these fair control variables. Th^ are yoir primary tool for imposing 
your will on the simplificatioi process, and any lade of understanding 
of their proper use will ultimately lead to frustration. Accordingly, 
why don't you try the ^tove and various other combinations an examples 
of your own choosing, until the usage becomes second nadire: % 

FDS: FALSE $ 

% Congratulations on completing CLE53.ALG. If the mathematical level 
was uncomfortably high, proceed to lesson FLESIJ^. Otherwise proceed 
to CLES4.ALG. In either event, it is advisable to initiate a fresh 
muM^ eiviroment, because our experiments have altered control values 
and made assignments which could interfere with those lessons in 
nefauaous ways. % 

ECHO: tECHOS 
EDS 0 5 
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%File: CLES4.AD3 (c) 


10/30/79 


mae Soft warehouse % 


LnCLElCra (78)$ #£(20; sao? EmD: TFUES 

% This is the fourth of a segueice of nuMAHH calculator-oode lessons. 

There are some other algebraic control variables besides PWREXPDr 
MDHNUMr DENDENr and DENNUM; and they are occasionally crucial for 
achieving a desired effect. One of these, named NOHDEN, provides the 
logical completitKi of the latter three, by ccxitrolling the distribution 
of factors in numerators over the terms of denominator msos. NUMDEN is 
initially 0, but integer numerators are distributed over denominator 
sums when NDMD£2 is a positive integer multiple of 2, monomial 
numerators are distributed over denominator sums when NUMDEN is a 
positive integer multiple of 3, and numerator sums are distributes over 
denominator rnnns when N0MD£2I is a positive integer multiple of 5. For 
example: % 

NOMNUM: Df2DEN: DI20XJM: 0 $ NOMDEW: 30 $ 

X / (X*3 + X + 1) / (T + 1) ; EG: (X+Y) / (l+X+Y) / (Y+1) ; 

% Isn't that intriguing? It yields a sort of ■continued-fraction" 
representation. Now for the reverse direction, which performs a 
denesting of denominators which is less drastic than a single common 
denominator: % 

NDMDEN: -€ $ Z + 1 / (lA + 1/Y) / (l-Hf) ; 

% See if you can devise examples exhibiting dramatic simplifications 
arising f rcmi either directioii for this novel transfocmatiou Tte fact 
that it so naturally complements NDMNUM, DENDEN, and DENNUM suggests 
that it must be useful for something! % EDS: FALSE $ 

% Another control variable named BASEXF controls distribution of a 
BASe over terms in an ESCPonent which is a sum, or controls the reverse 
process which is collection of similar factors. As might be e^^ected, 
integer bases are distributed over exponent sums when BASEXP is a 
positive integer multiple of 2, monomial bases are distributed over 
exponent asss when BASEXP is a positive integer multiple of 3, and base 
sums are distributed'over exponent sums when BASEXP is a positive 
integer multiple of 5. Morever, the corresponding n^ative values cause 
collection of similar factors having the corresponding types of bases. 
BASEXP is initially -30. However, distribution (followed perhaps by 
collection) is sometimes necessary to let some of the terms in an 
exponent sum confine with the base. For example: % 

EG: 2"(2'tX) / 4 ; BASEXP: 2 ; EVAL (EG) ,* 

% See if ^u can devise an example which requires evaluating an 
e^^ression first with sufficiently positive BASEXP, then reevaluating 
with sufficiently negative BASEXP, or vice-versa: % BDS: FALSE $ 

% Another control variable named EXPBAS oxitrols the distributic»i of 
Exponents over BASes whicn are PIOXJCTS. Integer ex^nents are 
distributed over base products when EXPBAS is a positive integer 
multiple of 2, monomial exponents are distributed over base products 
whei SCPSAS is a positive integer multiple of 3, and expwent sums are 
distributed over base products when EXPBAS is a positive integer 
multiple of 5. Naturally, the corresponding negative multiples request 
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collection of bases which hanra sinilcu: exponents of the indicated type. 
Tte initial value is 30# and here are some examples where distribution 
permits net simplificatioi: % 

(X'‘(l/2) * Y) “ 2 ; (X*Y) *2 - X*2*y*2 ; (4*X"2*Y) * (1/2) ; 

% However# the user should beware that as with manual computation# 
distribution of noninteger acponents is not alw 2 ^s valid. Consequently# 
c^ervative irsers may prefer to gaierally operate with SOPBAS being 2. 
Moreover# distribution of es^onents tends to make expressions more bull^ 
whei no cancellations occur. For example % 

(X * Y * Z) " (V2) ; 

% In fact# there are instances where negative settings of EXPBAS are 
necessary to acheive a desired result. For example: % 

BG: 2^ * 3*X + (1+X)*(V2) * (1-X)*(V2) - (l-X*2)*(l/2) ; 

EXmSi -6 ; NUMNUH: 30 7 E^TAL (BG) ; 

% See if you can devise an example which requires evaluating an 
eo^ression first with sufficiently positive EXPBAS# then reevaluating 
with sufficiently negative EXPBAS# or vice-versa# in order to sis^lify 
acceptably: % 6DS: FALSE $ 

% The variable neuned PBRCH# already discussed in conjunction with 
fractional powers of numbers# also controls transformations of the form 
u*v*w —> u*(v*w). PBRCH is initially TRUE# but when PBRCH is FALSE# 
the transformation occurs only for integer w. Otherwise the 
transformation occurs for any w. The user shixild be aware that in some 
circumstances the selected branch is an inappropriate one# so that it 
may sometimes be necessary to set PBRC2 to FALSE. See if you can devise 
such an instance: % RDS: FALSE* $ 

% New# try the ttuasples O'^ and X'^O# to see what bagpensi % 

RDS: FALSE $ 

% The reason that 0*X is not automatically simplified to 0 is that 0"^ 
is undefined for nonpositive values of X# so the transformation could 
lead to invalid results. Of course# sometimes users know that the 
exponent is positive# or they are willing to assume it is positive and 
verify the result afterwards. Cmsequ^tly# there is a control variable 
named ZERCBAS# initially FALSE# which permits 
the transformation whan nonFALSE. 

Why then do we automatically simplify X'^O to 1 even though X could 
perhaps take on the value 0# giving the undefined form 0'‘0? Well# we 
also have a control vari^le for that# named ZERCEXP of course# but we 
initialized it to TRUE because: . 

1. If we are thinking of polynomials in X rather than any one 
roecific value of X# thoi we are fz&a to regard the polynomial 
X 0 as being formally equivaloit to i. 

2. One cannot do effective simplification of rational 
functions without this widely aerated transformation. 

3. Since 1 is the limit of X"0 as X approaches 0 from either 
side of the real axis# 1 is a reasonable interpretation even 
for 0*0. 
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NevertbelesSf there is room for disagreemait, and anyone who wishes 
is free to run with ZEBOE2C? FALSE. Why don't you try itr using some 
rational expression examples, in order to see how you feel about this 
issue? % BDS: FALSE $ 

% It is easy to forget the current control-variable settings, and it 
is even e^ to forget the existence of certain oxitrol-^ariables, so we 
have provided a handy-dandy function named FLAGS which retunis the empty 
name ”* after printing a of all the flags and their values: % 

FLAGS {)? 

% If ^ ever get frustrated because you can't get an answer close to 
the desired form, try this command. It may reveal some inapprc^riate 
settings or remind you of some alternatives you forgot, or reveal the 
existence of potentially relevant flags of which you were unaware. 

Often a dialogue proceeds best under some control settings which are 
suitable for the majority of the computations, with an occasioial need 
for an eimluatioi under different control settings. Ea^h mich reception 
could iivralve new assignments to several control variables, followed 
an evaluation then assignments to restore the variables to their usual 
values. This process can become tedious, and baffling effects can 
result from inadvertently forgetting to restore a coitrol variable to 
its usual value. Consequently, as a conveiience, we have provided some 
functions which, for the most commonly desired sets of "drastic" control 
values, establishes these values, reevaluates its argument, than allows 
the control variables to revert to their former values before returning 
the reevaluated argument. 

One of these functions is called EXPAND, because it requests full 
expansion with fully distributed denominators, bases, and exponents. 
More specifically, it uses the following settings: 

PWBEXTO: 6; NDMDE27: 0; DENDEN: DEINNUM: BASEXP: EXPBAS: 30; 

To see its effect, try EXPAND (((1+X)/(1-X))*2); % !OS: FALSE $ 

% Another one of these convenience functions is called EXPD, and it 
fully expands over a common denominator. Thus the internal control 
settings are the same as for EXPAND, except that DE13NDM: -30. Try 
EXPD (1/(X+1) + (X+l)*?) ; % RDS: FALSE $ 

% Finally, there is a convenience function named FCTR, and it semi¬ 
factors over a common denominator. It evaluates its argument under the 
following ccxitrol-variable settings: 

IXmMi DEIOEN: DENNUM: BASEXP: EXPBAS: -30; SWPSXTO: NC24DEN: 0; 

Since semi-factoring is done termwise, it may be necessary to use 
EXPD before applying FCTR to an expression, in order to get the desired 
result. See if yea can devise an instance where this is true: % 

9DS: FALSE S 

% This brings us to the end of lesson CLES4.ALG. The next lesson is 
CLESSJULG, but as before, it is advisable to start a fresh muMATH to 
avoid conflicts with bindings established in the current lesson. % 

BC20; #ECHO $ 

HDS 0 $ 


3 


%File: aiS5.A£iS (c) 


10/30/79 


Tiie Soft wardiouse % 


UHELEZGIS (78)$ tGCSO: ECBO$ BC20: m[E$ 

% It is often desired to extract parts of an expression. Particularly 
frequent is a need to extract the numerator or denominator of an 
expression. Accordinglyr there are built>in SELECTOR functions named 
NOM and DEN for this purpose: % 

DENNCMi 0 $ S3: (1+X) / X ; NOM {S3) ; DQI (BG) i 
(1 + BG); DEN (1 + BG); 

% As the last two examples illustrate, NUM and DEN do not force a 
common denominator or any other transformation before selactioi, so the 
denominator is always 1 when the expression is a sum or when the 
expression is a product having no n^ative powers. Try out HIM and DEN 
oi a few examples of your own to gain some experience: % EDS: FALSE $ 

% The Programming-mode lessons will explain how to completely 
dismantle an expression to get at ai^ desired part, such as a specific 
term, coefficient, base, or exponent. 

muMAXS r^resents the imaginary number (>1)‘'(V2) as fl, and muMATH 
does apprqpriate simplificati(xi of integer powers of tl. For example: % 

#I " 7 ; EXPAM) ((3 + #1) * (1 + 2*#I)) ; EXPAID ((X + #I*y) " 3) ; 

% Try it, you'll like iti % KDS; FALSE $ 

% The definition of the operator in file ALGEBRA.ARI also 
implements two higher^level transformations which we mention here only 
in passing: 

omMAIE r^resents the baise of the natural logarithms as #E and the 
ratio of the circumference to the diameter of a circle as #PI. Using 
these, muMAIH performs the simplification 

#E " (n * #I * #PI / 2) —> #I*n, 

where n is ar^ integer constant, after which the power of tl is reduced 
aj;pr^riately. Also, if a coitrol variable called TPGDCFD is a negative 
multiple of 7, then co^lex exponentials are converted to trigonometric 
equivalents. (The o^osite transformation for sines and cosines to 
complex exponentials for TRGEXPD « 7, is implemented by file 
TRGPOSAL3.) If your mathematical backgrouixi includes these facts, you 
might wish to experience them here. Otherwise you can safely ignore 
this digression: % BDS: FALSE $ 

% You may have wondered whether or not an assignment to a variable, 
say X, automatically updates the value of a bound variable, say EG, 
which was previously assigned an repression containing X. Let's see: % 

X: 5 $ Y: 'Y $ B3: X + Y ; X: 3 ; BG; EVAL (BG) ; 

% Apparently the answer is ”no”, at least if X is bound when the 
assignment to B3 is made. ^Qiis should not be surprising, because after 
contritxiting its value to the sEpression X + Y, all traces of the name 
X are absent from this expression. Eowever, suppose that we do a 
similar calculation wherein X is initially unbound: % 
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X: 'X $ BG: X + Y; Xi 3; EG; 

% As when we change control variables, previously ev^d.uated 
expressi<^ are not automatically reevaluated when we bind an unbound 
varible therein. However, we can always use EVAL to force such a 
reevaluation: % 

EVAL (EG) ; 

% Since we did not assign the result to BG, reevaluation of EG after a 
different assignment to X still has an effect: % 

X: 7 $ BG: EVAL (BG); 

% Siiws this time we did assign the result to BG, further changes to X 
can have no effect on BG regardless of evaluation: % 

X: 9 $ EG: EVAL (BG) ; 

% If these examples are not entirely clear, you had better take the 
time to experimentally leam the principles by trying some examples of 
your owm % EDS: FALSE $ 

% It is often desired to reevaluate an expression under the influence 
of a tempora^ local assignment to one of the variables therein without 
disturbix^ either the existing value of the variable or else its unbound 
status. The built-in EVSUB function provides a convenient method of 
accomplishing this effect. E7SUB returns a reevaluated copy of its 
first argument, wherein every instance of its second argument is 
replaced by its third argument. For example: % 

NUl-lNUM: 6 $ M: 'M $ C: 'C $ V: 'V $ BG: M*C*2 + VS*^V2 $ 

EVSOB (BG, M, 5); EVSDB (BG, M, KL-HQ); M; 

% Play around with EVSUB for awhile until you are absolutely sure that 
you understand the difference between substitution and assignm«it: % 
EDS: FALSE $ 

% You may have discovered that EVSUB also permits substitution for 
arbitrary subexpressicns as its second argument. For example: % 

M: ’M $ C: ‘C $ E: 'E $ EVSUB (M*C"2 +7, M*<r2, E) ; 

% To keep the algebra package small, we have not endowed EVSUB with 
aiy sophisticaticn about finding aigdsraically IMPLIGCT instances of its 
second argument in its first. See if you can find examples where EVSIB 
does not do a substitutim that you would like it to do: % EDS: FALSE $ 
% Here is an example where a desired substitution doesn't fully o^ur:% 

NUMNDM: 6 $ C: 'C $ S: 'S $ EVSUB (1 - 2*S*2 + S*4, S"2, 1 - CTl); 

% The reason we did not get the desired simplification to C''4 is that 
if the second argument is a power, it matches (Xily the same power in the 
first argument. We can usually circumvent such problems by instead 
using an equivaleit substitution wherein the second argument is a name 
rather than a power. For example: % 

IWEEXPD: 2 $ EVSUB (1 - 2*S"2 + S*4, S, (l-<r2) * (1/2)); 

% Here is a somewhat different example wherein a desired substitution 
does not occur: % 

EVSUB (2*C*S, C*S, a); 

% The reason is that if the second argument is a product, it matches 
only the same COMPLETE product in the first argument. Again, the remedy 
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is to use an equivalent substitution wherein the second argument is a 
name. For example: % 

EVSOB (2*G*S, C, C2/S)i 

% Eere is a final example for which a desired substitution does not 
occur; % 

EVSOB (C"2 + r2 - 1 + C + S, <r2 + S*2, 1); 

% Similarly to products, if the second argument is a sum, it matches 
only the same G^lELErE sum in the first argumoit. As before, we could 
circumvent tiMi difficulty by making an equivalent substitution of 
(1-C*2) * (1/2) for S, or (1-3*2) * (1/2) for C, but that would leave 
an ugly s<^are root in the answer. If our goal is to delete the 
m:bexpression C*2 + S^2 - 1, the: we can use to our advantage the fact 
that powers must match exactly for a substitution to take place: % 

EVSOB (0*2 + S*2 - 1 + C + S, C*2, 1 - S*2) ; 

% See now if you can use such techni^es to get your examples to work: 
% BDS: FALSE $ 

% This brings us to the end of the calculator-mode lessons. There 
are, of course, higher-level math packages in muMATH, but the fact is 
that from a usage stai^aoint, we have alreac^ covered the hardest part, 
which is understanding evaluati(xi, substitution, and the ramificati(X)S 
of the varia;ts algebraic oantrol variables. You will find that if you 
know the relevemt math, use of the higher-level packages is quite 
straightforward, easily learned from stuffing the 
correspcnding DOC files. 

We suggest that before commencing the Programming-mode lessons, you 
explore calculator-mode usage of the higher-level packages as far as 
your math background permits. Math curriculum sequences differ, but 
probably most users will be most comfortable trying the higher-level 
packaes in the approximate order EQN, SOLVE, ARRAY, MATRIX, LOG, 
IBGMBG, TRGPOS, DIF, INS and Since ^ce becomes increasingly 

scuce as higher-level packages are loaded, you may have to reread file 
READ1ST.TXT to leam bow to CCNDEMSE and SAVE if you haven't already. 

Now for some parting advice about getting the most out of computer 
symbolic math: 

First, storage and time consumptioi tends to grow dramatically with 
the number of variables in the input expressions, even if the ultimate 
result is fortuitously compact. For example, the number of terms in the 
expanded form of 

(XI + X2 + ... + XM) * N 

grows outrageously with M and N. Consequently, it is important to make 
every effort to avoid needlessly introducing extra variables for 
g&»rality*s sake. Mathematical and physical problems are often stated 
using more variables than are strictly necessary, so it is also 
important to elicit every qgportunity to reduce the number of variables 
from the original problem. Bere are some general techniques for doing 
this: 
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. If members a set of variables can be made to occur only 
together as instances of a certain subexpressic^r consider 
replacing the subexpression with a single variable. For 
example: 

a) If K, Xf and XO can be made to occur only as instances 
of the subexpression K*(X-X0) / then consider replacing 
this subexpression with a variable named perhaps RDX. 

b) Similarly, perhaps a axnbination such as could be 

replaced with E, or FEO*V^2/L could be replaced with BE. 

These are respectively instances of absorbing an offset 
together with a proportionality coefficient, renaming a 
physically-meaningful subexpression, and grouping 
quantities into dimensionless quantities. Host ^igineering 
and science libraries have books describing a more 
systematic technique called DIMENSXa^ AIC^YSIS, and an 
article in the Journal of Computational Flimsies (June 1977} 
eaplains how computer algebra can automate the process. 

2. Even when a variable cannot be eliminated, the complexity 
of expressions may be reduced if the variable can be made 
to occur only as instances of a subexpression. For 
exanple: 

a) If only even powers of a variable X occur, consider 
replacing X"'2 with a variable named perhaps XSQ. 

b) If X only occurs as instances of 2''X, 2''(2‘*X), 
2'‘(3*X),..., consider replacing 2'‘X with a variable named 
perhaps TWOTOX, yielding mere integer powers of that 
variable. 

Some other advice is to avoid fracticxial powers and denominators as 
much as possible. They don't simplify well, they consume a lot of 
space, and they tend to be hard to decipher when printed one* 
dimensicxxally. Often a change in variable can eliminate a fracti(»ial 
power or a denominator. 

Sometimes, even when a problem cannot be solved in its full 
geierali^, solving a few special cases enables one to infer a general 
solution which can perhaps then be verified by substitution or by 
induction. Alternatively, perhaps the original problem can be 
simplified by neglecting some lower*order contributions, in order to get 
an analytic solution which will at least convey some qualitative 
informaticn about the solution to the original problem. 

Sometimes only part of a prdslem or perhaps even none can be solved 
symbolically, and the rest must be solved numerically. If so, the 
attraapt at an analytic solution at least allows ^e to proceed with an 
approximate numerical solution having more confidence that a concise 
analytical soluti(^ has not be&i overlooked. % 

ECHO: iECHO $ BDS () $ 
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%Pile: PLESl.m (c) 11/01/79 The Soft Warehouse % 


LINELENGTE (78) $ tECBD: SCBO $ BCBO: FALSE $ 

#CQCENS£: CCMDQISE $ CCHSQISE: FALSE $ KFFIBST: 'RFFIBST $ 

mm (wim, tmm) $ 

FONcncN pioNr (sa) , 

mEN AT3! (m), ♦ERIOT (m) DOT 

♦FRIOT (LEAR), PRIMT (FIRST (ECV ) , #PRIOT (" . "), 

PRIOT (REST (1X1)) , #PRINr (? , SQ, 

amai $ 

mm (PRBGLINE, tPRlNTLINE) $ 

FDNCnCN PRDTELINE (Da) , 

PRIOT (1X1), NEWLINE (), DQ, 

ENDFUN $ 

ECHO: TRUE $ 

% This is the first of a sequence of interactive lessons ^30ut louSIMP 
programming. It presumes that you have read files READ1ST.TXT and 
LESSCNS.TXT, and executed at least one of the calculator mode lessons. 
It also presumes that you have loaded packages through TRACEJIUS. 

muSIMP supplies a number of built-in functions and orators. The 
calculator-mode lessens introduced a few of these, such as RDS, REHACN, 
■f, *, etc. These progamming-mode lessons introduce more built-in 
finictiais md operators, but more important, the less^ reveal how to 
supplement the iMiilt-ln functions and aerators with definitions of y^or 
own, thus ext&iding ffluSIMP itself. 

It is important to realize that, until the last programming-mode 
lessons, we will not deal with muMATB. Instead we deal first with 
muSIMP, the language in which muMATB is written. The illustrative 
examples for these first few lessons are utterly different from muMATO, 
because we want to suggest a few of the many other applications for 
which muSIMP is especially well suited, and because we want these 
lessons to be coropreh^ible regardless of math training level. 

Data is what programs operate upon. The most primitive UNSTRUCTURED 
aoiSIKP data are integers and names, collectively called ATOt4S to suggest 
their indivisibility by ordinary means. Some programs must distinguish 
these two types of atoms, so there are two corresponding RECOGNIZER 
functions: % 

INTEGER (X76#) ; 

NAME (X76#) ; 

EG: -3271 $ 

INTEGER (EG) ; 

NAME (EG) ; 

% Do you suppose that "137", " ", and "X + 3", with the quotation 
marks included, are integers, names, or invalid? Find oat for y»irself% 
RDS: FALSE $ 

% Data can be stored in the computer's memory. The location at which 
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a data item is stored is called its AEDKESS. An address is analogous to 
a street address <xi tte outside of a mailbox. The data stored there is 
analogous to mail inside the mailbox, as with a cow of mailboxes, the 
contents of computer memory can change over time. 

There are useful programs which deal only with unstructured data, 
but the most interesting applicatiois involve aggregates of primitive 
data elements. One way to make an aggregate of 2 data elements is to 
use a structural data element called a IKDE, which stores the addresses 
of the 2 data elenents. Thus, a node is "data" consisting of addresses 
of 2 other data items. 

For example, su^ose that we wish to represent the aggregate 
consisting of the name BILBO and his age 31. We could store the name 
BILBO beginning at location 7, the number 31 beginning at location 2, 
and the node beginning at location 4. Then, begining at location 4, 
there would be stored the addresses 7 eund 2, as illustrated in the 
following diagram: 

Address: 1 2 3 4 5 6 7 

+■ .——- ' " I--H-1-!■ —.< -' I .. 

Contents: 1 I 31 I 17 12 11 BUflO 1 

4—— 1 ... I— -H -!-—H-1-1— 

Is that clear? 

The specific placement of data within memory is managed auto¬ 
matically, so all we are cot^emed about is the specific name and mmber 
values and the coinectivity o| the aggregates. Thus, for our purposes 
it is best to sufpress the irrelevant distracting deta^ aissociated with 
the specific addresses. The following diagram is one helpful way to 
portray only what we are ccxicemed about: 

1 / 1 \ 1 

+-/—»—\-^- 

/ \ 

BUBO 31 

This imagery suggests the word "pointers" for the addresses stored 
in nodes. 

If you have seen one bisected box you have seen them all, so to 
reduce the clutter and thus emphasize the essential features, we 
henceforth represent such nodes by a mere vertex in our diagrams, giving 
schenatics such as 

A 
/ \ 

BILBO 31 

Although most muSI!4P programs use such aggregates internally, many 
muSIMP programs are designed to read and print data in whatever 
specialized notation is most suitable for the application. For ecample, 
muMAIH uses operator and functional notation. Sippose howe/er that we 
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want to specify such aggregates directly in iiput and out^t. How can 
we do it? I£ we have a nice graphics terminal, then then we 
conveniently could use diagrams such as the above. Host o£ us do not 
have nice graphics terminals, so we must use another external 
representation. For this purpose muSiMP uses a representation 
consisting of the first data item, followed by the second data item, 
separated by a dot and spaces, all enclosed in a pair of matching 
parentheses. For example: 

{BII£0.31) 

We call this representation of a node a DOITED PAIR, However, this 
is a diferent use of parentheses and periods from bow they are otherwise 
used in muSBlP input, so we must preceed the dotted pair by the single* 
quote prefix operator to indicate to the parser that we are using 
dotted-pair notation rather than the usual operator or functional 
notatioi: 

' (BUBO . 31) 

Moreover, we must use an ampersand rather than a semicolon as the 
expression-terminator in order to inform the driver to print the 
expression as a dotted pair rather than attempt to print it using 
operator and functional notation. We say "attempt” because not all 
dotted pairs are appropriate for operator or £uncti(xial printing, as we 
will explain in the last lessons. Here then is an example of dotted- 
pair input and printing: % 

'(78 . TRCMBCNES) & 

% Try a few of your own, and note what happens when you forget the 
single-g^ote or use a semicolon rather than an aspersand: % 

IDS: FALSE $ ’ 

% What about wh&i we want to repres^t an aggregate of more than two 
atomic data elements? For example, what if we want to include BXIBO's 
last name, BAGGINS? Well, we can let one of the pointers of a node 
point to another node, whose first pointer points to BILBO and whose 
other pointer points to BA3GINS. For exanple: 


/ A 
A 31 
/ \ 

BILBO BAOdENS 

We can ii^t this as a dotted pair nested within a dotted pair: % 
'((BILBO . BfiGGINS) . 31) & 

% Mote that we only cpote the outermost dotted pair. 

MOW suppose that we want to also include BUBO'S species, structured 
as follows: 
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A HCBBET 
/ \ 

A 31 

/ \ 

BUBO BAGGIN5 
How would you ii^t that? 

Remorber/ your ii^t nust be tendnated by an an^^ersand. 

% RDS: FALSE $ 

% we would iJipit it as: % 

EG: *(((BIIBO . BAOGIMS} . 31) . BCBBIT) & 

% An alternative structure for this information is the one 
correi^cnding to the uput 

’((BILBO . BAGGINS) . (31 . HOBBIT)). 

On a piece of scratch paper, sketch the correspcnding diagram, then hold 
it close to my face so I can check it. 


/ 0^0 \ 

\\-// 

\—/ % IDS: FALSE $ 

I ^es must be getting bad, I couldn't see it. Oh well... 

Since either element of a dotted pair can be a dotted pair, they can 
be used to represent arbitraty "binary tree structures". Moreover, 
although perhaps unprintable using pure dotted-pair notation, linked 
networks of binary nodes can be used to represent any data-structure 
vrtiatsoever. 

In order to do anything interesting with data aggregates, a pm^ram 
must be able to extract their parts. Accordingly, there are a pair of 
SELECTOR functions namd FIRST and REST which respectively return the 
left and right pointers'in a node. For example: % 

REST (EG) & 

FIRST (EG) & 

FIRST (FIRST (EG)) & 

REST (FIRST (EG)) & 

% See if you can extract BILBO and BAGGINS from l^, using nested 
compositions of FIRST aixVor R&ST: % RDS: FALSE $ 

% Cur answers are: % 

FIRST (FIRST (FIRST (EG))) & 

REST (FIRST (FIRST (EG))) 4 

% Deeply nested function invocations become difficult to type and 
read, so let's define our first muSIMP function named FFFIRST, so that 
FFFIRST (EG) could be used as shorthand for the first of the above two 
exanples and for ai^ analogous example thereafter: % 

FUNCTION FFFIRST (U), 

FIRST (FIRST (FIRST (U))) 

QDFUN 4 
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% If you are not using a hard-^opy terminal, jot down this function 
definition and all subsegueit ones for reference later in the l e s son . 

Debits the word fSI^TSI, the fun has just begun: Now that FFFIBST 
is defined, we can apply it at any subsequent time daring the dialogue. 

example: % 

EETIRST (BG) & 

fPPIRSr ('{((BIG . MAC) . CKTSCP) . (FREia . FRIES))) & 

% using the definition of EEFIRST as a model, define a function named 
RFFIRST which extracts the REST of the FIRST of the FIRST of its 
argument, then test RFFIRST on B3: % RDS: FALSE $ 
t Our solution is: % 

FWCnCN RFFIRST (POO), 

REST (FIRST (FIRST (POO))), 

QCFUN & 

RFFIRST (EG) & 

% The name FOO in the definition is called a PARAMETER, whereas EG 
where the function is applied is an example of an ARGUMENT. We can use 
any name for a parameter even a name which has been bosid to a value 
or even the same name as an argument. The name is merely used as a 
"dummy variable" to help indicate what to do to an argument when the 
function is subsequently applied. A function definition is like a 
recipe. It is filed away, without actually being EXECUTED until apElied 
to actual arguments. 

As another sin^e example, since an atcns is defined as being either 
a name or an integer, it is convenient to hatve a recognizer fur^tion for 
atoms, so that we do not have t6 test separately for names and atoms 
when we do not care which type of atom is involved. We could define 
this recognizer as follows: 

FUlCnCN ATOM (U), 

MAKE (U) GR NUMBER (U) 

QI3FUN & 

ActueU.ly, ATOM is already built-into muSiMP, but the example 
provi<^ a good c^portiviity to introduce the built-in infix CR curator, 
which returns FALSE if both of its operands are FALSE, returning TRUE 
otherwise. Try oit ATOM oi the examples -5, X and BG % RDS: FALSE $ 
% Analogous to OR, there is a built-in infix AND operator which 
returns FALSE if either operand is FALSE, returning TRUE otherwise. 
There is also a ixiilt-in prefix NOT operator which returns TRUE if its 
operand is FALSE, returning FALSE otherwise. Knowing this, see if you 
can define a recognizer named NCDE, which returns TRUE if its argument 
is a node, returning FALSE otherwise: % RDS: FALSE $ 

% In programming there is rarely, if ever, one unique solution, but 
ours is: % 

FUNCnCN NCDE (U), 

NOT ATCM (U) 

EtmjN & 

ICDE (BG) & 

NCDE (5) i 

% So much for trivial exercises. Now let's write a function which 
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counts the number of atoms in its eurgument. We will count each instance 
of each atOTif evei if some atoms occur more than cnce. 

At first this mav seem like a formidable task, because a tree can be 
arbitrarily branched! Bow can we anticipate ahead of time all of these 
possibilities. Well, let's procrastinate by disposing of the most 
trivial cases even though we can't yet see the whole solution: If the 
argum^ is an atom, ti^ there is exactly 1 at^ in it. 

So nmch for trivial cases. We havei't yet solved the whole problem, 
but it builds our self-confidence to make progress, so that is a good 
psychological reason for first disposing of the eas^ cases. Also, with 
the easy cases out of the way, we can turn cxor full intellectual powers 
on the harder cases, unfettered by any distractions to trivial loose 
ends. 

We are left with the case wmre we know we have a node. Perhaps we 
cmrld scmdaow subdivide the problOB into smaller cases? 

Let's see ... Nodes have a FIPST part and a REST part, so perhaps 
that provides the natural subdivision. Bmmm 

If we knew the number of atoms for the left part azxi ti^ number for 
the right part, clearly the number for the whole aggregate is merely 
their sum. But txsw can we find out ^e ixrnber of atoms in these pari^? 
Why not RECURSIVELY use the very functi^ we are defining to determine 
these two contributioisl 

It may sound like cheating to refer to the function we are defining 
from with the definition itself, but remembering that the definition is 
not actually APPLIED until sometime after its definitio: is complete, 
perhaps it will work. We are working in a highly interactive 
oivironment, so the quickest way to resolve questions adXHit muSIMP is to 
try it and see! Here then is a formal muSIHP function definition 
correspondii^ the the above informal English "algorithm”: % 

FUNCnCN tAICKS (U), 

WHEN ATOH (U), 1 EXIT, 

tATGMS (FIRST(U)) tATCNS (RESr(U)} 

EtOFUN 6 

% Bert we introduce 2 new con^^ts: The BODY ai a function definitio: 
can consist of a sequence of one or more expressions separated by 
commas. A CCNDmCNAL-EXIT is an expression coxsisting of a sequoxre of 
o:e or more e^qpressions nested between the matching pair of words WEEN 
and EXIT. Who: a &nction def inition is APPLIED, the expressions in its 
bo^ are evaluated sequoitially, until perhaps a conditioud exit cuxises 
an exit from the procedure or until the delimiter named ENDFUN is 
reached. For a conditional exit, the first expression after the word 
WHEN is evaluated. If the value is FALSE, tho: evaluatio: proceeds to 
the point immediately following the matching delimiter named EXIT. 
Otherwise, evaluation proceeds sequentially through the remaining 
egressions in the conditional exit, if any, exactly as if the body of 
the conditional exit replaced that of the function. The value of a 
ccxiditional exit is that of the last expression evaluated therein, and 
the value returned by a function is that of the last expression 
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evaluated therein when the function is ^iplied. 

Thus, #A!rOMS ij&nediately returns the value 1 whenever the argument 
is an atom, and otherwise the function breaks the problems into two 
parts which are necessarily smaller, hence closer to being atoms. Let's 
test it, starting with trivial cases first: % 

♦ATOMS (fOC) & 

♦AKMS (5) & 

SS 6 

♦ATOMS (ES) k 

% It looks promising, 1:^ it is still perhaps mysterious how muSIMP 
and ♦ATOMS keep track of all of these recursive function invocations. 
Since the trace package is aipposedly loaded, to trace the execution of 
♦ATOMS, we merely issue the command: % 

TRACS UATOMS) & 

% Mow every time ♦ATOMS is entered, it prints its name and argument 
values, whereas every time it is exit^, it prints its name followed by 
an «iual si^, followed by the returned value. Moreover, the trace is 
indented in a manner which allows correspcnding entries and exits to be 
visually associated. Watch: % 

♦ATOMS (POO) & 

EG & 

♦ATOMS (EG) & 

% Try a few examples of your own, until these new ideas begin to gel: 
% RDS: FALSE ? 

UMERAGE (♦ATOMS) & 

♦ATOMS (POO) & 

% Here is a function which counts cnly the number of integers in its 
argument: % 

EONCnCN ♦INEEGESS (U), 

WHEN INTEGER (U), 1 EXIT, 

WHEN NAME (U), 0 EXIT, 

♦INTEGERS (PIRST(U)) +- ♦INTEGERS (REST(U)) 

ElOFUN ? 

EG & 

♦INTEGERS (EG) ; 

% Now, using it as a model, try writing a function named ♦NAMES, which 
returns the number of names in its argument. If your first 
syntactically accepted attempt fails any test, try using TRACE to reveal 
tl» reascxi why: % RDS: FALSE $ 

% Our solution is ... 

On second ti»u^t, we won't give you our solutiou Consequently, if 
you were lary and didn't try, you had better try now, because the 
examples will get steadily harder now. % RDS: FALSE $ 

% The HEIGHT of an atom is 1, and the of a node is 1 more than 

the maximum of the two heights of its FIRST and REST parts. 
Accordingly, let's first write a function named MAX, which returns the 
maximum of its two integer arguments. There is a built-in infix integer 
comparator named ">", so here is a hint: 
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PUNCTICN MRX (im/IOT2), 

WHiU INH > IME2r ... BtlT, 

$ 

Enter such a definiticn, with e^rqpriate substitutions for the missing 
portions, then test your function to make sure it works correctly: % 
R3S: FALSE $ 

% Now, with the help of our friend HAX, see if you can write a 
function named HEIGHT, which returns the height of its argument: % 

FDS: FALSE $ 

% Our solution is: % 

FONCnCN HEIGHT (0), 

WHEN A2DM (U) , 1 EXIT, 

1 + MAX (HEIGHT(FIRST(U)) , HEIGHT(KEST(U))) 

EtmiN $ 

% This brings us to the end of the first prograinmijig~mode lessons. It 
may be a good idea to review this lesson before proceeding to lesson 
PLES2.TRA. % 

ECHO: #ECH0 $ 

MOVD (tPRINr, PRINT) $ 

MOVD (#PRINELINE, PRINILINE) $ 

CONDENSE: #CQDENSE $ 

RDS 0 $ 
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% File: 7LES2,TBk (c) 


11/01/79 


TSae Soft Warehouse % 


LIMELE2CIB (78) $ tECSO: ECHO $ ECHO: FALSE $ 

#CCme]S£: CQCQISE S COCO 

MCWD (ERIOT, #PRINr) ? 

FaNGna,’ print (sxi) , 

WJiU MOM (ESQ), ♦PRINT (BQ) BtlT, 

♦PRINT (LPAR)r PRINT (FIRST(ESa)), ♦PRINT C . "J, 

PRINT (REST(Ba)), ♦PRINT (REAR), 1X1, 

QDFUN $ 

MCWD (PRINTLINE, ♦PRINUJNE) $ 

FUNCnCN PRINTLINE (EXl), 

PRINT (DQ), mwE 0, m, 

eCFUN ? 

ECHO : TRUE $ 

% This is the secoid of a segu^ice of nuSH'lP progranming lessons. 

EQ is a prinitive nuSIMP Comparator function which returns TRUE if 
its two arguments are the same address or equal integers, returning 
FALSE otherwise: % 

FIVE; 5 $ BQ (5, FIVE) ; 

% Names are stored unic^ely, so two oixurences of a name must involve 
the same address: % • 

ACTCR: 'BCX3AET ; BQ (ACTOR, ’BCGAKT) ; 

% Here is an example of two different references to the same physical 
node: % 


DATE: '(JULT . 4) S FOO: DATE $ BQ (FOO, DATE) ,* 

% However, watch this: % 

BQ (DME, '(JOLT . 4)) }. 

% What happened? The two aggregates are DUPLICATES, but since they 
were independently formed they do not start with the same node. In 
fact, only the name JUIX is shared among them, as shown below: 


second 

DATE argvment 

A A 
/ \/ \ 

I \ \ 

I / \ \ 

jms 4 4 


Clearly it is desirable to have a more comprehensive equality 
con^parator whidi also returns TRUE for aggregates which are duplicates 
in the sense of priixting similarly. Let's write such a function, ceilled 
DUP. Following the geieral advice given in PLESl, let's first dispose 
of the trivial cases: 
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If either argument is an atom, then they are duplicates if and only 
if they are EQ. 

Otherwise, th^ are both nodes, which is the nontrivial case. Now, 
let's employ our "divide-and-conguer" strategem, using FIBST and REST as 
the partitioning. Two nodes refer to duplicate aggregates if and only 
if the FIRST parts are duplicates and the REST parts are duplicates. 
Moreover, that can be tested with our beloved recursion, using DUP 


See if you can write a corre^^ponding function named OOP: % 

RDS: FALSE $ . 

% There are many possible variants, but here is one of the most 
compact: % 

FUNCnCN OOP (0, V), 

WHEN ATOM (U), BQ (U, V) EXIT, 

WHEN ATOM (V), FALSE EXIT, 

WHEN OOP (FIRST(U) , FIRST(V)), DOP (REST(O), REST(V)) EXIT, 

EH)FUN $ 

% An interesting challenge for your spare time is to see how many 
different but reasonable ways this function can be written. 

Actually, there already is a built-in infix operator named 
which is equivident to DOP: % 

DATE: '(JULE . 4) $ 

DATE - '(JDLSf . 4) ; 

% Do you feel DUl^ to learn that an exercise duplicated an existing 
facility? 

It is crucial to understand exactly what the existing faciliti^ do, 
and the best way to learn that is to understand how they work by 
creating them indeptod^tly. 

Here is a good exercise: See if you can write a comparator function 
named SAMESHAPE, which returns TRJE if its two arguments are similar in 
the sense of having nodes and atoms at similar places. For example, 
SAMESHAPE {'((KINGS . R(X)K) . 5), '((QUEENS . 3) . PAWN)) 
is TRUE: % H>S: FALSE $ 

% This is one of those instances where we will not give the answer. 

Now, using the infix operator named see if you can write a 
function named CC^AINS which returns TRUE if its first argument is a 
duplicate of its second argument or contains a duplicate of its second 
argument. For example, 

((JULY . 4) . (1931 . FRIDAY)) 

contains (1931 . FRIDAY). It is at least as hard as DUP, so take your 
time and don't give easily. % RDS: FALSE $ 

% Here is a h^er exercise: The two aggregates 

A A 

/ \ / \ 

CARBOJ A CARBCN A 
/ \ / \ 

SULFUR IBCN IRCN SULFUR 
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are ISOMEPS because th^ are either the same atom or at every level 
either the left branches are isomers and the right branches are isomers, 
or the left branch of me is an isomer of the right branch of the other 
and vice-versa. Write a corresponding comparator function named 
Ia^OS. (It's similar to DUP, with a twist.) % !a)S: FALSE $ 

% Our answer is: % 

FDNCnCN ISOMERS (0, V), 

WHIN ATOM (0), BQ (U, V) EXIT, 

WHIN ATCti (V), FALSE EXIT, 

ISOMERS (FIRST(O) , FIRST{V)) AND ISOMERS (REST(U), raST(V)) 

OR ISC14ERS (FIRST(U), RES?r(V) ) AND ISOMERS (RESr(O), FIEST(V)) 
QCFUN $ 

% Because of all the combinaticns which might have to be chedced, the 
execution time for this function can grow quite quickly with depth. Try 
tracing a few examples of moderate depth: % RDS: FALSE $ 

% So far our functions have merely dismantled or analyzed aggregates 
given to them as arguments. None of our examples have constructed new 
aggregates. The dot of course results in aggregates, but this occurs as 
the dot is read. Moreover, since the single quote necessuily 
preceeding an outermost dotted pair prevents evaluatioi, bound variables 
in a dotted pair contribute merely their names ratl:»r than their values. 
For example: % 

EG: 7 $ '(EG . 3) & 

% What we want is a function which evaluates its two arguments in the 
usual way, then returns a node, whose two pointers point to those values. 
There is such a function, named AEJOIN: % 

ADJOIN (EG, 3) & 

% A dotted pair within a function definition is a static entity, 
frozen at the time the functicxi is defined. In cmtrast, a reference to 
ADJOIN within a function definition is dynamic. The node creation is 
done afresh, with the current values of its arguments every time that 
part of the function is applied. As an example of the use of ADJOIN, 
let's write a functi(mi named SKELETON, which constructs a new tree which 
is structurally similar to its argument but has the name of laigth zero, 
wherever its argument has an atom. Thus, when printed, the new 
a^regate will display the skeletal structure of the aggregate without 
visually->discernable atoms. For example, 

SKELETON {'((HALLOWEEN . GHOSTS) . WITCHES)) & will yield (( . ) . ) 

CK, let's recite the litary: What canes first? 

trivial cases. 

So, if the argiment is an atom we return what? 

' im ‘ . 


otherwise we have a node, which is the most general case. However, 
nodes have a FIRST and a REST, so can we sonehow recurse, using SKELETON 
an these parts, then combine them? 
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Yes, as foUowsi % 


FONCnCN aCELETON (C), 

WHEN ASOd (U), "" ESCIT, 

ADJOIN (an-JElXM (FIRST(U)), SSELEICN (BEST(U)}) 

ENEra $ 

SKELETON (' ((NCO . GOO) . (GUY . SAN))} & 

% Easy. Yes? 

Now it is your turn. Write a functioi named TREHiEV, which produces 
a copy of its argument in which every left and right branch are 
interchanged at every level. For example, 

TKEEREV {'{{MOO . GOO) . (GUY . (PAN . CAKE)))) & 
should yisld 

(((CAKE . PAN) . GUY) . (GOO . MOO}) 

% RDS: FALSE $ 

% If you didn't get the following solution, you may groan when you see 
how easa it is: % 

FUNCnCN TREERE7 (U), 

WHEN AJXai (U), U EXIT, 

ADJOIN (TREERE7 (FEST(U)), TREEREV (FIIIST(U))) 

ENDFUN & 

TREEREV {'(("Isn't" . that) . ea^)) & 

% Here is a somewhat harder exercise: Write a function named SUBST, 
which returns a cofy of its first airgument wherein every instance of its 
second argument is replaced ky its third argument. For example, if 


PHRASE: 

'((CIHIS . (GOSH . DARN)) . CAR) . (IS . {{GOSH . DARN) . BAD))) $ 

then SUBST (FHRASE, ' (GOSH . DARN), ' (expletive . deleted)) yields 

({(THIS . (expletive . deleted)) . CAR) 

. (IS . ((expletive . deleted) . BAD))) % RDS: FALSE $ 

% That's all folks. 

The next lesson deals with a special form of tree called a list. 
Mary people find lists more to their liking, and perhaps you will too.% 
ECHO: FALSES 

mm (fPRIOT, PRIOT) $ MOVD (»RINILINE, PRIOTLINE) $ 

CC2I3ENSE: fCODENSE $ EOiO: fECBO $ RDS () $ 
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%FiIe: SLES3 (c) 


11/01/79 


Ibe Soft Warehouse % 



UNELENSTH (78) $ 

«E<20: ECHO $ tCCNDEtlSE: CGtDEIiSE $ CCUlQiSE: FALSE $ ECBO: TRUE $ 

% This is the third of a sequence of interactive lessons about muSIMP 
pcograotning. 

Oftenr it is most natural to represent a data aggregate as a 
sequence or LIST of items rather than as a general binary tree. For 
example, such a sequence is quite natural for the elements of a vector 
or of a set. We can represent such a sequence in terms of nodes by 
having all of the FIRST cells point to the data elenents, using the PEST 
cells to link the sequence together. The last linkage node can have a 
REST which is FALSE to indicate that there are no further linlu^e nodes: 

A 
/ \ 
iteml A 

/ \ 

item2 . 



e 

A 

/A 

itemN FALSE 

Viewed at a 45 degree rotation, this diagram is analogous to a 
clothes line with the successive data elaaents suspended from it, thus 
more clearly suggesting a sequence. The simple regularity of the 
structure permits correspondingly simple function definitions for 
processing such structures. Moreover, the linear structure suggests an 
external printed representation which is far more readable than dotted 
pairs. In response to an ampersand terminator, muSlMP prints the above 
aggregate in the more natural LIST notaticxi: 

(iteml, item2, iteoN} 

rather than the equivalent dot notation 

(iteml . (item2 . ... (itecAl . FALSE) ... }} 

Conversely, the reader accepts list notation as an alternative 
ii^t form to dot notation. Naturally, any of the items in a list can 
themselves be either lists or more general dotted pairs. The printer 
uses list notation as much as possible. Thus, a structure of the form 
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A 
/ \ 
itaol A 

/ \ 
iten2 . 


• 

A 

/A 

iceiSl atoD 

where "atom” Is not the atom PALSEr is printed in a adxed notation as 

(itesBlf ite&2r (itemN . atom)) 

Similarly, the reader can a^rqpriately read such mixed notation. 

Ihus the last iten in a list is implicitly dotted with FALSE, and a 
blank between two items is equivalent to " . (", toqether with a 
matchii^ adjacent to the next right parenthesis. You may wonder why 
you never noticed such printing conventions during lessons PLESl emd 
FLES2. Tte reason is that we purposely redefined the printer for those 
lessons so that it did not use the list-abbreviation convention. 

It is important to fully understand the connection between dotted 
pairs and lists, so take 5 minutes or so to type in some lists, nested 
lists, nested dotted pairs, and mixtures, noting carefully how they 
print. % BDS: FALSE $ 

% Did your exanples include: % > 

'{) & 

% Is that surprising? 

Since FALSE is used to signal the eiui of the list, FALSE and the 
empty list must be equivalent. 

OLearly the trivial terminal case in processing lists will involve 
an equality test against FALSE. Since this test is so common, there is 
a correspcnding built-in recognizer defined as follows: 

FUNcncN Sim (lis) , 

BQ (LIS, •()) 

E2EFUN; 

Using EMPTY, see if you can define a function named #IT£MS, which 
returns the number of (top-level) items in its list argument. For 
example, #ITEMS ('(FIOG, (FBDIT . BAT), NEWT)) should yield 3. Here is 
an incomplete solution. All you have to do is enter it with the 
portions marked ”..." apprcpriately filled. 

FUNCnCN #ITn^S (LIS), 

WHIN n4PTY (LIS) , ... EXIT, 

1 + #IIENS ( ... ) 

ENDFUN $ % BDS: FALSE $ 
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% Actually, tl^e is alreac^ a built-in function caUeci LfSlGTH, which 
raturns the length c£ a list. It is somewhat more general in that it 
returns the number of characters necessary for printing when giv^ an 
atom. 

liote that with lists it is typical to recur only on the REST of the 
list, whereas with geieral binary trees it is typical to recur on both 
the FIRST and the REST. 

So far, the examples and exercises have been relatively isolated 
ones. Now we will focus on writing a collection of functions which 
together provide a signifi«Bit a^licaticns package: 

A list provides a natural representation for a set. For example, 
(MANGO, (CHOCOLATE . FUDGE), (ALFALFA, SPROUTS}) can represent a set of 
three fo^s. Using this representation, let's write functions which 
test set membership and form unions, intersecti(xis, etc. 

First, write a function named ISIN, which returns TRUE if its first 
argument is in the list which is its second argument, returning FALSE 
otherwise: % RDS: FALSE $ 

% Our solution is: % 

FUNCnCN ISIN (U, LIS) , 

WHE21 WPn (LIS), FALSE BCIT, 

WHEN U - FIRST (LIS), EXIT, 

ISIN (U, REST(US)) 

ENDFUN $ 

ISIN ('FROG, '(SALAMAICER NEWT-TOAD)) ; 

% Actually, there is already a built-in version of ISIN called MEMBER. 

A set contains no di^licates, so we really should have a recognizer 
function named ISSET, which returns TRUE if its list argument ^stains 
no duplicates, returning FALSE otherwise. Try to write such a function: 
% RDS: FALSE $ 

% Here is a hint, in case you gave i^: 

FUNCnCK SET (LIS), 

WHEN ... EXIT, 

WHEN MEJBER (FIRST(LIS), ... ), FALSE EXIT, 

ENDFUN; % RDS: FALSE $ 

% In case it isn't clear by now, a rule of this game is that you are 
free (and encouraged) to use any functions we have already discussed, 
whether th«^ are built-in, previous examples, or previous exercises. 
That is one reason it is adviseable for you to actually do the 
exercises. 



NOW write a function named SUBSET, which returns 'H^E if the set 
which is its first argument is a subset of that which is its second 
argument. (Remcanber that every set is a subset of itself and the empty 
set is a subset of every set.) % RDS: FALSE $ 

% Here is a hint, in case you gave up or bad a less cccpact solution: 
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FUNCnCN SUBSET {SEU, SES2) , 

WHEN ... Bcrr, 

f«iQl member (FIPSTtSm), ...), SOBSET( ...) ESOT 
E2SSTJN; % BDS: FALSE $ 

% Two sets are equal if and only if they contain the same elements. 
However, the elements need not occur in the same order. Write a 
corre^nding comparator function named EQSET: % BDS: FALSE $ 

% Ah yes, a hint perhaps?: 

FONCnCN EDSST (SETl, SEE2), 

• • • 

EIGFUN; % BDS: FALSE $ 

% Do you think that's not much of a hir&? 

Well, tile bod^ of the function really can be written with one modest 
line, so t^ harder: % BDS: FALSE S 

% Remember the rules of the game: You are encouraged to use any 
function discussed previwjsly: % 

FONCnCN EDSET (SEU, SET2), 

SUBSET (SEn, SET2} AND SDBSEI (SEI2, SETl) 

ENDFUN; 

% Our ecamples so far have merely analyzed sets. We can use ADJOIN to 
construct lists, just as we used ADJOIN to construct binary trees. As 
an example of this, write a functioi named t<lAKESEr, which returns a copy 
of its list argument, except without duplicates if there are any: 
% IDS: FALSE S 

% If you need a hint, here is one, but it is all you will get: 

FONCnCN MAKESET (LIS) 

WHEN ..., '() HOT, 

IBEN ME2BER ( ... ), ... SaCIT, 

ADJOIN ( ... ) 

ENDFUN; % EDS: FALSE $ 

% Let's see if your solution works correctly: % 

MAKESET ('(FROG, FROG, FROG)) & 

% If there is a duplicate in the answer, then back to the computer 
terminal: % BDS: FALSE $ 

% (It helps to think of nasty test cases BEFORE you start 
prograianix^}. 

Now for the crowning glory of our set package: The UNION of two 
sets is defined as the set of ail elements which are in either (perhaps 
both) sets. Give it a try: % BDS: FALSE $ 

% A hint perhaps? Well, the fiaiction body can be written in 3 lines, 
each of which b^ins just like the corresponding line in our hint for 
MAKESET. % BDS: FALSE $ 
t Here is our solution: % 

FUNCnCN ONICN (SEU, SEr2), 

WHEN EMPTY (SEU), SET? EKIT, 

WHEN ME 2 BER (FIRST(SETl), SET2), UNICN (BEST(SEn), SET2) EXIT, 
ADJOIN (FIRST(SETl), UNION (REST(SETl), SET2 )) 

ENDFUN $ 
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% Cie Intersection of two sets is the set of all elements which are in 
both sets. Using our definition of UNION as inspirationr write a 
corresponding functi(» for the intersection: % BDS: FALSE $ 

% so farr our set algebra package has been developed in a so-called 
BOnOH-UF inaner, with the most primitive functions defined first, and 
with the more sophisticated functions defined in terms of them. The 
opposite ^roach is TOP-DOWN, where we define the most con^rehmisive 
functions in terms of more primitive ones, then we def ine those more 
primitive ones in terms of still more primitive ones, until no undefined 
functions remain. 

As an example of the top-down attitude, let’s write a SYMMETRIC 
DIFFERENCE function for our set-algebra package. The ^’mmetric 
difference of two sets is the set of all ele&ents which are in exactly 
one of the two sets. This is in contrast to the ordina^ dif erence of 
two sets, which is all of the elements that are in the first set but not 
the sec(x^ However, if an ordinary difference function was available, 
we could write the symmetric difference as the union of the ordinary 
difference between setl and set2, with the ordinary difference betwe«i 
set2 and setl. We have already written UNION, but an ordinary set 
difference is not yet available. Nevertheless, let's bravely proceed to 
write the symmetric difference in terms of the ordinary difference, then 
we will worry about how to write the latter: 

% 

FUNCnCN SXMDIF {SEU, SEI2 ), 

UNION (CREOIF (SETl, SEK), QRDDIF (SEK, SETl)) 

ENDFON $ * 

% Now you try to write OHDDIF. It may help you to know that it can be 
written very similarly to QHCN: % RDS: FALSE $ 

% Some programmers are initially uncomfortable with the top-down 
approach because it makes them nervous to refer to undefined functions: 
there are obvious loose ends during the writing process. However, it is 
not necessary to understand how an auxiliary function can be written 
before daring to refer to it. All that is necessary is that the duty 
relegated to the auxiliary function be stxnehow more elen^tary than the 
overall duty performed by the function which refers to it. 

There are necessarily loose ends during the writiiig of a program in 
any sequential order. With the bottom-up ^roach, the loose ends are 
neither written nor referred to until lower-level functions have been 
written. Unfortunately, as such hidden loose ends are revealed they 
often make apparent the need to completely reorganize and rewrite all 
aibordinate functions into a more suitable organization. In contrast, 
the obvious loose during a top-down development provide invaluable 
clues about how to organize the remaining functions. Moreover, any 
subsequent changes tand to be easier, because communication between the 
functions is more localized, more independent, and more hierarchial. 
For example, we know that in the definition of SYMDIF we are taking the 
union of two DISJOINT sets, because from the definiticn of CRDDIF it is 
clear that ORDDIF (SETl, SET2) and ORDDIF (SET2, SZTl) cannot have 
elements in common. Hence it would be more efficient merely to append 
the secawi ordinary set difference to the first ordinary set difference, 
or vice-versa. Unfortunately, ADJOIN does not accomplish the desired 
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effect. 


For example, ADJOIN {'(5, 9), '(3, 7)) yields ((5, 9), 3, 7) 
rather than the desired (5, 9, 3, 7). What we must do is ADJOIN 9 to 
(3, 7), then adjoin 5 to that result. See if you can generalize this 
process into a fuixtion naned AFPEND, which returns a list caisisting of 
the list which is its first argumait appended onto the beginning of tt^ 
list which is its second argument:% EDS: FALSE $ 

% BOW about: % 

FUNCnCN APPEND (LISl, LIS2), 

WHEN EMPn (LISl), LIS2 EXIT, 

ADJOIN (FIEST(LISl), APPEND (EEST{LIS1), LIS2)) 

QCFUN $ 

% You may not be getting tired, but my circuits are weary, so let's 
bring this lesson to a close, % 

£G!0: «ECEO $ CCNDQISE: tCCNDENSE $ 

EDS {) § 
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%Fiie: PLES4.'IRA (c) 11/02/79 The Soft Warehouse % 


LINELE2CIH (78) $■ 

fCODStiSE: CCCCEMSE $ Cai}E21SE: FALSE $ tECHO: SCSO $ EC20: TROE $ 

% This is the fourth is a series of ouSIMP pcogranining les8c»is. 

Often within a functioi definition it is desired to form a list of 
values OTNAMIG^LUf. For example, si^ipose that we wish to form a list of 
the VALUES of the variahles FIBSTNAME, LASTNAME, and MAILADDRESS. It 
Will not do to use '(FIBSTNAME, LASTNAME, MAILADDRESS), because the 
(^»3te prevents evaluation of the variables. 

We can accomplish the desired effect by writing 

ADJOIN (FIBSTNAME, ADJOIN (LASTNAME, ADJOIN (MAILADDRESS, '()))). 

However, this rather unreadable construct is tedious to write. 
Fortunately, muSIMP provides a convenient function named LIST for this 
purpose: We can accomplish the desired effect by merely writing 

LIST (FIBSTNAME, LASTNAME, MAILADDRESS). 

Qilike most functions, LIST uses any mmber of arguments. As specific 
examples: % 

FIBSTNAME: 'JCEN & 
lASTNAME: 'DCE & 

MAILADDRESS: 'TIMBUKTU & 

% New, ccopare usir^ a quote with using LIST: % EDS: FALSE $ 

% Reversing a list is an occasional ne^, and it is somewhat tricky to 
write a function for this. The following partial definition reveals 
that our friends APPEND and LIST can help: 

FDNCTICN REVUST (LIS), 

WHEN ... EXIT, 

APPEND ( ... , UST (FIRST (LIS))) 
a®FUN § 

See if you can successfully complete this definition. Naturally, you 
also have to reenter APPEND if a correct version is not around from the 
previous lesson. (Remember also to jot down all function definitions if 
you are not using a hard-cc^ terminal.) % RDS:FALS: $ 

% A well-written APPEND necessarily requires execution time which is 
approximately proportional to the length of its first argument. The 
REVLIS fundticxi outlined above invokes APPlZaD n times if n is the length 
of its original argument, and the average length of the argument to 
APPEND is V2. Thus, the time is approximately proportioial to n*(n/2), 
which is proportiomd to n^^l. 

Fortunately, an important technique called a COLLECTION VARIABLE 
permits list reversal in time proportional to n, yielding tremendous 
time savings for long lists: % 
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FUNCnCN REVLIS (LIS, ANS) , 

WHQJ EMPEf (LIS), ANS EXIT, 

REVLIS (REST(LIS), ADJOIN (FIRST(US), ANS}) 

QOFUN $ 

TRACE (REVLIS) ; 

REVLIS (’(1, 2 , 3)) & 

% A collection variable accumulates the answer during successive 
recursive invocatiois. Then, the resulting value is passed back through 
successive levels as the returned answer. 

As is illustrated here, we can invoke a function with fewer 
arguments than there are parameters. When this is done, the extra 
parameters are initialized to FALSE, and ti^ are available for use as 
LOCAL VARIABLES within the function body. Quite often, as in this 
example, the initial V2due of FALSE is exactly what we want, because it 
also represents the empty list. (Miai we want some other initial value, 
either the user can supply it, or the function can supply it to an 
auxiliary function which d<^s the recursion.) 

Of course, if a user of REVLIS supplies a second argument, then the 
function returns the reversed first argument appa:ided onto the secoxvi 
argument, which is eilso occasionally useful. 

What if the user supplies more arguments than there are parameters? 
The extra arguments are evaluated, but ignored. This is also 
occasionally convenient. 

The style of programming exemplified so far is the so-called 
"applicative" style polarized by the influential Turing lecture of J 
Backus, published in the August 1978 issue of the Communications of the 
ACM: The emphasis is on expressions, functional composition, and 
recursion. 

muSIMP also supports the alternative "Von Neumann" style emphasizing 
lo^)s, assignments, and other side-effects. To illustrate this style, 
here is an alterative definition of REVLIS which introduces the LOOP 
construct: % 

FUNCnCN REVLIS (LIS, ANS), 

LOOP 

WHEN EMPTf (LIS), ANS EXIT, 

ANS: ADJOIN (FIRST(LIS), ANS), 

LIS: REST (LIS) 

ENDLOQP 
E2DFDN $ 

% ffluSIHP has a primitively defined function named REVERSE, which has 
an equivalent machine language definition. 


An iterative loop is an expression consisting of the keyword LOOP, 
followed iiy a sequence of one or more expressiois separated by commas, 
followed by the matching delimiter named ENDLOOP. The bo(^ of a loop is 
evaluated similarly to a function body, exc^: 

1. When evaluation reaches the delimiter named ENDLOOP, 
evaluation proceeds back to the first e^q^ression in the loop. 
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2. When evaluation reaches an QOT within the locpr evaluation 
proceeds to the point imsediately following EII3L0QP/ and the 
value of the loop is that of the last expression evaluated 
therein. 

There can be any mmber of conditional exits emywhere in a loop. 
Ordinarily there is at least (»e exit unless the user plans to have the 
loop repeat indefinitely until perhaps interrupted by typing ESCape, 
ALTbode or CIHr’Z. interrupt can stxxreed only if the loop invokes 
at least one function which is not built-into muSIMF.) 

Now consider the following sequence: % 

U: '(THE OR33GDJAL ) $ 

L2: '(TAIL) $ 

REVLIS (LI, L2) & 

% The above definition of REVLIS makes assignments to its parameters 
LIS and ANS. For this example, the final assignments are LIS: '(} and 
ANS: '(ORIGINAL, THE, TAIL!. So, what do you guess are the 
corresponding current values for LI and L2? See for yourself: % 

ICS: FALSE $ 

% The assignments to parameters LIS and ANS have no effect on 
arguments LI and L2! This "call-by-value" mechanism permits function 
definitions to freely utilize their parameters without fear of dashing 
the values of user's argument variables outside. Thus, ordinary 
functioi parameters are never employed for passing information back to 
the user. If we wish to return mgre than one piece of informaticxi, the 
most well-disciplined way to. do so is to return an aggregate of the 
pieces as the returned value. However, another way is to make 
assignments within the function body to verifies which are not among 
its parameters — so-called "fluid" or "global" variables. 

As is often the case for iteratioi versus recursion in muSIHF, the 
iterative LOOP version of REVLIE is slightly faster than the recursive 
collection-variable version, but the latter is more con^ct. When there 
is such a trade-off betweei speed and compactness, a good strategy is to 
program for speed in the crucial few most-freguently invoked functions, 
and program for coi^ctness elsewhere. 

However, losing does have another advantage wt^ it is applicable: 
Recursion entails a "stack" of information which grows with the d^th of 
recursion. C^egueitly, evea though the allocated to the stack 
is quite generous, excessively deep recursion can abort a computation by 
exhausting this space. 

For practice with loops, use one to write a nonrecursive reosgnizer 
named ISSET, which returns TRHE if its list argument contains no 
(iplicate elements, returning FALSE otherwise. (C^pare your definition 
with the recursive version in lesscx: PLES3.) % RDS: FALSE $ 

% Here is our soluticxi: % 
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FUNCnCN ISSET (LIS) , 

LCX3P 

WHEN SJim (LIS), DCIT, 

WRQJ MEMBER (FIRST(LIS), REST(LIS)) , FALSE EXIT, 

LIS: REST (LIS) 

EH3L00P 
a®FUN $ 

% Another good exercise adapted from PLES3 is to use a lo<^ to write a 
nonrecursive function named SUBSET, which returns TRUE if its first 
argument is a subset of its second argument, returning FALSE otherwise: 
% RDS: FALSE $ 

% A BLOCS is another control construct which is sometimes convenient, 
particularly in conjuction with the Von Neumann style. As an 
illustration of its use, the following iterative version of the MAKESEI 
function from FLES3 returns a set composed of the unique elements in the 
list which is its first argument: % 

FUNCnCN MAKESET (LIS, ANS), 

LOOP 

WHQJ EMPTY (LIS), ANS EXIT, 

BLOCS 

WHEN MEMBER (FIPST(LIS), ANS), EXIT, 

ANS: ADJOIN (FIRST(US) , ANS) 

EIDBLCXX, 

LIS: REST (LIS) 

ENDLOOP 
ENDFUN $ 

MAKESET ('(FROG, FROG, FROG, TERMITE)) & 

% When evaluation reaches an BCIT, it proceeds to the point following 
the next ENDBLOCS, ENDLOCP, or ENDFUN delimiter — whichever is nearest. 
Thus, BLOCS provides a means for alternative evaluation paths which 
rejoin within the same function body or loop body, without causing an 
exit from that body. The first expression in a block must be a 
conditional-exit (arything else can be moved outside anyw 2 ^), but since 
there can be ana number of other oxiditicnial exits or other expressions 
within the block, the block provides a very general structured control 
mechanism. For example, the CASE-statem^t and IF-TBEN-ELSE construct 
of some other languages sure essentially special cases of a block. 

You may not have noticed, but the loop version of MAKESET has the 
effect of reversing the order of the set elements. Using ADJOIN in a 
loop generally has this effect, which is why it is so suitable for 
REVERSE. With sets, incidental list reversal is perhaps acceptable, but 
for most applications of lists it is not. We could of course use a 
preliminary or final invocation of REVERSE so that the final list would 
emerge in the original order, but that would relijxjuish the speed 
advemtage of the loop approach, while further increasing its greater 
bulk. Thus, recursion is usually preferable to loops when ADJOIN is 
involved. For example, recursion is used aULmost exclusively to 
implenent mu^]ATH, because its syndic &eprassior)S are r^resented as 
ordered lists. 

Loops are also less applicable to general tree structures than to 
lists, but it is often possible to loop on the REST pointer while 
recursing on the first pointer, or vice-versa, particularly if ADJOIN is 


4 


not involved. For exaapler compare the following semi-recursive 
de£initic»i of fATOMS with the fiiUy-recursive one in FLESl: % 

FDNCnCN #AaEMS (U, N) , 

N: 1, 

LOOT 

msn imDM (U), N DOT, 

N: N + lAICMS (riRSr(U)), 

0: REST (0) 

ENDLCOT 
aOFUN $ 

♦ATOMS ('((3 . FOO), BAZ)) ? 

% If the answer surprises you, don't forget the FALSE which BAZ is 
implicitly dotted with. 

See if you can similarly write a semi-recursive function named IXJP 
which does what the infix operator named does: % RDS: FALSE $ 
t Those of you with previous exposure to only Von Neumann style 
programming undoubtedly feel more at hone now. The reason we ^stponed 
revealir^ t±ese features until now is that we wanted to force the use of 
applicative programming long enough for you to appreciate it too. 
Naturally, one should employ whichever style is best suited for each 
application, so it is worthwhile to become equally conv^sant with both 
styles. 

Ttus endeth the sermon. % 

ECHO: ♦ECHO $ COCENSE: #CaCENSE $ 

RDS () $ 
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tFile: PLESS.m (c) 


11/01/79 


The Soft Warehouse % 


LINELESGTE (78) $ 

fCQCENSE: CQCQiSE 8 COtDBUSEi FALSE 8 fECEO: ECSO 8 ECHO: TRUE 8 
% This is the fifth in a sequence of nuSIMP progr amin g lessons. 

In the previous lesson our original version of REVEBSEr called 
ISlVLlSr required time proportional to n'^lf where n is the length of tte 
first argument. We then showed how a collection variable or a loop 
could yield a much faster technique using time proportional only to n. 
Kow, let's consider the i^eed of some of the other set functicris that we 
defined: 

Whether iterative or recursive, MEMBER can require a number of 
equality comparisons equal to the length of its second argument. 
Whether defined iteratively or recursively, SUBSET, EQSET, UNION, and 
INTERSECTIC^ all require a membership test for each element of one 
argument in the list which is the other argument. Thus, these 
definitions can all consume computation time which grows as the product 
of the lengths of the two arguments. By similar reasoning, the one- 
argument functions ISSET and MAKESET are seen to require time 
proportioial to the square of the leigth of their argumoit. Data-baise 
applications and others can involve thousands of set operations on sets 
having thousands of elements, so it is worthwhile to seek methods for 
which the computation time grows more slowly with set size. 

In muSIMP, every name has, an' associated PROPERTY LIST which is 
immediately accessible in an amount of time that is intepeident of the 
total number of names in use. Provided the elanents of the sets are all 
names, this permits techniques for the above set operations requiring 
time proportional merely to the length of the one set or to the sum of 
the lengths of the two sets. 

A property list is a list of dotted pairs. The first of each dotted 
pair is a atom called the KEY or nciG^3^ and the rest of each dotted 
pair is an expression called the associated INFORI-IATICR For example, 
in a meteorological data-base application, the name HONOLULU might have 
the property list 

((RAIN . 2), (HUMIDITY .40), (TEMPERATURE, 58, 96)) 

The function used in the form GET (name, key) returns the 
information which is dotted with the vaU.ue of *k^ m the property list 
of the value of "name”, returning FALSE if no such key occurred on the 
property list. 

A command of the form REMPROP (name, key) has the side effect of 
deleting from the property list of "name” the first dotted pair 
beginning with the value of "k^, if ary. REMPROP returns FALSE if no 
such indicator occurs on that property list, returning TRUE otherwise. 

A command of the form PUT (name, key, information) causes a 
conmumd of the form REMPROP (name, key) to be executed, after which the 
value of "key" dotted with the value of "information" is put on the 
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property list o£ the value o£ *naae”. PUT returns the value of 
"information”. 

Using property lists, the basic technique for accomplishing our 
various ^eratio^s on two sets of names is: 

1. For each name in one of the two sets of names, store 
TPUE under the key SEESiL 

2. For each name in the other set, check to determine 
whether or not the name has already been seen, and act 
accordingly. 


3. For each name in the first set, remove the property 
SEEN so that we won't invalidate subsequent set operations 
which utilize ai^ of the same elemat^s. 


A simpler variant of this idea is applicable to the one~argument 
functions named ISSST and MAKESEi: 

As an example, here is UNI^ defined using this technique togetlMir 
with the applicative style: % 

FUNCnCN UNION (SEIl, SEI2), 

msK ( sm ), 

UtV^HK (SEn, UNIGNAOX (SET2}} ENDFUN $ 

FUNCTION MARK (SEEL), 

WHEN EMPTY (SETl), EXIT, 

POT (FIBSr(SEIl), •SEEN, TPDE)f 
MARK (REST (SEU)) E213FUN S 
FlACnON inUCNAUX (SEr2) , 

WHEN E!<1PTY (SEEZ), SEEL EXIT, 

WHEN GET (FIRST(SETZ), ‘SEEN), UNICNAUX (REST(SEr2)) EXIT, 

ADJOIN (FIRST(SEr2), UNI0NADX(REST(SEr2))) ENDFUN $ 

FUNCTION DNMARK (SETl, ANS) , 

WHEN EMPTY (SEEL), ANS EXIT, 

REI-JPROP (FIRST(SEEL), ’SEEN), 

Uia-IARK (REST(SEEL), ANS) EM5FUN $ 

% Each time any function is invoked, the outside values of its 
parameter names, if ai^, are ”stad(ed" away to be restored later, just 
prior to return from that invocation. If a function refers to a 
variable which is not among its parameters, the: the most recent value 
of the variable on the stack is used. Thus, when UNICNAUX is invoked 
from within UNION, SEU in the definition of UNIONAUX refers to the 
argument value associated with that parameter of UNION. This treatment 
is called "dynamic binding”, and a reference such as to SEU in UNIONAUX 
is called a "fluid reference". We could have avoided this by making 
SETl be an argument and a parameter to UNIONAUX, but that would have 
made the program slightly slower and more bulky. However, fluid 
variables make programs much harder to debug and maintain, especially if 
assignments are made to them in functions other than the ones which 
establish them. Consequently, we recommend generally avoiding fluid 
variables. The (xily reasoi we used one here is to intro^ce the concept 
to issue this advice. 
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Values assigned at the top-level of muSIMP, outside all function 
definitiois, ue called 6LCBAL values. Examples are the initial values 
of muSIHF coitrol variables such as RDSr EOiO ai^ CCtCENSEr or of muMAIE 
control variables aich as FBPCH or FWPBCFD. Reference to a global value 
from within a function definition is not quite as confusing as refereice 
to a fluid value, and it is indeed cmerous to creat numerous long lists 
of p 2 irameters in order to pass such eiviroimental control values through 
a long sequence of functioi definiticxis for use deep within. 

However, here too it is at the very least considered bad 
programming s^le to unnecessarily modify aich glc^sal values from within 
a f'unction without restoring the values before exiting from the 
function. In fact it is generally bad manners for any program file to 
modify global values if the modification is merely incidental to the 
central purpose. That is why these lessons carefully save the 
prevaiiling v^ues of the control variables named ECHO and CCIOE21SE, th&i 
restore these values just prior to the end of the file. (It is truely 
annoying to have someone else's program litter your environment 
unnecessarily.) 

The property-list technique for set operations is one which we think 
is more naturally implement^ using the Von Neumann programming style. 
Try to write such a version of UNI^: % RDS: FALSE $ 

% Now, using either style, write an INTERSECTION function using the 
prc^r^list technique: % ICS: FALSE $ 

% One does not usually take the FIRST or REST of an atom 
intentionally, but th^ do in fact have well-defined valu^: ISat FIRST 
of an atom is its value, and the I&ST of an atom is its property list. 
For example: % 

WEA3HER: 'FOOL $ 

FIRST (NIATSER) ; 

POT {'WERIHER, 'TEPIPERMORE, -3) & 

POT CWEftSHER, 'WD©, '((NOPra . WEST), 30)) & 

REST ('LEATHER) 6 

% T^ is true of integer atoms too, though it is usually pointless to 
put anything on the prc^rty list of an integer, because integers are 
not stored uniquely: % 

FIRST (7),* 

REST (7); 

NINE: 9 $ 

R3T (NINE, 'TESTUG, '(1, 2, 3)) & 

GET (NINE, 'TESTING) & 

GET (9, 'TESTING) & 

% Since all nodes and atoms have a FIRST and a REST which are either 
nodes or atoms, misuse of these selectors can't accidently give access 
to the machine language, stack, print names, or anywhere else which 
could inadvertently compromise the integrity of muSIMP. Thus, 
inadvertent omission of a termination test in functions which follow 
chains of pointers is likely to be revealed by stack exhausti<mi in the 
case of recursion, and by an infinite loop in the case of iteration. 
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It is conmion practice to use EMFET to test for the end conditicai as 
a function proceeds down a list. If such a function is inadvertently 
given a iKR-Iist (i.e. a Non-FALSE atom or a structure whose final BEST 
cell points to a NckhPALSE atu»a), the function will use the FIBST cell 
of that atom (i.e. its Value cell) as an element of the list and the 
REST cell of the atom (i.e. its Property List cell) as the REST of the 
list. Generally the Property List is a well defined list so the EMPIY 
test will ultimately cause termination with no ill affects. 

We prefer to have non-list arguments give more predictable 
results confined to the argument. Thus, our internal implemaitations of 
MEMBER, REVERSE, and any Other fiinctiois ordinarily applied to lists use 
ATOM rather than EMPTY as the termination test. This is slightly faster 
too, so you may wish to generally avoid EMPTY in favor of ATOM. 
Alternatively, you can redefine EMPTY to print and return an error 
message when given a nonFM«SE atom: % 

FUNCnCN EMPTY (US) , 

WBEM Kim (US), 

WHEM EQ (US, FALSE), EXIT, 

PRINT ("*** Warning: EMPTY givei nonlist ■) EXIT 
EMDFUM $ 

EMPTY (5) $ 

% This is our first example illustrating the fact that conditional 
exits can be nested arbitrarily deep. The same is true of loops or 
blocks. This example also illustrates the fwcticn, which prints 

its one argument the same way that expressions terminated with an 
ampersand are printed. There is an analogous function named PRTMATE 
which prints its one argument same wey that e^qpressicxis terminated 
with a semicolcxi are printed. ' 

When functions are called with fewer actual arguments than the 
function has formal arguments, the ranaining formal arguments are 
assigned the value FALSE. This provides a convenient mechanism for 
automatically inserting default values for these extra argvanents. Whei 
an argument evaluates to FALSE, the function can assign the appropriate 
default value. For example, if the user omits the drive as the third 
argumait of IDS, that functi<»i uses the currently logged in drive (i.e. 
the drive indicated by the last operating system prompt given before 
altering muSIMP). 

There are instances where it is desirable to permit a function to 
have an arbitrary number of argumoits. This is accomplished by making 
the formal parameter list of a function definition be an atom or non¬ 
list rather than a list. The argum&its are passed to the function as a 
single list of argument values, from which the function can extract 
the values. For example, it is ctxxvenienb to have a function named MAX 
which returns the largest of one or more argument values. We can 
implement this as follows: % 

FQNCnCN MAX ARGUS, 
mXAUX (FIRST(ARGLIS), REST(ARGLIS)) 

EZCFUN $ 

FUNCnCN MAXAOX (BIGGEST, UNTRIED), 

WHEN a^PTY (UNTRIED), BIGGEST EXIT, 
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WHEN BIGGEST > FIESTCUNnUH)), MAXAIK (BIGGEST, REST (UNTRIED}} SIT, 
MAXMTC (FIRST(CNIRIED), REST(UinSIED)) 

ENDFON $ - 
MAX (7) ; 

JftX (3, 8, -2) ; 

% This collection of arguments into a list is called NOSPREAD, to 
distinguish fri^ the SPREAD brand of peanut butter. 

More generally, nuMATB permits a coidsination of the 2 techniques: If 
a parameter-list is a dotted-pair of two names or a list whose last 
element is a dotted pair of two names, then the last parameter name 
accumulates a list of any excess arguments beyond those spread to the 
other parameter names. Thus, we can simplify our ^inition of MAX to: 

% 

FUNCnCN MAX (FRST . GIBERS), 

MAXAUX (FRST, OIHEBS) 

ENDFUN $ 

% Would you like to try this technique? Appropriate candidates 
include MIN, UNICN, and INTER^XmCN. % RDS: FALSE $ 

% Now, suppose that for some reasai we alreac^ have a list of integers 
such as % 

NUMBLIS: '(18, 3, 7, 91, 12, 2) $ 

% and we want to find their maximum. The expression MAX (NUMBLIS) 
will not work, because MAX is designed for numeric azqm&ntSf not for a 
list of numbers. We could of course extract the elanents and feed than 
individually to MAX, but this is awkward, especially if we are rt^erring 
to MAX inside a function and we* do not know ahead of time how many 
integers are in NUMBLIS. Fortunately there is a convenient function 
named APPL!?, which aq^lies the function whose name is the value of its 
first argument to the argument list which is the value of its second 
argument. Consequently, we need merely write % 

APPL^ ('MAX, NUMBLIS) & 

% APFLX works on either SPREAD or tICSPREAD functions. Why don't you 
try out a few examples: . % RDS: FALSE $ 

% A function written in muSIMP-79 is stored internally as a nested 
list, and the function named GETD returns a pointer to this list. 
Consequently, to see what the internal representation of UNION looKs 
like: % 

GEXD (UNICM) & 

% GEZD returns TRUE if the def initicx: is in machine language, and GEXD 
returns FALSE if there is no function definition for its argument. 
Those who are curious may wish to use this functio: to experimentally 
determine the correspixidaice between the external and internal forms of 
a function definition. This can be useful for revealing bugs arising 
from misconceptions about how the parser regards certain constructs. 
All we want to point out here is that since function definitions are 
represented as lists, muSIMP functions can easily operate upon other 
muSIHP functions. This makes it easy to write muBIMP programs which 
service other muSIMP programs. Examples include muSIMP-oriented 
editors, cross-reference programs, debuggers, verifiers, statistics- 
gatherers, pretty-printers, file comparators, emd compilers. The 
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internal r^resentatioi also makes it possible for functicns to modify 
each other dynamically, as they execute. Ihe implications for artifical 
intelligence are intriguing to contemplate. 

SEMD is a related command which clears any function definition 
existing under the name which is the value of its argument. For 
example, % 

REMD (UNXGM) & 

GEm (CNION) & 

% One good use of REMD is to free space occupied by functions which 
are no longer needed immediately, in order to provide enough space for a 
more urgent need. For example, suppose that in muMATE a problem 
requires the SdiVE package followed by the MAIEIX package, but there is 
not room enough for both packages to coexist in the amount of memory 
presait on the machine. Then, after usii^ the SOLVE package but before 
reading in the MATRIX package we could remove function def initials for 
SOLVE by commands such as 

REZ^ (SOLVE) $ REMD (SOLSCP} $ ... 

Less typing would be involved if we defined a command named 
MULUREMD, which for an argument which is a list of names, »iccessively 
allies REMD to each name. In this example it is the side efects rather 
than the returned value which is of interest, so MULTIREMD can return 
whatever is the least trouble. MULTIREMD is trivial to write, using 
either recursion or iteration, because the same "program scheoa" occurs 
so often: Walk down a list, successively applying a function of one 
argument to each element of the list, then return anything. This 
observation leads to the following idea: Let's write a functioi which, 
given the name of any function of 1 argument, together with a list, 
successively applies the function to the elements, thcsi returns anything 
conveiient: % 

FUNCnCN (FUNNMffi, LIS), 

LOOP 

WHEN e®Ty (LIS), EXIT, 

APPLY (FUNMAME, FIRSr(LIS)), 

LIS: REST (LIS) 

QOLOOP 
ENDFUN $ 

% Then, for example, we could write 

HAP ('REND, '(SOLVE, SCLBCP, ...)}$ 

What we have done is to separate the general-purpose control- 
sequence from the specific tasks which can use it. This division of 
labor accomplishes two useful things: 

1. Program space savings can accrue for each use of MAP with a 
different function, beyond the first, because essentially 
diiplicate coitrol sequoices are avoided. 

2. Once the meaning of MAP becomes familiar, the program is 
more readable, because MAP ('REMD, '(SOLVE, SOLEXP, ... )) is 
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then instantly understood to mean R£MD all of SOLVE, SOLEXP, 
etc.. In contrast, the altennative form MULTIHEMD ('(SOLVE, 
SOLEXP, ... I) requires the user to check the definition of 
HULTIBEMD to be sure the purpose is correctly understood. 

Another frequent need is to walk down a list, applying a function of 
me argument to each element, but return the list of results. Write a 
mapping function of this kind, called MAPLIST since it returns a list. 
Then, try 

MAPLIST (•-, '(3, 8, 14)), and MAPLIST ('NOT, '(TRUE, FALSE, MAYBE)) 
%RDS: FALSE S 

% MAP and MAPLIST are the most widely applicable mapping functions, 
but if you grow to like mapping functions you may develop a large suite 
of them. For example; 

1. For functions of two arguments you could have a map 
function of the form MAP2 (function name, listl, list!) or 
MAP2 (function name, list of pairs) . Since much of muMATE is 
stored on property lists, this could be used to apply REMPROP 
appropriately to help delete high-level muMATH packages in 
order to make ^ce. (Here is an idea: for each muMAIH file, 
write a corresponding file of type DEL, which has an 
appropriate command of the form MAP CREMD, ... ) , together 
with one of the form MAP2 ('REMPROP, ... ). Then, to delete 
the SOLVE package froa memory, one merely issues the command 
RDS (SOLVE, DEL, drive).) 

2. For functions of two arguments you could have a mapping 
function used in the form ^^2LIST (function name, listl, 
list2) orMAP2LIST (function n^uae, list of pairs), which is 
like MAP2 but returns a list of results. 

3. E^r general trees you could have a mapping functioi called 
TREEMAP which applies a function to the atoms in a tree, and 
there could be a similar one called TREEMAPTREE which is 
similar but returns a tree. % 

ECHO; #ECHO $ CCtDENSE: #Ca©ENSE $ 

RDS 0 $ 
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wuSIMP-79 Primitive Data Structures 
fhe Soft Warehouse 11/26/79 


I. Dm SZPOCIUBES. 

muSIHP'data is coisprised of names, numbers, and nodes. Each type 
is recognizable and consists of a fixed number of "pointer* cells 
containing memory addresses. The cells can either ^int to other 
objects or to special-purpose entities outside the pointer space of 
objects. However, all three types have a PIPJT cell and a REST cell- 
Moreover, these PIRST/REST cell pairs can only point to other objects 
within the pointer space. This eliminates the need for time-consuming 
run-time type-checks in the crucial selector functions which fetch these 
pointers. 


i ... - I . . I. I II I ■ I 

A. Names | value | Property | Punction | Fnames i 
I.. . . "■ I . .. - 1 - 1 . ' . .— I' 

A name is a recognizable, structured object consisting of four 
pointer cells. Names are uniquely stored so that di^licate names camot 
coexist in storage. Here are the uses of the four cells: 

1. The FIRST or value cell contains a pointer to the 
name's current value which is used by the evaluation 
functions. The value of a name is initialized to a self¬ 
reference of the name; however, it is modified by the 
assignment functions and when the name is used as a formal 
parameter in a function definition. 

2. ihe REST or pr^r^ list cell contains a pointer to 

the name's property list which is used by the property 
functions. El«Dent8 of this list are indicators dotted with 
the corresponding values. Property lists are initially set to 
the list. 

3. The Punction cell contains a pointer to the name's 
function definition if any. The contents of this cell can't 
be accessed except as functi^ applications, and the contents 
can't be modified except by means of the functicxi definition 
primitives. When a new name is first created, its function 
cell is initialized to the undefined-function trap routine. 

4. The Pname cell ccntains a pointer to the name's ASCII 
print-name string, which can be of arbitrary length. Acc e ss 
to this cell is restricted to the I/O and sub-atomic 
primitives. Print names are defined when a name is first 
used, and they cannot be modified or expunged. 
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B. NUnbers 


—- 1 —^—-- 1 --- 1 - 

I Self i FALSE | Vector | 

+----1-H- . - . I 

A number is a recognizable, structured object oxisisting of three 
pointer cells. Numbers are not uniquely stored, so duplicate numbers 
might coexist in storage. The cells are used as follows: 

1. The FIRST cell contains a pointer to itself. 

2. The REST cell is initialized to FALSE. 

3. The Number Vector cell contains a pointer to the 
actual number, which oxisists of a signed vector of up to 254 
bytes. Thus, the magnitude of numbers is limited to 256^254, 
which is approximately 10*611. 


. 1 I 

C. Nodes I FIRST I REST | 

---^ 

Bina]^ trees are the primary data structure in muSIMF. Internally 
they are implemented as a network of cell pairs called nodes. Each node 
consist of a FIRST cell and a REST cell. As mentioned earlier, the 
node's cells can only point to other bonified muSIHP data objects; 
either a name, a number, or a node. Nodes are often called "dotted* 
pairs”, because of their linearized extemed. notation produced by PRINT 
or accepted by READLIST: notation 

(X . Y) 

represents a node whose FIRST cell points to the object X, and whose 
REST cell points to the object Y. Although the dot notation is more 
general, it is often more convenient to think of data as a linear list 
than as a deeply nested binary tree. For this purpose, lists are 
recursively defir^ as follows: 

1. The emp ty list is demoted by the name FALSE. 

2. If Y is a list and X an object, then (X . Y) is a list. 

A list of objects is printed by the function PRINT as a sequence of 
its elements separated by commas and delimited by parenthesis. The 
functio: REALIST recognizes this notation for ii^t. For example, if Y 
is the list (Yl, Y2, ..., Yn) then the dotted pair (X . Y) is printed as 

(X, Yl, Y2, ..., Yn) 

Conversely, the input of the form (X, Yl, Y2, ..., Yn) is recognized as 
(X . Y) by the READLIST function. 
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II. mm MANAGaen. 


P^namicr transparent aenory management gives muSIMP much of its 
inherent power. Ideally, at any given time during the execution of a 
program, all of tht memocy not actually required to decribe the state of 
the machine should be available for any subsequent program use. This is 
approximated in muSIKF-79 by first partitioning the available resources 
into ti» ^^iojs data-spaces and then recycling storage within each of 
these ^aces as required. Normal stack operations continuously reclaim 
the stack space; whereas, an automatically invoked garbage collector 
reclaims the remaining ^s^es. 


A. Initial Data-space Partition 

During the initialization ^lase of muSIMP, the amount of reai^write 
meaory available to the int^reter is first computed. Hemocy is then 
partitioned into four distinct data^spaces using the following 
proportiOTS: 

4/32 Atom j^ce Name and number pointer cells. 

3/32 Vector ^ce Print-name strings and number vectors. 

23/32 Node ^ce Node cell pairs. 

2/32 Stack Sp&c^ Control/value stack. 

Based on our experience, these proportions provide a reasonable 
balaix:e between the spaces for most applicati<m^. 


B. Garbage CoUectioi 

New data structures are generally constructed dur^ the execution 
of a muSIHF program, while others are implicitly discarded as they 
become un-referenced. When the construction process uses up all 
available resources, a garbage collector routine is called to r^aim 
the storage space vacated by discarded data structures, so that the user 
prograun can continue, Jn muSIMP-79 the exhaustion of resources in 
either the atom, vector, or node spaces will cause collectiai to o^r. 
Those data structures accessible by means of chaining through pointer 
cells beginning either from a name cell or from a value stack entry are 
marked. Th^ durii^ the secotnd pass all the unmairked r^es and numbers 
are collected for re-use, while simultaneously removing the mark on the 
j^c^ible nodes. 

Although garbage collection is automatic, it is not entirely 
invisible to the user since it periodically causes a pause in the 
execution of a program. About 1.5 seconds is required for the 
collection process in a 48K byte muSIMP-79 system using a 2MHz CPU 
clock. Normally this is of no concern to the programmer; however, it 
shoild be considered in the design of real time systems. A phmomenon 
know as thrashing occurs when the system is forced to spend an 
inordinate amount of time garbage collecting for a very small amount of 
nodes. This can be resolved by increasing the computer's menory size or 
decreasing the amount of program and data storage requirements. 
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m. ERFGR And HTTERRUPT TRAPS 


I£ there is a reasonable interpretation for a construct, auSIM? 
generally uses it. Consequently, error traps are induced only by 
situations for which there is no satisfactory recovery. Exaa^les are 
the exhaustion of available data ^ace or disk I/O errors. Cn the other 
hand, a software internjqpt is caused by ah internet character (i.e. an 
ESC, ALT, or Ctrl'Z) receivfd from the terminal, which can be sent at 
ai^ time, vihesi a particular trap occurs, the ^^rqpriate diagnostic and 
the following "cations” message ue sent to the terminal: 

EXEOJnVEi ESC, ALT, Ctrl-Z; RESTART: RIB, DEL; STSIEl: CtrlHl? 

The user may then type one of the appropriate alternative option 
characters. The "ElCECtJTIVE" ^icn is the least drastic siixse it merely 
causes control to retum to the muSIMP executive driver lo^, withexjit 
changing function definitions, property values, or name values, from 
what they were just prior to the interrupt. The second option destroys 
all non>primitive muSIHP functions, pre^rty values, and name values, 
then restarts muSIMP afresh. Finally, the "SYSTEM" option terminates 
muSIMP, and returns coitrol to the operating sytem. 


A. Data ^ce Overflow 

As discussed in Section II, there are four distinct data ^aces in 
muSIMP to accommodate the various data types. Normally, aut^oatically 
invoked garbage collecticxis will provide sufficient q^e in each area 
to contimously satis^ the demaix3s of user programs. However, in the 
event all of the available resources in an area become exhausted, an 
error trap will occur and one of the following diagnostic messages will 
be displayed on the terminal: 

NCDE Space Exhausted 
ATCM ^ace Exhausted 
VSCIQR ^ace Exhausted 
STACK Overflow 


B. Disk Pile I/O Errors 

Disk errors may be caused by inaifficieit disk space, attempts to 
read past the «xi-of-file, or hardware malfuncti<xis. The read and write 
disk error diagnostics are respectively: 

E:^ of Pile or READ Error 
No Disk ^pace or tfRTIB Error 
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C. undefined Nymerical Operations 

If the second argument to any of the functions QUOTIEKTr HOOr or 
DIVIDE is Cn a zero-divide trap occurs with the following diagnostic: 

ZEEC Divide Error 

D. Input §^tajt Error 

The only syntax error trap caused by the function READ is when a 
closing right parenthesis is not found when using the 'dot' notation. 
The diagnostic is: 

Input Syntax Error 

Fimction PARSE can produce syntax error traps together with 
diagnostics of the following forms: 

*** snmx ERRCE: expression USED AS NAME/ 

*** SXNIAX ERRCR: expression USED AS PREFIX OPERATOR, 

*** SYNTAX ERROR: expression USED AS INFIX OPERATOR, 


*** SYNTAX ERROR: delimiter NOT FOUND, 

where "e^^jressiai'' is the apparent offending portiai of the input, and 
where "delimiter" is an apparently missing right delimiter such as a 
right parenthesis, ENDFUN, ENDSUB, EXIT, ENDLOOP, or ENDBLOCK. In any 
event, the remainder of the input from the point of confusion through 
the next terminator, such as or is output to the terminal 

to help itxiicate the probable neigi:^rhood of the cause. Examples which 
provoke the above four types are respectively: 

5 (X); {perha^ 5*(X) was intended?} 

X Y; {perhaps x*Y was intaxaed?} 

X*A; {perhaps X/Y was intended?} 


WHEN ATOM(X, EXIT {perhaps WHEN ATOM(X) , EXIT was 

intended}. 
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17. PRUCCnVEIZ DEFINED FONCTICNS. 


'Ike muSIMP C&tructured IMP lementation^ Language is a high level 
computer language ideally suited for symbolic and semi-numerical 
processing. Currently, it is impl^ented by means of a bootstrap file 
MUSMOBEJiUS which is automatically loas^ prior to using the language. 
For an interactive introducticn to the features available in muSIMP, the 
tutorial lesson files, beginning with FLES1.IRA, nu^ be executed. See 
the description on how to take the programming lessons in LESSOiS.'DCI. 

Every language must be described in terms of some language, which 
must be described in terms of some language, etc. Thus it is clear that 
at some point we must s^peal to etssumed inborn or culturally acquired 
understanding. This unnecessary sequence of "buck passing" can be 
avoided by using a somewhat circular Ascription of muSXNP. In other 
words muSIMP-79 can be described in terms of muSIMP supplemented with 
English where necessary. Use of such a description requires some 
prerequisite knowledge of imiSINP gained by other means, just as use of 
an English dictionary requires some prerequisite knowleA® English. 

After one has initially learned the basics of muSIMP from the 
lessons, this type of reference manual has the advantage of being 
compact while requiring mastery of no auxiliary notations. In addition 
it provides excellent, nontrivial examples of structured programs 
written in muSIMP. 

The following is a description of all of the primitively defined 
user-level functions, operators, control constructs, and control 
variables in muSIMP-79. For descriptive purposes only, we introduce 
some fictitious functions which are unavailable to the user. Th^ are 
indicated by being unnumbered and are also unindexed. Lower-case type 
is enployed where English is used rather than legitimate muSIMP program 
constructs. 
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A. Selector Functic»is 


1. FUNCnCN FIRST (X), 

tae contents of the FIRST cell of X, 
BBFUN; 

Interpretatlcns: 

9, ^e first item of a list Xr 

b. The left element of a dotted-^air X, 

c. The value of an atom x. 


2. FONCTICN REST (X), 

the extents of the REST cell of X, 
ESCFON; 

Interpretations: 

a. The tail of a list X, 

b. Tte right element of a dotted-pair Xf 

c. The pr^rty list of an aton X. 


3. FUNCnOi SBCCND (X), 
FIRST (REST (X)), 


4. FUNCTION RREST (X), . 
REST (REST (X}}, 
QDFUN; 


5. FUNCTION THU© (X), 

FIRST (REST (REST (X))), 
E1©FUN; 


6. FUNCTICN RRREST (X), 

REST (REST (REST (X)}}, 
ENDFUN; 



B. Constructor Functions 


1. FUNCTION ADJOIN (X, Y), 

a new cell-pair wnose FIBST cell is X and whose REST 
cell is Y, 

ENDFUN; 

Interpretations: 

a. A list whose first elonent is X and whose tail 
is Y, 

b. A dotted-pair whose left element is X ai^ whose 
right elessent is Y. 


2. SUBRCOTINE LIST (Xl, X2, Xn), 

NHEN n « 0, FALSE EXIT, 

ADJOIN (EVAL (XI), LIST (X2, X3, ..., Xn)), 
ENDSUB; 

Interpretatiai: The list (XI, X2, ..., Xn). 


3. FUNCTION REVERSE (X, Y), 

HHEN ATOM (X), Y EXIT, 

REVERSE (REST (X), ADJOIN (FIRST (X), Y)), 
ENDFUN; 

Interpretation: The reverse of the list X. If a 
second argun^t Y is givei, the reversed list is 
appended to the beginning of the object Y. 


4. FONCnCN OBLIST (), 

a list of the current built-in and user-introduced 
names, 

ENDFUN; 

Interpretation: The object (name) list. 
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C. Modifier Functiois 


1. PUNCTICN REPLACEF (X, Y) , 

PISSr cell of X; Y, 

X, 

ENDFUN; 

Interpretations: 

a. Replace the first elesi^t of a list X by Y/ 

b. Replace the left element of a dotted-pair X by Y/ 

c. R^lace the value of an atom X Y. 


2. FUNCnON REPIACER (X, Y), 

REST cell of X: Y, 

X, 

E2DFIJN; 

Interpretations: 

a. Replace the tail of a list X by Y, 

b. Replace the right element of a dotted-pair X 
by Yf 

c. Replace the property list of an atom X by Y. 


3. FUNCnCN CCNCfilEN (X, Y), 

WHEN ATOM (X), Y EXIT, 

WHEU POm (REST (X)f, REPLACER (X, Y) EXIT, 

CCtJCAIEN (REST (X), Y), 

X, 

ESCFUN; 

Interpretatim: Concatenate, without adjoining, the list 
Y onto the right end of the list X. 
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D. Recognizer Functions 


1. FUNCnON NAME (X), 

X is a name, EXIT, 

BDFON; 

Interpretation: Recognize objects which are names. 

2. FDNCnON INTEGER (X), 

ViQQi X is an integer, EXIT, 

SESFUN; 

Interpretation: Recognize objects which are integers. 

3. FUNCnCN ATOM (X), 

(X) OR INTEGER (X), 

ENDFUN; 

Interpretation: Recognize objects which are atoms. 

4. FmcnCN E^PIY (X), 

X * FALSE, 

2NDFUN; 

# ' 

Interpretation: Recognize the empty list. 

5. FUNCnCN POSITIVE (X), 

X > 0, 

£20FUN; 

Interpretation; Recognize positive integers. 

6. FUNCnCN NEGATIVE (X), 

X < 0, 

ENDFUN; 

Interpretation: Recognize negative integers. 

7. FONCnCN ZERO (X), 

X « 0, 

ENDFUN; 

Interpretatioi: Recognize zero. 

Note: All recognizers return TI^ or FALSE. 
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E. Comparator Functions and Operators 


1. PONCnCN BQ (X, Y), 

WEBi INTEGER (X) AID IKESGER (Y), X - Y EOT, 
VSEN X and Y point to the same abject, EXIT, 
EIDFUN; 

Interpretation: The identity comparison of X and Y. 


2. PBDPERry RBP, », 80; 

PROPERTY I£P, *, 80; 

FUNCnCN * (X, Y), 

WHQJ ATOM (X), BQ (X, Y) EXIT, 

HBEH ATOM (Y), FALSE EXIT, 

WHEN FIRST (X) ■ FIRST(Y), REST(X) * REST(Y) EXIT, 
ENDFDN; 

Interpretation: The infix equality cperatx^r, «, treats 
X and Y as being equal if and only if th^ have 
isomorphic tree structures with identical atomic 
terminid nodes. 


3. FUNCTION CRDERP (X, Y), 

WHEN the address of the object X is less than the 
address of the'object Y, EXIT, 

ENDFUN; 

Interpretation: A generic ordering function for system 
- names based on their order of introduction. 


4. PEQPERTY RBP, >, 80; 

PROPERTY LBP, >, 80; 

FUNCnCN > (X, Y), 

WHEN INTEGER (X) AID INTEGER (Y), X > Y EXIT, 

msem; 


5. KDPEKPY ms, <, 80; 

PROPERTY LBP, <, 80; 

FUNCnCN < (X, Y), 

WHEN INTEGER (X) AND INTEGER (Y), X < Y EXIT, 
ENDFUN; 

Note: All cooparators return TRUE or FAL^. 
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F. Logical Operators 


1. PROPERTY RBP, NOT, 70; 

FONCnCN NOT (X), 

X * FALSEf 

msFuti! 

Interpretation: NOT is a prefix operator with right 
binding power 70. 


2. PROPERTY RBP, AND, 60; 

PROPERTY LBP, AND, 60; 

SBROOnNE AND (XI, X2, Xn), 

MBEN n « 0, EXIT, 

WHEN NOT EVAL (XI), FALSE EXIT, 

Ate (X2, X3, ..., Xn), 
qOSUB; 

Interpretation: AND is a logical infix operator with a 
left and right binding power of 60. 


3. PROPERTY RBP, CR, 50; 

PROPERTY IfiP, OR, 50; 

SUBROUTINE CR (XI, X2, ..., Xn), 

WHEN n - 0, FALSE EXIT, 

WHIN EVAL (XI), EXIT, 

OR (X2, X3, ..., xn), 

eesuB; 

Interpretation: CR is a logical infix operator with a 
left and right binding power of 50. 


Note: All logical operators return TRUE or FALSE, and ai^ 
nonFALSE logical operand has the same effect as TRUE. 
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G. Assignment Functioxs 


1. FUNCTION ASSIGN {X, Y) , 

FIRST cell of X: Y, 

Yr 

BDFUN; 

Interpretation: Set the veOue of the atom X to Yf 
and return Y. 


2. PROPERIY RBP, 20; 

PROPERTY IBP, ISO; 

SUBRCXJTINE : (X, Y), 

ASSIGN (Xr EVAL (Y)), 

QCSUB; 

Interpretatioi: Set the value of 'X to Y, and return Y. 
*:'* is an infix <^rator with left binding power 
180 and right binding power 20. Evaluation of this 
form returns the v 2 due of the expressioir after 
achieving the side effect of assigning the value 
of the expression to the nesne. 

Note: Assignmeits to non-names are allowable, having an 
effect similar to REFLACEF/ but returning a different 
pointer. 


13 



H. Proper^' Functions 


1. FUNCnCN ATSOC (X, Y) , 

WHEN ATCM (Y), Y EXIT, 

WHEN ATa-l (FIRST (Y)) , ATSOC (X, REST (Y)) EXIT, 
WHEN BQ (FIRST (FIRST (Y)) , X), FIRST (Y) EXIT, 
ATSOC (X, REST (Y)), 

E2DFUN; 



Interpretation: The first non-atomic object on the 
"asSCCiation" list Y whose i^Toroic FIRST cell is X. 


2. FUNCnC3N GET (X, Y), 

X: ATSOC (Y, REST (X)), 

WHEN psm (X), FALSE EXIT, 

REST (X), 

EtOFQN; 

InterpretaticKi: The pr(^rty value associated with the 
indicator Y on the proper^ list of X. 


3. FUNCnCN PUT (X, Y, 2), 

WHEN EMPTY (GET (X, Y)), 

REPLACER (X, ADJOIN (ADJOIN (Y, 2), REST (X))), 
2 EXIT, 

REPLACER (ATSOC (Y, REST (X)), 2}, 

Z, .. ' 

ENI^UN; 

Interprstatical: Place <xi the property list of the &tm 
X under the indicator Y the property value 2, 
destroying aiy previous value laider the same 
indicator. 



4. FUNCTION REI4PRDP (X, Y), 

WHEN Aral (REST (X)), REST (X) EXIT, 

WHEN EQ (FIRST (SECCND (X)), Y), 

Y: REST (SECCND (X}), 

REPLACER (X, RREST (X)), 

Y EXIT, 

RS'lPROP (REST (X), Y), 

ENDFUN; 

Interpretation: Remove from the proper^ list of X 
the property value associated with the indicator Y. 
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. PBOPBFTSf namer atour value 

read two names X and Y, then parse an ca^ressicn Z. 
then without evaluating any of thra, place Z on the 
property list of X, under key Y. 
th«i return Z; 

PROPERTY is a data-base construct which returns a list of name 
and atom after accomplishing the side effect of storing the 
value on the property list of the first operand, under the key 
which is the second operand. Any previous value on that 
property list under the same key is deleted, with a 
correspcr^ij^ warning message. The three (grands of PROPERTY 
are automatically quoted, so that, for example, they can be 
unquoted aerators. 
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1. Defiiiitiai PuxK±ions 
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1. FUNCnCN GEID (X), 

WHEN ICT (NAME (X)), FALSE BCIT, 

NHEN X is not a defined function, FALSE EXIT, 

WHQ{ tile functicn cell of X points to a machine 
language function, EXIT, 
the object pointed to by the functicn cell of X, 
QOFUN; 

Interpretation: The definiticxi of the function named X. 


2. nacnoN puiD (x, y), 

WEEN NOT (NAME (X)), FALSE EXIT, 
functicxi cell of X: Y, 

V, 

ENDFUN; 

Interpretation: Place a pointer to the definition Y in 
the function cell of X. 


3. FUNCTION MDVD (X, Y), 

WHEN NOT (NAME (X)) OR NOT (NAME (Y)) , FALSE EXIT, 
function cell of Y: functicai cell of X, 

GE2D (Y), 

ENDFUN; 


Interpretatioi: Copy the ctefinitial of the functicn 
named X to Y. 


4. FO^XCTION REMD (X), 

WHEN NOT (NAI^ (X)), FALSE EXIT, 
functioi cell of X: undefined, 

GEID (old definition of X), 

ENDFUN; 

Interpretation: Remove the functioi definition froa X. 


5, PBOPERry PREFIX, FUNCTION, 

parse a function definitioi, then use ECTD to put it 
in the function cell of the function name, 
then return the function name. 

Interpretation: FUNCTICN is ti« l e a d i n g keyword of a 
control construct which has the gei^ral form: 
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FUNCTION name parameters, 
taskl, 
taslc2f 

taskn 

QDFUN 

The name can be omitted when there is no need to refer to it, 
such as when a nonrecursive function is stored on a property 
list for use by APPLY, "parameters" can be an arbitrary 
symbolic expression. V?hen "parameters" is any name, except 
"FALSE", the function ma’- ; absequently called with an 
arbitrary number of argi." .d passed to the function as a 

list assigned to the argu.; ..... If "parameters” is not a name, 
the first argument in a call to the function is assigned to the 
FIPST of "parameters" and the REST of the arguments is assigned 
to the REST of the "parameters" in an identical manner. 

Task evaluation within the function is performed successively 
until either the end of the tasks is reached or a non-FALSE 
predicate is evaluated. In the latter case evaluation proceeds 
as before except down the predicates task list. In either case, 
the value of the function is the vcdue of the last task 
evaluated. 
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J. Sub-atonic Functions 


1. FUtOICN CCMPPESS (X), 

WHEN A2CM (X), "" EXIT, 

WHEN NAME (FIRST (X)), 

coicatenate the print name of FIRST (X) onto the 
beginning of COMPRESS (REST (X)) th&i return the 
corresponding muSS'lP name, 

COMPRESS (REST (X)), 

ENDFHN; 

Interpretatioi: The atoc whose print name is the 
padced version of the names in the list X. 


2. FtJNCnON EXPLCDE (X), 

WHEN NAME (X), 

a list of n a m es whose print names correspond 
to the chairacters in the print name of X, EXIT, 

ENDFON; 

Inte^retation: A list of the characters, in order, 
in the print name of X. 


3. FUNCTION LEM31H (X), 

WHEN NAt4£ (X), 

WHEN EMPTY (X), 6 EXIT, 

the number of characters in the print name of X EXIT, 
WHEN INTEGER (X), 

the number of bytes in the vector of X EXIT, 

1 + LENGTH (REST (X)), 

ENDFUN; 

Interpretations: 

a. The nuE^r of characters in the name X, 

b. The number of bytes in the number X, 

c. The number of top-level itons in the list X. 
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K. Numerical Functiois 


1. FUNCnCN M32IJS (X), 
WB£2i mrESm (X), 
-X EOT, 

QDFUN; 


2. FDNCnCN ELDS (X, Y), 

WEE7 ir’iSGER (X) AND INXESER (Y), 
X + Y EOT, 

EEFDN; 


3. PUlCnCN DIFFERENCE (X, Y), 

WHIN INTEGER (X) AND INTEGER (Y), 
X - Y EXIT, 

EOFUN; 


4. FUNCnCN TIMES (X, Y), 

WHIN INTEGER (X) AND INTEGER (Y), 
X * Y ECIT, 

ENDFON; 


5. FUNCnCN QOCTIENT (X, Y), 

WHEN INTICER (X) AND INTEGER (Y), 

WHEN Y ■ 0, z6ro-divide error-tra^ ECTT, 
WHEN PCSrnVE (Y), floor (X/Y) BCIT, 
ceiling (X/Y) EXIT, 

ENDFDN; 


Note: file integer quotient which is ^isistent 

with MCD being a periodic nonnegative remainder. 


6. PUNCnCW MCD (X, Y), 

X - (Y * QDCTIENr{X,Y)), 
ENIXLN; 


7. FUNCHCN DIVIDE (X, Y), 

WHEN INTEGER (X) AND INTEGER (Y), 

ADJOIN {QOOnENT (X, Y), MCD (X, Y)) EXIT, 

ESOFQN; 

Note: All muSIMP-*79 numerical functions return FALSE if either 
of their eurguments is non-numeric. 


19 


8. HMCnON + (X, Y), 

ELDS (Xr Y), 

ENDFDN; 

EROPEKTY +, PREFIX, PARSE (SCAN, 130) ; 



PROPERTY +, RBP, 100; 
PROPERTif +, LBP, 100; 


9. FUNCnCN - (X, Y), 

WHEN Q®TY (Y), MINDS (EXl) BCIT, 

DIFFERENCE (X, Y), 

ENDFDN; 

PR OPE RT Y PREFIX, -, LIST PARSE (SCAN, 130)); 

PROPERTy RBP, -,100; 

PROPERTY LBP, -, 100; 


10. FUNCTION * (X, Y), 

TIMES (X, Y), , 

ENDFDN; 

PROPERTY RBP, *, 120; 

PROPERTY LBP, *, 120; 

U. FDNCnON / (X, Y), • ^ 

QDCTIENr (X, Y), 

ENDFDN; 

PROPERTY RBP, /, 120; 

PROPERTY IBP, /, 120; 
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L. Iteader Functions 


1. ETOCriCN REWXHAR (), 

read cne character from the current input file and 
return the corresponding nuSBl? atom. Integer atoms 
aze returned only if the character is a decimal 
digit less than the curr^t base. 

ENKDN; 


2. FUNCTION SCAN (), 

read ana atom from the current iiiput file and 
return the corresponding muSIMP atom. Atoms are 
delimited either separator or break characters, 
however, the latter also are returned as atoms 
themselves. 

£ND; 

Separator characters: space, carriage return, line feed, 
and t 2 do (control-I). 

Break characters: ! $&'()* + ,- ./ 

@ \ I }' 


3. FUNCnCN BEAD (), 

read from the current input file oie caii>lete 
expression writtei in dotted-^>air and/or list 
notation, then return the corresponding generated 
dsject. (Atoms are delimited by either separator 
or break characters, but the latter also are returned 
as atoms themselves). 

EJDFUN; 


Separator characters: space, conma, carriage return, 
line feed, and tab (control*!). 

Break characters: . ) ( 

Note: Extra right parentheses aixi dots are ignored. 


4. FUNCnCN PARSE (EXl, RBP, EX2), 

from the current iapat file, read the complete 
eapressicm including the already*read token QQ, where 
RBP is the right bindii^ power of the aerator to the 
left of £X1, if any, thmi return the resulting 
unevaluated object. (DC2 is a local variable w'nich 
accumulates the parsed representation.) 

QDFUN; 
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Note: When two cperatocs compete for an operand between 
them, the operator with higher binding power towards 
the operand acquires the operand. In the case of a 
tie, the operator on the left acqpjires the ^erand. 

. PCNCnCN SXN13« X, 

print " *** SYNTRX ERRCR: ", thei print each elemeit 
in the list of arguments, X, timi read up through the 
ne^ terminator, echoing the input beginning on a new 
line, tbm return FALSE, 


6. FUNCTXCN HA3!C2 (DSLXH), 

% Uses ti» fluid variable SCAN set by SCAN () % 

WHEN SCAN > DELIM, SCAN (}, FALSE EXITr 
WEEN SCAN » coRina, SCAN (}, MAICH (DEUH) BUT, 
WHEN DELIMITER (), SSmX. (DELIM, "NOT FOa©") EXIT, 
ADJOIN (PARSE(SCAN,0), MAOXS(DELIM)), 

ENDFUN; 


7. DELIMITER: '(EXIT, ENDFUN, ENDLOOP, ENDBLCCK, 
right parenthesis, cooma) & 

HacnCN DELIMITER (), 

1ERMIN2m3R () OR MB4BER (SCAN, DELIMITER), 
E2BFUN; 


8. FUNCnCN TERtTOOOR (), 

SCAN - OR SCAN - •$ OR SCAN* 
ENDFUN; 


10. ECHO: FALSE;. 

FUNCnON ECHO (), 

NOT BDS OR ECHO, 

ENDFUN; 

Interpretation: ECHO () is a function which returns TRUE 
if input is being echoed to the terminal. 


11. RDS: FALSE; % Device ReaD Select % 

FUNCnCN PDS (X, y, Z), 

WEEN EMPTY (X), RDS: FALSE EXIT, 

WHEN NAME (X) AND NAtC (Y), 

WHEN EMPTY (Z), 

if there exists a file named X.Y on the 
currently logged disk drive, then ap^ 
that file and RDS: X, EXIT, 

WHEN NAt^E (Z), 
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i£ there exists a file named X.Y on drive Z, 
then open that file and SDS: X, EXIT EXIT, 

QDFUN; 


Notes; 

1. Normally control of the current input file is> dene through 
the use of the function EDS as described above. However, after a file 
has been opened and made current, control can be returned to the 
console without closing the input file, simply setting the value of 
EDS to FALSE. A subsequent non-FALSE assignment to RDS will thei return 
control to the point in the opened disk file at which reading was 
suspended. 


2. If the console is the current input file and all the 
characters have been read from the current line, the operating system's 
line-edit routine is called for further input, laius the system's normal 
line-edit features will be used until a carriage return is typed, at 
which point muSII4P will regain controL 

3. If a disk file is the current input file and the BCF (end- 
of-file) character is read, an error message is sent to the console, 
the console is made the current input file, amd an error-options trap 
occurs. 

4. If a disk file is being read and the value of the name 
ECHO is non-FALSE, the characters being read are also echoed to the 
current output file. 

5. Comments in an input source file must be delimited by 
matching percent signs. The text of the comment will then be ignored by 
the functions SCAN and READ, except to possiblly echo the comment as 
described in note 4 above. 

6. Special characters such as the comment, separator, and 
break characters can be. read in as names or parts of names by means of 
quoted strings. Such strings are delimited by double quote marks. The 
double quote can be included within the string by using two adjacent 
double ^otes for each desired internal double quote. 

7. As an added programming convenience the muSII'lP name SCAN 
is alws^ set to the most recently read atom. 

8. Lower-case letters are legitimate and distinct from their 
upper-case counterparts. The only exertion to this is file names and 
types giv^ to the functions RDS and WRS. They eze always converted to 
upper case in order to eliminate ca' 4 flicts with the operating system's 
file naming convention. 
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M. Printer Functions 


1. FUNCnCN PRINT (X), 
whet; name (X), 

output the print name of X to the current 
output file, EXIT, 

. WHEN INTEGER (X) , 

output to the current ou^t file the digits of X 
stressed in the current base, preceeded by a 
minus sign if X is negative, EXIT, 

PRINT (LPAR), 

PRINLIsr (X), 

X, 

Q1E3F0N; 

FONCnCN PRINLIST (X), 

PRINT (FIRST (X)), 

WHEN EMPTX (REST (X)), PRINT (REAR) DOT, 

PRINT (" "), 

WHEN ATOM (REST (X)), 

PRINT (". "), 

PRINT (REST (X)), 

PRINT (RPAR) EXIT, 

PRINLIST (REST (X)), 

ENDFDN; 

Interpretation: Print the stand 2 u:d list notation 
of the object X tO'the current output file. 


2. FUNCTION NES€JNE (), 

output a carriage return and line feed to the 
current output file, ti^ return FAL^, 

Interpretation: Terminate the current output line. 


3. FONCnCN PRINTLINE (X), 

PRINT (X), 

NEWLINE 0, 

X, 

E2I3FUN; 

Interpretati(»i: Print the expression X, terminate the 
last line and return X. 


4. FUNCTION SPACES (X), 

WHIN X > 0 AND X < 256, 

PRINT (" "), 

SPACES (X-1) EXIT, 
return the current cursor position, 
ENDFUN; 
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Interpretation: Output X spaces to the current output 
file and return the resulting cursor position. 

5. FUNCTION FRn-JAIH (EXl, FBP, UBP, PKTSPACE), 

Taking account of declared binding powers and any 
special princ rules on the pr<^rty list of PRn-iAIE, 
print a deparsed representation of BXl, assuming 
a) the curator to its left, if any, has right 

binding power RBP, the operator to its right, if 
any, has left binding power LBP, 
b} appropriate ^ces are to be printed if PRISFACS 
is nonFALSE. 

ENI^TJN; 

Interpretation: PREMAIS (expr, PBP, I£P) is a function 
which prints expr in standard mathematical form 
and surrounds it within parenthesis if the leading 
operator in expr has a left binding pwer less 
than or equal to RBP, or a right binding power less 
than ££P. 


6. WBS: FALSE; 

FUNCnCN \ms (X, Y, Z) , Ofrite Select) 

WHfN NCT B-im (WES),, . 

write out the ^inal record of WES and 
close the file, 

WES: FALSE, 

VJES (X, Y, Z) EXIT, 

WH3r EI-IPTY (X), \WSi FALSE EXIT, 

WHEN NAi-E (X) AND NAIE (Y), 

^flHEN BETY (Z), 

on the currently logged disk drive, 
delete any previous file named X.Y and 
make a new directory entry for X.Y, 

WES: X EXIT, 

WHIN NAIE (Z), 

on drive Z, delete any previous file 
named X.Y and make a new directory entry 
for X.Y, 

WES: X EXIT EXIT, 

a®FUN; 


7. FUNCTiaiLINELB331H(X), 

WHEN X > U AND X < 256, 

set maximum line-lengtn to X, 
return the previous line-length EXIT, 
return the current line-leigth, 

E5DFUN; 
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Interpretatioi: Set the length at which output lines 
will autonatically be tesninated. The line-*leigth 
is initially set to 72. 

8. PUNCTIGN RRDIX (X) , 

WHE2I X > 1 AIX> X < 37, 
set base to x, 
return the old base EXIT, 
return the current radix base, 

EX3DFUN; 

Interpretation: Set the base in which nunbers 
are expressed for both ii^t and output. The 
base is initially set to ten. 


Notes: 

1. Normally control of the current output file is done 
through the use of the function WBS as described above. However, after 
a file has been opened for output, output can be directed to the console 
without closing the disk file by simply setting the value of WBS to 
FALSE. A subsequent non>FALSE assignment to WBS will then redirect 
output to the disk file and append data onto the end of the file. 

2. If there is insufficient disk space or a hardware write 
error prevents correctly writing output onto the disk, an error message 
is sent to the console, the console is made the current output file, and 
an error-options trap occurs. 
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N. Evaluation Functions 


1. EROPERiy PREFIX, LIST (READLIST (SCAN)) & 

SUBROUTINE ' (X), 

X, 

ENDSUB; 

lnterpretatic»i: Stress evaluation and return the 
object X itself." 


2. FOICTCN EVAL (X), 

WHENAKHCX), FIRST (X) EXIT, 
raEN NAI'ffi (FIRST (X)), 

WHiH UNDEFnCD (GEIF (FIRST (X))), 

WHEN BQ (FIRST (X), EVAL (FIRST (X))), 

EVLIS (X) EXEC, 

EVAL (ADaon'J(EVAL(FIRST(X)), REST(X))) EXEC, 
I'JHIN FUNCTICNP (GEIF (FIRST (X))), 

APPLy (FIRST (X), EVLIS (REST (X))) EXIT, 
WHEN SUBROUTINEP (GETF (FIRST (X))) , 

APPLY (FIRST (X), REST (X)) EXIT, 

EVLIS (X) EXIT, 

WHEN FUNCTICNP (FIRST (X)), 

APPLY (FIRST (X), EVLIS (REST (X))) E5IT, 

WHEN SUBROUTINEP (FIRST (X)), 

apply (FIRST (X),'REST (X)) EXIT, 

EVLIS (X), 

ENDFUN; 

Interpretation: Evaluate tne object X. 

FUNCnCN EVLIS (X), 

WHEN ATOM (X), FALSE EXIT, 

ADJOIN (EVAL (FIRST (X)), EVLIS (REST (X))), 

ENKUN; 

FUNCTION GETF (X), 

contents o£ the function cell of the name X, 

EIDFUN; 

FUNCTION UNDEFINEI) (X), 

return FALSE if X is a pointer to tne undefined 
function trap, TRUE otherwise, 

ENDFUN; 

Interpretation: The recognizer for undefined functicns. 

FUNCncai FUl'JCTICNP (X), 

SUBR (X) OR EXPR (X), 

ENDFUN; 

Interpretation: The recognizer for call by value 
functions. 
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FONCnCN SBIOTTINEP (X), 

FSIBR (X) OR FEXPR (X), 

Q®FUN; 

Interpretation: The recognizer for call oy name 
functions. 

FONCnCN SDBR (X), 

return TRUE if X is a pointer to a FUNCTICN subroutine, 
FALSE otherwise, 

QIDFUN; 

FUNCnCSJ FSUBR (X), 

return IRUE if X points to a SIBPOUTZNE subroutine, 
FALSE otherwise, 

ENDFUN; 

FUNCnCN DCPR (X), 

FIRST(X) » 'D£PR, 

QDFUH; 

FONCnCN FDCPR (X), 

FIRST(X) ■ 'FEXPR, 

Q3DFUN; 


3. FUNCnCN APPLY (X, Y) , 

WBQI NAME (X), * 

WHIN UNDEFINED (GETF (X)), 

WHIN X » EVAL(X), FALSE EXIT, 

EVAL (ADJOIN (EVAL (X), Y)), 

WHIN SUBR (GEIF (X)), 

WHEN ATOM (Y), X (Y, FALSE, FALSE) EXIT, 
WHEN ATOM (REST (Y)), 

X (FIRST (Y), REST (Y), FALSE) EXIT, 
WHEN ATOM (RREST (Y)), 

X (FIRST(X), SECCND(X), RREST(Y)) EXIT, 
X (FIRST (Y), SECOND (Y), THIRD (Y)) EXIT, 
WHEN FSUBR (GEIF (X)), X (Y) EXIT, 

WHEN EXPR (GEIF (X)) OR FEXPR (GETF (X)), 

BIND (SECOND (GEIF (X)), Y), 

Y; EVAIBCDY (FALSE, RREST (GETF (X))), 
UNBIM) (SECOND (GEIF (X))), 

Y EXIT EXIT, 

WHEN EXPR (X) OR FEXPR (X), 

BIND (SECCtD (X), Y), 

Y: EVALBCI3Y (FALSE, RREST (X)), 

UIBIND (SBCCND (X)), 

Y EXIT, 

ENDFUN; 

Interpretation: A^ply the function X to the list 
of arguments Y. 
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FUNCTiaJ EVMSXiy (X, Y), 
mm ATOM (Y), X EXIT, 

1®EN ATOM (FIEST (Y)) OR ATOM (FIRST (FIRST (Y))), 

evaibchy (eval (first (y) ), rest (y) ) exit, 

'VfflEK ATOM (FIRST (FIRST (FIRST (Y)) )), 

X; EVAL (FIRST (FIRST (Y))), 

WHEN NOT X, EVALBCDY (X, REST (Y)) EXIT, 
EVAIBCDY (X, RE35T (FIRST (Y) )) EXIT, 

EVAIBCDY (EVALBCDY (X, FIRST (Y)), REST (Y)) , 
ENDFUN; 

ETOOICN BIND (X, Y), 

WHEN ATOM (Y), 

WHEN ATOM (X) , 

l^HEN eiPTY (X), FALSE EXIT, 

ARGSTACK: ADJOIN (EVAL (X), ARGSTACK), 
ASSIGN (X, FALSE) EXIT, 

ARGSTACK: ADJOIN (EVAL (FIRST (X)), ARGSTACK), 
ASSIGN (FIRST (X), FALSE), 

BIND (REST (X), Y) EXIT, 

WHEN ATa4 (X), 

WHEN n-IETY (X), FALSE EXIT, 

ARGSTACK: ADJOIN (EVAL (X), ARGSTACK), 

ASSIGN (X, Y) EXIT, 

ARGSTACK: ADJOIN (EVAL (FIRST (X)), ARGSTACO , 
ASSIGN (FIRST (X), FIRST (Y)), 

BIND (REST (X), REST (Y)), 

ENDFUN; 

FUNCnCN UNBD® (X) 

WHEN Aral (X), 

WHEN EMPTY (X), FALSE EXIT, 

ASSIGN (X, FIRST (ARGSTACK)), 

ARGSTACK: REST (ARGSIACK) EXIT, 

UNBIND (REST (X)), 

ASSIGN (FIRST (X), FIRST (ARGSTACK)), 

ARGSTACK: REST (ARGSTAa), 

ENDFUN; 


4. SUBRCUrnC CCND (XI, X2, ..., xn), 

EVALCCUD (LIST (XI, X2, ..., Xn)), 

EM)SB; 

FONCnCN EVALCCND (X, Y), 

WHENA3EM(X), FALSE EXIT, 

Y: EVAL (FIRST (FIRST (X))), 

WHEN NOT Y, EVALCCK3 (REST (X)) EXIT, 

EVALBCDY (Y, REST (FIRST (X))), 

SOFUN; 

Interpretation: Successively evaluate tne FIRST of 
XI, X2, ..., Xn until either a non-FALSE value is 
encountered or all have evaluated to FALSE. In the 
fomer case the REST of that argument is evaluated 
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as a function body (see the interpretation of 
AFPL2f for details). In the latter case FALSE is 
retunied by CCND. 


5. PRDPERIY PREFIX, LOOP, ADJOIN ('LOOP, HATCH(ENDLOOP)) ; 

SOBRCUTINE LOOP (XI, X2, ..., Xn), 

EVALLOOP (LIST (XI, X2, ..., Xn), 

LIST (XI, X2, ..., Xn)), 

BOSDB; 

FUNCTION E';ALLOOP (X, y, Z), 

WHEN ATOM (Y), EVALLOOP (X, X) EXIT, 

WHEN ATOM (FIRST (Y)) OR ATOM (FIRST (FIRST (Y))), 

EVAL (FIRST (Y)), 

EVALLOOP (X, REST (Y)) EXIT, 

WHEN ATOM (FIRST (FIRST (FIRST (Y)))), 

Z: EVAL (FIRST (FIRST (Y))), 

WHEN NOT Z, EVALLOOP (X, REST (Y)) EXIT, 

EVALBCDY (Z, REST (FIRST (Y))) EXIT, 

EVAI£a3Y (FALSE, FIRST (Y)), 

EVALLOOP (X, REST (Y)), 

ENDFUN; 

Interpretation: The LOOP construct evaluates its 
argument in a manner ideitical to the evaluatioi 
of the clauses in a function body. However, 
if all the arguments are evaluated without a 
conditional having been satisfied, evaluation 
begins again with the first argum^t. 

LOOP is the leading keyword of a control construct 
having the fom: 

LOOP 

taskl, 

- task2, 

• • • 
taskn 
ENDLOOP 

This construct parses to the internal representatioi (LOOP taskl 
task2 ... ). Usually at least one of the tasks is a conditional 
exit. Evaluation repetitively cycles through the sequence of 
tasks until a conditional exit causes control to proceed 
directly to the point following the matching delimiter ENDLOCP. 
The value of a LOOP construct is that of the last task evaluated 
therein. Since the LOOP construct parses to a function 
invocation, this construct can be used outside function 
definitions. 
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6. PBOPERiy PBEFIX, WHEN, MftTCH (EXIT); 

WHEN is the leading keyword of the conditional'exit 
control construct/ which has the general form 

WHQI expression!, expression!, ... SUT 

This construct parses to the internal representation 

((expression! expression! ... ) 

If expression! evaluates to FALSE, then evaluation proceeds 
directly to the point immediately following the matching EXIT. 
Otherwise, the expressions between expression! and the matching 
EXIT, if any, are successively evaluateo, and the last 
evaluated, after which evaluation proceeds to the point 
immediately following the next delimiter ENDLOCP, ENDBLOCK, 
EIDFUN, or QIDSIB. 


7. PK3PEm K?EFIX, BLOCK, I^RTCH (ENDBLOCK) ; 

BLOCK is the leading ke^rd for the control construct 
of the form: 

BLOCK 

WHEN ... EXIT, 

ENffiLOCK,-. 

As indicated, the first task within a block must be a 
conditional exit. Since other tasks within the block can also 
. be conditional exits, blocks provide a generedization of the 
"case” construct of some other languages, which includes the 
*if~then-else* constn;tct as a special instance. The evaluation 
of tasks within a block proceeds sequentially unless a 
conditional exit therein causes evaluati^ to proceed directly 
to the point following the matching delimiter Et^BLOCK. The 
value of a block is that of the last expression evaluated 
therein. 


8. FUNCnCN DRIVER (EXl, EX!), 

RDS: EXl, 

WRS: FALSE, 

NEI'^LINE 0, 

NEI JLINE 0 , 

LCQP 

ERR: FALSE, 

BLOCK 

WHEN ECHO {), 

PRINT ("? "), 

WHEN HOT RDS, PROTT {"") EXIT EXIT, 
ENDBLOCK, 

EXl; FALSE, 
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EXl: PARSE (SCAN(), 0), 

DC2: SCAN, 

n y/yK 

WEEN BCSO 0, NEWLINE () EXIT, 

EWDBLOOC, 

BTi OC ff 

WHEN EEOl OR NCT TEI^MINATQR (}, 

SmAX 0, 

NEWLINE (} EXIT, 

WHEiJ EX2 * •$, 

tANS: EVAL (EXl), 

WHEN BOiO (}, NEI^JNE (} EXIT EXIT, 
PRINT (@), 

*ANS: EVAL (Sa), 

PRIOT (" "), 
fffiEN BQ - 

PR3MA2H ("ANS, 0, 0, TRUE), 

NB€.INE 0, NEWLINE (), NEWLINE () EXIT, 
PRINTLINE (#ANS), 

NEWLINE 0, NB^INE (), 

ENDBLOCK, 

ENDLOOP, 

ENDFUN; 


DRIVER is a function which controls the interaction cycle. 
After establishing the console as the current input and output 
file by setting RDS and WRS to FALSE respectively, the main 
read, evaluate, and print driver loc^ is entered. An expression 
is first read by PARSE. If the terminator was the character 
the result is printed in mathenatical notatioi by PRn-lATE. 
If an *&", it is printed in List notation. And if a it is 
not printed at all. However, in all cases it is assigned to the 
variable IANS unless an error occurred during the parse ^aase in 
which case an error message is displayed. For some 
applications, it may be desirable to (perhaps dynamically) 
replace this driver with another one or to recursively call 
DRIVER. 



O 
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P. storage Functions 


1. FONCTIOl RECLAn4 (), 

reclaim all wreferenced noaes by generating a 
free node list from tnem, and ccnpact tne atom space 
and vector space. Ttie total resulting number of free 
nodes is returned. 

ENDFUN; 


2. CCNDQiSE: FALSE; 

FUNCTICN CCNDENSE (X, Y, 2), 

WHQl A!Ea4 (FIRST (X)) , 

HHEM ATOM (REST (X)} r FALSE EXIT, 

Z: SUBEXPN (REST (X), Y), 

WHEN aiPTY (2), CCNDEKSE (REST (X), Y) EXIT, 
REPLACER (X, 2), 

FALSE EXIT, 

2; SUBEXPN (FIRST (X), Y), 

WHEN EMPTY (2), 

CONDENSE (FIRST (X) , Y), 

WHEN ATOM (REST (X)), FALSE EXIT, 

Y: ADJOIN (FIRST (X), Y; , 

2: SUBEXPN (REST (X), Y), 

mm EMPTY (2), CCNDEMSE (REST (X), Y) EXIT, 
REPLACER (X, 2), 

FALSE EXIT, '' 

REPLACEF (X, 2),' 

WHEN ATOM (REST (X)), FALSE EXIT, 

2: SUBEXPN (REST (X), Y), 

WHEN EMPTY (2), CCNDEMSE (REST (X) Y) EXIT, 
REPLACER (X, 2}, 

FALSE, 

EMDFUN; 

FUNCTION SUBEXPN (X, Y, 2), 

WHEN CQ'iPARE (X, Y}, 

WHEN X » Y, Y EXIT EXIT, 

2! SUBEXPN (X, FIRST (Y)), 

WHEN EJIPTY (2) , SUBEXPN (X, REST (Y)) EXIT, 

2, 

ENDFUN; 

FUNCnOJ Ca-!PARE (X, Y), 
t4HEN Aia<l (Y), EXIT, 

WHEN AT0I4 (X), FALSE EXIT, 

WHEN CaffiARE (FIRST (X), FIRST (Y)), 

COMPARE (REST (X), REST (Y)) EXIT EXIT, 

EMDFUN; 
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Q. ^stem Functions 


1. FUNCnOJ SAVE (X, Y), 

WHEN NOT E^IFTY (WES) , 

write out the final record of and 
close the file, 

WBS: FALSE, 

SAVE (X, Y) SOT, 

WHUI NAME (X) AID NAME (Y), 

WHO! EMPTY (Y), 

save a binary oeniocy ioiaige of the current 
suSIMP system as a file named X of 
type "SYS" OR the current drive, 

TRUE EXIT, 

save a binary mcnory Image of the current suSSiP 
system as a file named X of ^pe "SYS” on 
drive Y, 

TRUE BilT 

ENDFUN; 



Note: SYS files occu^ about 15 kilobytes less than the 
memory size for which the (grating ^stem is 
generated. 


2. FUNCnCN LCftD (X, Y), 

raiEM NACC (X) AND NAME (Y), 

WHEN EI-IPTY (Y),, 

load a memory image file named X of type "SYS" 
frcm the current disk, 

return control to the executive DRIVER loop EXIT, 
load a memory image file named X of type "SYS" 
from drive Y, 

returo oxatrol to the executive DRIVER loc^ EXIT, 

EM^W; 



Interpretation: Restore the muSIMP enviroraoent present 
at the tiine of the SAVE. 
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V. The niuSIMP79 PRA2T Parser 


A. Operators 

1« INFIX is a name on whose property list is stored 
expressions specifying bow to parse infix curators for which mere left 
and right bindiitg powers do not suffice. It is used for the assignment 
operator "i” since a checic is made on its left operand to make sure it 
is a name. Also the INFIX property is used for "(" to correctly parse 
function calls written using mathematical notation. The respective 
operator's left-hand operand is passed to the expression as the fluid 
name "EXl". 

2. PREFIX is a neune on whose property list is stored 
expressions sp^ifying how to parse prefix (^rators for which mere left 
and right binding powers do not suffice. The matchfix operatorsr which 
include WHEN, UOOB, WXK, FUNCTION, SUBPOUTINE, PBOPERTy, and "(" when 
used to delimit a functions argument list, are examples of the use of 
the PREFIX prcperty. 


6. Bindi:^ Powers 

1. ISP is a name on whose property list is the integer left 
binding powers of infix and postfix operators. When two (^rators are 
competing for an operand, the operator with higher binding power toward 
the operand obtains the operand. In case of a tie, the left operator 
obtains the operand, so that infix operators with the same left and 
right binding powers associate left, as is usually desired. 

2. RBP is a name on whose property list is the integer 
right-binding powers of infix and prefix operators, for use as described 
for LBP. 


C. Constants 

1. COMMA is a global constant having the value Because 
of conflicting parse properties associated with its use as a separator 
character, the name O^lMA should be used for the literal 

2. LPAR is a global constant having the value Because 
of conflicting parse properties associated with its use as a separator 
character, the name IFAR should be used for the literal 

3. RPAR is a global constant having the value ")". Because 
of conflicting parse properties associated with its use as a s^aracor 
character, the name RPAR snould be used for the literal ")". 
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D. Delimiters 


1. DELIMITER is a name which is initialized to the list 
(Em, ENDLOQP, EIlDBLOac, EbJOFUN, QmB, RPAR, COMMA). When a matchfix 
operator is establisi^, adjoining the matching delimiter to this list 
aiables the j^^ser to give more informative diagnostics hy recognizing 
when a delimiter is used out of place. Bowever, it is not necessary to 
adjoin delimiters to this list, and adjoining a delimiter has the effect 
of precluding its use out of context for other purposes. 

2. DELIMITER () is a predicate which returns FALSE if the 
curreit val ue o f the name SCAN is neitimtr a terminator nor an the list 
named DELIMITIR. 

3. MATCH (delim) is a fwction which parses zero or more 
expressions separated by commas and delimited by the value of its 
argument. MATCH returns a list of the parsed representations of these 
repressions. 

4. MATCHNOP (expr, deHm) is a function used to verify that 
a matching "delim” was found following the PARSE of an expression within 
delimiters. 


E. Parsing 

1. PARSE (expr, rbp) is a function used to read a muMATH 
expression and convert it to List notation according to various rules 
established by operators LBP' and RBP binding powers and/or PREFIX or 
INFIX property rules as described earlier. 

Errors specific to PARSE include a member of DELIMITER "USED AS AN 
INDETERMINATE", an infix operator "USED AS AN PREFIX OPERATOR", and a 
prefix or postfix operator "USED AS AN INFIX OPERATOR?'. 

2. SYNTAX exprs is a function which takes an arbitrary 
lumber of arguments. If the value of the Global variable ERR is FALSE, 
then the message "*** SYNTAX ERROR; " is printed followed by the 
arguments to SYNTAX separated by spaces. Unless input echoing from a 
file, the remainder of the expression is printed until a TERMINATOR 
character is reached. Finally, in order to return control to the 
console, the control variable BDS (i.e. ReaD Select described in Sectioi 
IV. L.) is set to FALSE. 
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mTRix Package Docunisntatign 


PURPOSE; 

File MAIRIXJyRR provides the following matrix operations on arrays: 
transpose, multiplication, division, inverse, and other integer 
powers. Elementwise operations such as additicxi are provided by 
the prerequisite file ARRAYJ^. 

PREREQUISITE PILE; ARRAy.ARI 

US?CE; 

EDHAT (positiveinteger), 
array ' , 
arrayl . array2, 
arrayl \ array2, 
array “ integer 

E5CAMPLES: 

IDMftT(2) —> {[1], 

[ 0 , 1 ]}, 

If A* {[1, 2], and B = {P, then: 

[0, 3]} 

B' —> [P, 61, 

A' —> [{1, 

2 }, { 0 , 

3}], 

A' . IDI'lfiT(2) —> {[1, 0] , 

: [2, 3]}, 

A . B —> {P+12, 

18 }, 

A \ B —> {P-4, 

2 }, 

A - 2 —> {[1, 8], 

[0, 9]}, 

A —> {[1, -2/31, 

[0,1/31} 

RH’IARKS: 

1. The function named IDMAT returns a (left-triangular) 
identity matrix with the number of rows indicateo by its positive 
integer argumeit. 
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2. The postxix operator named having a left binding power 
of 160 r the same ais ”1"/ requests the transpose of its operand. 
(This ‘‘backward accent” characterr ASCII code 60 hex, different 
from an apostropne or single quote, is usually found on the same 
key as the character The trauispose of a scalar is a scalar, 

the transpose of a row is the column of the transposes of its 
elements, and the transpose of a column is the row of the 
transposes of its elenaits. These rules eire recursively employed 
so that the transpose of a ragged wd/or nested matrix is 
e^ropriately performed. These rules also caivert a column of rows 
into a row of columns, which does not print attractively. However, 
multiplication by an appropriate sized identity matrix always 
yields the attractive column-of-rows form of a matrix. 


3. The matrix-product infix operator designated ty a period 
has left and right binding powers 120, the same as for The 
interpretations are: 


scalarl . scalar2 —> scalarl * scalar2, 

scalar . array —> scalar * array, 

array , scalar —> array * scalar, 

row , col —> row .col + row ,col + ..., 

1 . 1 2 2 


col . row —> {[col .row , col .row , ...], 

1112 
[col .row , col .row , ...], 

2 12 2 

]}/ 


rowA . rows —> [[rowA .rowB , rowA .rov® , ...], 

1112 
[rowA .rov® , rowA .rows , ...], 

2 12 2 

]], 


colA . colB —> {{colA .colB , 

1 1 

colA .colB , 

1 2 

•.. } , 

(colA .colB , 

2 1 

colA .colB , 
2 2 


Consistent with the interpretaticxi described in ARPAy.D0C, when a 
row and column are of unequal length, the shorter is treated as 
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having implied trailing zero elements when forming "row . col" . 
These interpretaticws of matrix proouct are recursively anplc^ed so 
that matrix products of nested and/or ragged arrays are appro- 
priatedly performed. 


4. For a matrix A: 

A"0 IDMAT (LEN3TH(A)-1), 

A * 1 —> Ar 

A * -1 —> A inverse, 

For integer n > 1: 

A * n —> A . (A " (n-l)), 

A * -n —> (A * -1) . (A * (n+D). 

When a matrix is singular, raising it to a negative power yields 
warning messages about divisions by zero, and the offending 
sube^qpressions are eicapsulated in a question-mark form according 
to the usual muMATH-79 computational error treatment. 


5. When a matrix A is square and nonsingular, then A\B is 
equivalent to (A * -1) . B. However, WE STRONGLY RECOMMEND using 
A'^ unless the inverse is of indep^ent interest or must be used 
many separate times, because A\B is more efficient and because, 
provided B is ccxisistent, A\B will yield a parameterized solution 
even when A is singular. In this case, the parameters are 
designated by the forms ARB(l)., ARB(2) , ..., starting with 1 when 
file MATRIXJ^ or SQLVE.EQ(^ is most recently loaded. 


6. Comments in file MATRIX.ARR indicate now to save space by 
omitting the matrix transpose, division, or power packages. 
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lcas£ PsgKaqe D<?g«nentatiQn 


PURPOSE: 

File TRACE.MUS provides a trace package to help debug prograias. 
PRBQOISITE FILE: M0S1MP79.CC1M 
USAGE: 

1. TRACE (namely naine2/ ,,,), 

2. UWTRM:e (namel, nanie2, ...). 

EXAMPLE: 

FUNCnCN mm (EXl, EX2), 

WHEN EMPTY {EX2), FALSE EXIT 
WHEN EXl « FIRST (EX2), TRUE EXIT 
MEMB (EXl, REST (EX2)) 

ENDFUN; 

TRACE (Mil®); 

MB© ('DOG, '(CAT, COW, DOG, PIG)); 

MB© (DOG, (CAT, COW, DOG, PIG)] IThis is ccnputer gaierated% 

MB© [DOG, (COW, DOG, PIG) ] 

MEl© (DOG, (DOG, PIG)] 

MES© * TRUE 
MB© * TRUE 
MB© » TRUE 
§ TRUE 

UNIRACE (MEMB) ; 

RB©RKS: 

1. The trace of a function during the execution of a program 
provides an invaluable ddsugging tool. 

2. Whenever a function is called it arguments are first 
evaluated and thei printed following the function name. 

3. After the function has been applied it's value is printed 
following the functicxi name. 

4. Indention is used to more easily pair corresponding calls 
and returns. 

5. The functicsi is restored to rotmal by UNIRACE. 
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Basic ALGEBRA Package Documentation 

FQBFQ6E: 

File ADSEBRAJ^ provides for the basic algebraic simplificatic» ^of 
expressions using the elementary operators ■*", “/"r ana "*"• 

Simplifications may be categorized as either automatic or user 
controlled by means of CCtnSCL VARIABLES. 


AUTOMATIC SIMPLIFICATICNS: 

1. Rational arithmetic is used to combine numerical operands, 
(see A RTTHjy g for a conplete descripti<xi) 

2. Identities and zeros are appropriately €^lied to en^ressions. 

0+X —> X; 1*Y —> Y; 0*2 —> 0; 

3. Sums and products are flattened and uniquely ordered to 
facilitate expression comparisons. 

X+(y+Z) —> X+Y+Z; Z*(Y*X) —> X*Y*Z; 

4. Similar tenns and prodacts are conbined. 

3*X + 2*X —> 5*X; X*5 / X*2 —> X"3; 

5. Powers of #I (i.e. the square root of -1) are reduced. 

♦I^ —> -#I; 


CaORQL VARIABLES: 

Ti:^ control variables described in this section enable ti^ mui4ATB 
user to have complete control over the rules used to simplify an 
expression. Boweverr they are rather difficult for the novice to master. 
Therefore the utility functions EXPAND, EXPD, and FCTR (described 
below) have been included in muMATH to make it easy to obtain the most 
common forms of an expression without the need to individually set 
control variables. We recommend these functions be used until more 
precise control of the control variables is required. 

1. NUMNUM controls the distribution (factoring) of factors in the 
NUMerator of an expression over (from) a sum in the NUMerator. 

Identity: A * (B+C) <—> A*B + A*C 


2. DQDEN controls the distribution (factoring) of factors in the 
OENominator of an e:^res8ion over {from) a sum in the OENominator. 

11 1 

Identity; - * —- <—> . . 

A B + C A*B + A*C 
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3. D£l]NCi4 controls distribution (factoring) of factors in fchg 
DENominator of an expression over (from) a sum in the NUMerator. 

1 B C 

Identity: - * (B4C) <—> — + — 

A A A 

4. NUI’IDEN controls the distribution (factoring) of factors in the 
NDI-ierator of an expression over (from) a sum in the DENominator. 

1 1 

Identity; A * ——- <—> .— 

B + C B/A + <yA 

5. BAEEXF controls the distribution (factoring) of the B^ of an 
expression over (from) the Exponent. 

Identity; A(B+C) <—> aB * A^ 

6. EXFBAS controls the distribution (factoring) of the Exponent of 
an expression over (from) the BASe. 

Identity; (A*B)C <—> aC * bC 

7. PWREXPD controls whether or not integer PoWeBs of sums are 
EXPanDed in numerators anc^'or denominators. 

8« ZEBOEXPT controls the use of the following identity which is 
valid for all A not equal to 0. 

Identity; . A^ —> 1 

9. ZEROBASE controls the use of the following identity which is 
only valid for positive A. 

Identity; OA —.> o 


USAGE: 

1. For ti:» first six of the above control variables, the icinds of 
factors, bases, or exponents which are distributed or factored from the 
expression can be precisely controlled by assigning a^ropriate values 
to the respective control variable. Positive integer values will cause 
distribution; whereas, negative values cause factoring. The exact type 
of e;^r^sion which will be distriixited or factored can be determined 
from the following table: 

Prime 'lype Examples 

2 Numerical expressions 4, **1/3, 5/7 

3 Other non-sums X, SIN(Y), Z*3 

5 Suns RfS, X"2-X, LN(X)+Z 

Therefore, if a control variable is a multiple of one or more of 
the above primes, then that type of expression will be distributed or 
factored in accordarue with that control variable's identity transform. 
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2. For example, since differences are internally represented as 
suns involving negative coefficients, evaluation of 


3 * X * (1+X) * (l-X) —> 

3 * X * (1+X) * (l-X) 
X * (3+3*X) * (l-X) 

3 * (X+X*2) * (l-X) 

3 * X * (1-X"2) 
(3*X+3*X*2) * (l-X) 

X * (3-3*X"2) 

3 * (X-X*3) 

3*X - 3*X^ 


if NUMNDM is 0, 
if liMUM is 2, 
if is 3, 
if NIMCM is 5, 
if I04NQH is 6, 
if 104^ is 10, 
if is 15, 
if NUMNUM is 30 . 


3. As another exanple, if DE^EN is 15, then 

Y / 3 / X / (1+X) / (l-X) —> Y / (3*(X-X"3)). 

4. AS another example, if DE3MJM is 6, then 

(X+3) / 3 / X —> 1/3 + 1/X. 

5. When PWREXFD is a positive integer multiple of 2, then 
multinomial expansion occurs in numerators. When FWHEXTO is a positive 
integer multiple of 3, then multinomial expansion occurs in 
denominators. Thus, when PWPEXID is 6, 

(1+X)^ / (1+X+Y)*2 —> 

(1+3*X+3*X*2+X*3) / (1+2*X+2*Y+2*X*Y+X"2+Y"2). 

6. The importance of becoming thoroughly familar with the use of 

PWPDCID, DENOQ}, and DE1M)M cannot be over-emphasizedi muHAlB- 

79 cannot read a user's mind, so these control variables are the major 
means of specifying which of the mai^ alternative transformaticns are 
desired at each stage in a dialog. 


7. The remaiiiing control variables are of less frequent concern, 
but changing their settings is occasionally crucial to acheiving a 
desired effect. Since they follow the same general scheme, they are 
easy to use after the more important control variables have been 
mastered. For example,- 

(3+x) / (1+X) —> 

1 / (3/(l+X) + X/(l+X)) if NDMDEN is 5, 

1 / (l/(l/3+X/3) + 1/{1/X+1)) if NDMDEN is 30. 


Thus, this transformation yields a kind of "continued-fraction" 
expansion. 

8. BASQCP is set in an analogous fasioi as follows: 

2 * (1+N) —> 2 * 2^N 

if BASDCP is a positive integer multiple of 2, 
X * (1+N) —> X * X^ 

if BASEXP is a positive Integer multiple of 3, 
(A+B) " (l+N) —> (A+B) * (A+B)*N 

if BASQCP is a positive integer aultiple of 5. 

■Rm opposite of these transformations is more often appropriate, 
and is acccnplished by setting BASEX? to be negative. 
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ITTILITX FCNCnCNS: 

1. EVAL (expc) returns the evaluated and simplified expression 
resulting from expr operated on under the current control variable 
environm&it. 

2. SOB (^rl, expc2, expr3) returns the es^ression which results 
from SQBstituting all occurrences of expr2 by exDr3 in exorl. 

3. EVSUB (exprl, expr2, expr3) is defined as 

EVAL (SUB (exprlr expr2, expr3)). 

4. NUM (expr) returns the NUIterator of expr. 

5. DEN (expr) returns the DENdminator of expc^ 

6. FLAGS 0 prints the current value of the system control 
variable. 

7. EXPAND (expr) evaluates expr to yield a fully expanded 
denominator distributed over the terms of a fully expanded numerator. 
The following temporary assignments are made: 

FWREXPD: 6; NUMDEN: 0; NUMNUM: DENDEN: DENNUM: BASEXP: EXPBAS: 30; 


8. EXPD (expr) evaluates expr to yield a fully expanded numerator 
over a fully expa^ed denominator. The following temporary assignments 
are made: 

PWHEXPD: 6; NUMDEN: 0; DENNUM: -30; NUI-MJH: DENDQI: BASEXP: EXPBAS: 30; 


9. FCTR (e^^r) evaluates expr to yield a seai-factored numerator 
over a semi-factored denominator. The following tmnporarv assignments 
are made: 

PWPEXFD: NUHDES]: 0; NU!4NUM: DENDE^: -6; DE^INUM: BASEXP: BCPBAS: -30; 
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CGNmX VKSI?ELE 


Control Initial Positive Negative 

Var. Value Transfocsation Transfosnation 


NOr-fiVM 

6 

a*(bk:) 

-»> A*B + A*C 

A*B + A*C 

««> A*(B4C) 

DENDEN 

2 

11 
' m it MMM 


1 

1 


11 
•• it mmmmrn 

A B4C 


A*B + A*C 

A*B + A*C 

A hK 

DEtlNUr^ 

6 

m lj^ 

A 

“> 

B C 
- + - 
A A 

B C 
- + - 
A A 

-> 

sn\0 

A 

NDt-lDEN 

0 

A 


1 

1 


A 

n+r 


B/A + C/A 

B/A + C/A 


UTrV 

BASE^CP 

-30 

A"(B+C) 

-> 

A*B*A"C 

A'^*A*C 

-> 

a*(&k:) 

DCPBAS 

30 

(A*B)"C 

-> 

A*C*B"C 

A*C*B"C 

-> 

{A*B) *C 


A 

(A+B) ^ 

—> A*l«-...46"N 

(A+B) *-N - 

-V _ 

1 

0 

. . "" ' ' " 

A*Nf...+B"N 
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Equation Package 




itatifln 


FORFOSE: 

File BQliAljG provides a facili^ wherety equations are treated as 
expressions which can be assign^, added, multiplied, squauted, etc. 

FREREQOISHE FILE: PUSEBStuASl 

USAGE: expression! « expression! 


QOMFLES: 

BQNl: 5*X - 3*X - 7 — 2 + 4; 

—> 2*X - 7 — 6 

then 

EQNl + (7 — 7); 

—> 2*X — 13 

thei 

#AHS/2; 

—> X — U/JL 


REMARKS: 

1. The two sides of the eqpaation are ind^iendeitly simplified 
according to the current control settings. However, there is no 
attempt to automatically shift terms from one side to the other, 
etc. Moreover, there is no attempt to verify or disprove that the 
equaticxi is an id^tity or has a soluti^ 

2. This use of the ■* sign to indicate equations should not 
be C(mifused with the use of « within the conditional EXIT construct 
in muSIMF functi^ definitions. When used in this more active role 
the result is always either TRUE oc FALSE d^ending i;^n wi^ther or 
not the left and right sides have identical (as distinct from 
equal) values. 


3. The left and right binding powers of are 80, which is 
the same as for >. 

4, As Illustrated by the above example, whei a non-equaticxi 
is combined with an equation, the non-equation is independently 
combined with both sides. 


5. Although the above example illustrates how equations can 
be solved st^ise, file SCLVE.EQN autcxnates this process. 

6. Frovided file ARRAY.ARI is loaded, sets of simultaneous 
equations can be represented as an array of equations. For 

[2*X — 6 , 4*y —81/2; —> [X — 3, 2*Y — 4] . 

7. As with manueQ computation, operations such as squaring 
both sides or clearing non-numeric denominators can enlarge the 
solution set, so the user should exercise caution and verify 
ca nd i d ate solutions generated by such means, 

8. If FOO is an ecjiation, then SBtXND (FOS) returns the left- 
hand side of the equation, and THIRD (POO) returns the right-hand 
side. 
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Solve Bouation Package Documentation 


PURPOSE: 

File SOLVE.EQN provides a function for the exact solution of an 
algebraic equation. 

FREREQOISIIE FILE: BONJXX 

USAGE: SOLVE (equatioor unkncwn) 

ESSAMPLES: 

SaVE {X"2 — 4*A, X); —> {X — 2*A*(V2) 

X —-2*A*(1/2)} 

SOLVE (Ill(AmN(X-l)) — B, X); —> {X — 1 + TAN(#irB)} 


REMARKS: 

1. SOLVE returns a column of solutions, where columns are as 
described in file ARRAYJ}OC. The functions FIRST, REST, SECOND, 
and IBIRD can be used to extract individual solutions fraa a column 
of solutions. Alternatively, subscripts can be used for this 
purpose provided file ARRAYJ^ is loaded. 

2. Forgetting the second argument of SOLVE is a frequent 
mistake. 

3. As a c^ivanience, whai either side of an equation is zero, 
the 0 can be omitted. 

4. Wbai no soluticm exists, SOLVE returns the mpty column. 


5. When degenerate equations have an entire locus of 
solutions which require parameterization to r^resent coo^letely, 
SOLVE introduces the parameters, 

ARBd), ARB(2), ARB(3}, ... 

Their indexes start at 1 every time SOLVE.EQN or MATRIX.ARR is 
loaded. The following is an example of the iu:SIMP-79 solution to a 
degoierate station: 

SOLVE (X — X, X); —> {X-« ARB(1)}. 

6. SOLVE expands the difference in two sides of an equation 
over a common derarainator, then multiplies by the deiominator to 
clear it. This multiplication can introduce spurious solutiois if 
a zero of the denominator coincides with one of the numerator. 
Similarly, this multiplication can suppress a solution associated 
with an infini^ of the denominator. 'n»js, the returned set stould 
be regarded as caiviidates for some of the solutions rather than the 
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complete verified solution set if an elation has a denominator 
which could be zero or infinite for finite values of the unknown. 
When these possibilities are present, it is the user's 
re^cnsibility to verify his solutions fy substitution or perhaps 
by t£iking limits. It may be helpful in such instances to also use 
SOLVE to find any zeros of the common dmominator in order to see 
if they coincide with any of those in the numerator. 

7. After clearing the denominator, SOLVE attempts moderate 
factorization, then ind^^peidently attenpts to determine ti» zeros 
of each resulting factor. SOLVE r^ursively employs appr<^riate 
formulas for the inverses of the eleaentary functicms and for the 
zeros of linear, quadratic, and binomial factors. When SOLVE 
axxHmters a factor which it cannot treat, it returns a "solution” 
of the form "factor ■■ 0”. Since the factor may be simpler than 
the original equation, it might serve as a useful point of 
^parture for an a^r^imate numerical solution. 

8. A careful study of the source listing for the file 
SGLVE.EQN reveals bow additional inverse functions can be employed. 

9. File MA33UXARR contains a matrix division (^ratioi which 
can be used to solve simultaneous linear alg^raic equations. 
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Array Package 




itation 


File AHBMJVRl provides a facility for establishing generalized 
arraysr for extracting their components, and for performing 
elementwise operations between arrays or between arrays and 
scalars. 

PREREQUISITE FILE; ARITH.i'lUS 
USAGE: 

1. Fozmation of a ralumn vector: 

{e:qpressionl, e:^ression2r expressionN}. 

2. Fonnation of a row vector: 

[expression!, expression2r expression!!]. 

3. Extraction of coRiponents: 

array rowvector 

4. (^rations having fonus such as: 

arrayl qperator array2, 
scalar curator array, 
functionname (array). 


EXAl-IPLES: 


[0, X] + [5, X, y]; 

— > 

(5, 2"X, y] 

2 * {X, LN(y)}; 


r2"x, 

2*LN(y)} 

siN ([x, y]),- 

—> 

(siN(X), sii!(y)] 

IX, [y,Z], [W]][2J; 

—> 

[y^z] 

[X, [y,ZJ, [W]][2,l]; 

— > 

y 

[X, [y,Z], [Wll[2]UJ; 

— > 

y 


REMARKS; 

1. Arrays can be nested to any desired depth. The elements of 
a row or column can be any arbitrary expressions, including perhaps 
another row or column. 

2. Columns are printed starting each element on a new line. 
Thus, 2'>dimensional arrays generally look better as a column of 
rows than as a row of columns. Higher dimensional arrays generally 
appear best as a column of rows of rows ... of rows. 

3. When rows or columns of unequal length are combi.nec 
elementwise by an arithmetic operation, the shorter of the two 
arrays is treated as having implied zeros corresponding to the 
extra elements of the longer array. (Consistent with this 
interpretation, a subscript value larger tnam the numoer of 
solicit elements in a row or column yields a zero as the value of 
the element.) Thus, upper-triangular, left-triangular, and otner 
"ragged" arrays are efficiently represented. 
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4. When an array is combined with a scalatr the latter 
distributes over the elements of the surray. 

5. Functions of one argument such as SIN, ATANr etc., which 
employ the general rule-application function named SIMPU, 
distribute over the elements of an array. 

6. Subscripts can be recursively employed to any level, and 

ttey can be symboUa For example, [Y, Z] [2] [N] —> Z[N]. 

7. FIRST(row) —> [, FIRST(column) —> {, SEC0ND(row or 
column) —> first element, etc. 

8. Comments in file ARBAYJVRI indicate how to save space by 
omitting the column anc^of subscri^ packages. (Bows together witn 
FIBST, BEST, SECOND, etc. are sufficient for mar^ purposes.) 

9. File MATBIX.ABR implements matrix operations on arrays, 
including matrix transpose, multiplication, divisioi, and power, 
including inverse. 
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PURPOSE: 

File ARITH.MUS can perforin exact rational arithmetic operations 
including sumsr differencesr products^ quotients, and powers to 611 
digits of accuracy in any desired radix base, 

PREREQUISITE FILE: MUSIMP79.COM 

DCAt-IPLES: 

5/9 + 7/12; 

POO: (236 - 3*127) * -13; 

FOO * 16; 

RADIX (2); 

POO; 

1011000101 + 111010001; 

RADIX (1010); 

GO) (436 , 582); 


CCNTROL VARIABLES: 

1. PBRCH is a control variable which, when TRUE, permits 
selection of a branch of a multiplyrbranched function. For ARIULliUS, 
PBRCH nonFALSE permits the simpli|icati<xi 

(exprl * e^c2) * e;q»r3 —> exprl * (ex?r2 * expr3) 

evei when expr3 is not an integer. 

2. ZEROBAS is a control variable whichf when TIUJE, permits the 
simplification 0 * expr —> 1 even when expr is nonnumeric. 

3. ZEROEXP is a control variable which, when TRUE, ^rmits the 

simplification expr 0 —> 1 even when expr is nonnumeric. 


PRIMinVLF DEFINH) FUNCTICNS: 

1. ASS (expr) is a function which returns the absolute value of 
its argument when the argument is a rational number. Otherwise, the 
rule-application function SIMPU is invoked, so the unevaluated absolute- 
value form is returned if no applicable rules are present. 

2. ARSX (expr) is a helper function used ^ SII4PU and elsewhere 
to s^r(^riately partiti(xi an expression for application of a rule. 

3. ARGLIST (expr) is a helper function used to appropriately 
groi^ the operands of an expression for application of rules to varyary 
operators such as and 


% Exact rational arithmetic % 

% Make assignments to variables % 

% Raise nunbers to integer powers % 

% The radix base can be set from 2-36 % 
% Convert nunbers between radix bases % 
% Do binary arithmetic % 

% TO return to base 10 % 

% Ccnpute the GCD of two numbers % 
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4. BASE (expc) is a selector function which returns the base of 
an expressioi of the form base exp; otherwise it returns expr itself. 

5. QCDVf (expr) is a selector function which returns the 
codivisor (i.e. the non-numeric factors) of an expression which is a 
product; X i£ NUI1BER (expr); otherwise it returns expr itself. 

6. COEFF (expr) is a selector function which returns the 

coefficient (i.e. the numeric factors) of an expression which is a 
prodjct; the exor if (expr); otherwise it returns 1. Note that 

in all cases 


expr « CCEET (expr) * CCD2V (expr) 

7. DEN (expr) is a selector function which returns the 
denominator of its argument^ returning 1 when there is none. 

8. DENOH (expr) is a recognizer function which returns TIUIE iff 
its argument has the internal form C bas exp), with exp being negative 
or having a negative coefficient. 

9. EVSDB (expr, subexpr, r^lacement) is a function which returns 
the result of evaluating a copy of its first argument, wherein each 
syntactic occurrence of its second argument is replaced by the third 
argument. 

10. EXPON (expr) is a selector function which returns the 
exponent of an expression of the form base exp; otherwise it returns 
1. Note that in all cases 

expr » BASE (expr) * SCPCN (expr) 

11. GCD (intgrl, intgr2) is a function which returns the positive 
greatest common divisor of its integer arguments. 

12. IDENTITY (ej^r) returns its argument. This trivial functicxi 
is used for applying inverses and accomodating conditioial exits having 
atomic conditions. 

13. LCM (intgrl, intgr2) is a function which returns the positive 
least common multiple of its integer arguments. 

14. MIN (intgrl, intgr2) is a function which returns txm minimum 
of its two int^er arguments. 

15. MULTIPLE (intgrl, intgr2) is a function which returns FALSE 
if its second integer argument is not an integer multiple of its first 
integer argument. 

16. NS3COCEF (e:^r) is a recognizer function which returns TPUE 
iff its argument is negative or has a negative coefficient, returning 
FALSE otherwise. 


2 



17. NESMCLT (intgcl,. int9r2) is a predicate which returns TRDE 
iff its second integer argument is a negative integer multiple of its 
first integer argument. 

18. NDM (exgv) is a selector function which returns the numerator 
of its argument/ returioing the entire argument when there is no 
denominator. 

Id. NUMBER (expr) is a recognizer functioi which returns 'SSJE iff 
its argument is an integer or a rational number. 

20. POSMULT (intgrl/ intgr2} is a predicate which returns TRUE 
iff its first integer argument is a positive multiple of its second 
integer argument. 

21. FCW£R (expr) is a recognizer function which returns TRUE iff 
its argument is of the form exprl * expr2/ returning FALSE otherwise. 

22. PRODUCT (expr) is a recognizer function which returns TRUE 
iff its argument is of the form exprl * expr2/ returning FALSE 
otherwise. It is important to realize that quotients are represented as 
products involving negative powers. 

23. BBCIP (e:q>r) is a recognizer function which returns TRUE iff 
its argument is a rational number of the form l/d, returning FALSE 
otherwise. 

24. SIMPU (name, expr) is a function which applies any 
j^ropriate established rules for' the unary functioi or aerator whose 
name is the first argument of SIMPU and whose operand is the second 
argument of SIMPU. 

25. SUB (^cpTr subes^r, r^lacement) returns a copy of its first 
argument/ wherein every syntactic instance of its second argument is 
replaced by its third argument. In general this will produce an 
unsimplified result, so the similar E7SU3 functi(xi uses SUB, then EVAL. 

26. SUM (expr) is a recognizer function which returns TRUE iff 
its argumwt is of the form e:q>rl 9xpt2f returning FALSE otherwise. 
It is important to realize that differences are represented as sums 
involving terms having negative coefficients. 
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^ B3ffiR Package 


PORPC6E; 

Provides the facilities for the simplificatioi of fractiaxal powers 
of numbers and complex expoisntials. 

USAGE: 

nunijer * (fraction ), 

*E " (intgr * *1 * #PI / 2). 

EXAMPLES: 

(-24) ^ (1/3) -> -2 * 3 * (1/3), 

(-4) ^ (1/2) -> ♦! * 2, 

#E " (3 * #I * tPI / 2) -> - #1. 

(XmSCSL VARIABLE: 

PBRCHr which if FALSE, prevents Picking a BRanCH for fractional 
powers, (e.g. 4 * (1/2) will not simplify to 2.) 

REMARKS: 

1. #E represents the base of the natural logarithms, #I represents 
the positive ^are root of minus one (+ (-l)^(l/2)), and #PI represents 
the ratio of the circumference of a circle to its diameter. 

2. Simplification of fractioneU powers takes place only if the 
control variable named PBRCH is not FALSE. The positive real branch is 
selected if one exists. Otherwise, the negative real branch is selected 
if one exists. Otherwise, the branch with smallest positive argument is 
selected. 

3. As in manual computations, Picking a BRanCH of a fractional 
power involves an arbitrary choice which can yield invalid results. 
Thus, the user is cautioned to verify results obtained by such 
(^rations. 

4. The global variable named PRIMES oxitains a list of successive 
primes, beginning with the integer 2. For fractional powers, the 
radicand is factored into a product of powers of the ixmobers in HUMES, 
perhaps times a residual having no factors in PRIMES. The fractional 
power is then distributed over this prodxt, with a discrete variant of 
Newton's method being used to determine if tm fractional power of any 
residual is an integer. Thus, simplification of fractional powers of 
large integers might be ineanplete if PRIMES is not long enough. 

5. As in manual computations, reduction of complex e^^nentials 
modulo (2 * #PI * #1) is inconsistent with the identity LN(Z*W) « III(Z) 
* LN(W). Thus, the user is cautioned to verify results obtained using 
both transformations together. 
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National FACTORIAL Package 


PURPOSE: 

Provides the factorial |»stfix operator "I*. The factorial of a 
non-negative integer is recursively defined as follows: 

01 • 1 , 

N1 ■ N*(N-1) I, for N > 0. 


USAGE: 

N1 where N is a non-negative integer. 


DCAMPLE: 

5i —>120. 


REMARKS: 

1. The left binding power of ■!" is 160. Thus -51 parses to 
- (5!) aaid 3*51 parses to 3'‘(51). 

2. When not given a nonnegative integer operand, ”1” calls upon 
the SIMPU rule-application function, thus returning the unevaluated 
factoriad fonn if no apprc^riate rules are established. 
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FJRPOSEs 

File L0G.AL3 provides for logaxithaic siicplificatians. 

PBSRBQUISITE FILE: ALGEBE2V.ARI 
OMIROL Vl^RIABLSS: 

1. LOGBAS, which is the default LCXSarithro BASe when LOG is given 
only ^e argumait. 

2. FBBCHr which if FALSE prevents Picking a BRanCH of logaritxsns. 

3. lOGEXIDr which controls expansion or collection of logarithms, 
and base conversion. 


USAGE: 

LN (expr) , 

LOG (expr), 

LOG (e:qpr, base). 


Pa>IABKS: 

1. Since the emphasis of muJlATH is on exact results, there is no 
attempt to a^roximate irrational logarithms. 

2. The unbound variable represents the base of the natural 

logarithms. / 

3. Although all logarithms are stored internally as two argument 
functions, LN (expr) is used as an aobreviation for LOG (expr, »£} on 
input and ou^t. 

4. LOG (expr) is used as an ^reviation for LOG (expr, LOCBAS) on 
input and output, where LOGBAS is a control variable initially set to 
#E. 

5. base * LOG (ej^r, base) —> expt, 

6. Provided PBPCB is IKJE: 

LOG (1, base) —> 0, 

LOG (base, base) »> 1, 

LOG (base^expr, base) —> expr. 
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7. Provided LCGESffO is a positive integer multiple , of 2: 
LOG (expr,base) —> LOG (expr, #E) / LOG (base, #E) 


when base is not #E. When LOGEXPD is a negative integer multiple of 2, 
the opposite transformation of combining appropriate ratios of 
logarithms occurs. 


8. Provided LOGDCPD is a positive integer multiple of 3: 

LOG (esqpr^exp, base) -—> exp * LOG (expr, base). 

A negative integer multiple of 3 causes the epposite transformation. 


9. Provided LOGEXPD is a positive integer multiple of 5: 

LOG (exprl*expr2, base) —> LOG (exprl, base) LOG (expr2, base), 
LOG (exprl/expr2, base) —> LOG (exprl, base) - LOG (expr2, base). 

A negative integer multiple of 5 causes the opposite collection 
transformation. 
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Simplification Package Documentation 


PURPOSE: Pile TRGPOS.ALG provides the following trigonometric trans¬ 
formations: 

1. exploitation of symmetry to simplify trig arguments 

2. replacement of other trig functions 1:^ sines and cosines 

3. replacement of integer powers of sines and cosines by linear 
combinations of sines and cosines of multiple angles 

4. replacement of products of sines and cosines by linear 
combinations of sines and cosines of angle sums 

5. replacement of integer powers of sines by those of cosines or 
vice-versa 


PREREQUISITE FILES: ALSEBRUARI 

(Note: Loading TRGNEG.ALG after TRGPOS.ALG preserves the full 
capabilities of both files. Loading TRGPOSJUlXs after TBGNEG.AI/3 
destroys the angle-reduction capabilities of the latter, thus 
saving some space.) 


CCNIRQL VARIABLES: 

1. TRaEXPD controls replacement of trig functions by sines and 
cosines and replacsmient of powers and products of sines and cosines by 
linear combinations. Only positive values of IRGE2CFD are significant 
when TRGPOSJUiG is loaded without TRGNBGJUjG. 


2. 'HUSSQ controls the conversion of integer powers of sines to 
cosines and vice-versa. 


USIGE: 

SW (expression), 
COS (expressicxi), 
(expression), 
CSC (expression), 
SEC (expression) , 
COT (expressicxx). 


REIAHRS: 

1. SIN(O) —> 0, and COS(O) —> 1. 

2. Symmetry is exploited to simplify the arguments of sines and 
cosines. For excample, SIN(-X) *> -SIN(X) and COS(-X) —> COS{X). 
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3. When TRGESCTO is a positive multiple of 2, then tangents, 
cotangents, secants, and cosecants are replaced by corresponding 
expressions involving sines anl/or cosines. For example, when TBGEXH) » 
30, CSC(X) —> 1/SIN(X). 

4. When TEGEXPD is a positive multiple of 3, then integer powers 

of sines and cosines are expanded in terms of sines and cosines of 
multiple auigles. For example, when TBGEXPD * 30, COS(X)'‘2 ~> 
(l'KX6(2*X))/2. These transformations usually give the most attractive 
results if NCJMMJM and perhaps also eure positive multiples of 6. 

5. When TPGEXPD is a positive multiple of 5, then products of 
sines and cosines are expanded in terms of angle sums. For example, 
when TRGEXPD is 30, SIN(X)*S1N(Y) —> (COS(X-y) - C0S(X+Y))/2. These 
transformations usually give the most attractive results if NUMIAJH is a 
positive multiple of 30 and DEltUM is a positive multiple of 2. 

6. Expanding over a common denominator with TFGDCFO * 30 yields a 
normal form for a large class of trigonometric-*rational expressions. 
Thus, the most straightforward way to prove most trig identities is to 
evaluate the difference in the two sides with TBGDCPD: NUMHJM: 

30, FWBEXED: €, and OEtJNUM: >30. 

7. TRGEXPD » 30 has the effect of "linearizing" trigonometric 
polynomials, thus facilitating harmonic or Fourier analysis. 


8. For integer n with |n'| >1 and for all u, when TPGSQ is a 
positive integer, then 

COS(u)*n —> COS(u)"REI'!AINDER(n,2) * (1 - SIN(u)*QUC7nEOT(n,2))"2. 

Conversely, when TPGSQ is a negative integer, then 

SIN(u)"n —> SIN(u)"REMAINDER(n,2) * (1 - COS(u)"QUCTIEWr(n,2))*2. 

These transformations are sometimes useful for transforming a 
trigcxiometric polynomial to a more compact equivalent trig polynomial. 


9. Even when a trig polynomial is preferred for the final form, 
net simplification is often achieved by evaluating with ^n^GEXFD • 30, 
then -30, then perhaps again with 138330 • 1 or -l according to the 
appearance of the result produced by -30. 

10. File TR3NBG.AIiG provides for the negative settings of TPGEXPD 
to yield the converse of the above transfonnatiozxs. 


2 



TRGNEG.DOC (c) 


02/09/80 


The Soft Warehouse 


HsgatJLYg 2SISC 




RIC Simplification Package Documentation 


PUBPOSE: File TBGNEGJUiG provides the following trigoiometric trans- 
fOQoations: 

1. exploitation of syniaetries to sini>lify trig argvaneits 

2. angle reduction 

2, sultiple-angle expansicxi 

4. angle>sum expansion 

5. elminaticn of reciprocals of trigcnonetric forms 

6. eliminatiox of certeun products of trigoncmetric fosns 

7. simplificatioi of trig functiois of their own inverses 

8. replacement of sines and cosines by complex exponentials 


PBEBEQUISZTE FILE: ALSIBBA.ARI 

(Note: Loading TRGPOS.ALG after TRGNEGJ^LG destroys the angle- 
reduction capabilities of the latter, thus saving some s^ce. 
L oa di n g TBGNEG.ALG after TEX30S.ALS preserves the full capabilities 
of both files.) 


CCNEROL VARIABLES: 

1. TRGEXPD controls the use of multiple angle and angle sum 
expansions and r^aceaent of trig functions by complex exponentials. 
Only negative values of TI^EXPD are significant when TRiGN£G.ALG is 
loaded without TRGPOS.ALG. 


USAGE: 

SIN (expression), 

COS (e^gpressi^), 

TAN (expression), 

CSC (expressicxi), 

SBC (expression), 

COT (expression), 

mSEXH} (expression, integer). 


RQ1ARKS: 

1. Since the emphasis of muHATH-79 is on exact results, there is 
no attempt to a]^roximate irrational trig expressions. 

2. The ratio of the circumference to the diameter of a circle is 
repr^ented bY the unbound variable #PI. The user is of course free to 
assign a rational approximation to #PI and use series agproximations to 
the trig functions. 



3« Angles axe assumed to oe measured in radians. Those who would 
prefer some other unit such as degrees may wish to define additional 
functiois named SCOr OOBD, etc. 

4. Sines and cosines of aisles which are raomeric multiples of *PI 
eure reduced to ^ivalent sines or cosines in the range [Or tPV4lr then 
sines emd cosines of the special angles Or #PI/6r and tPI/4 are 
evaluated exactly. For example, 

SIN (20*#PI/7) —> SIN (tPI/7)r 

and SIN (7*#PI/3) —> 3*(l/2)/2. 

5. Symmetry is exploited to simplify the arguments of sines and 
cosines. For example, 

SIN (-X) —> -SIN (X), 

and COS (-X) —> COS (X). 

6. Trigonometric functions of the corresponding inverse trig 
functions simplify. For example, SIN(ASIN(X+5)) —> X+5. The inverse 
trig functions are named AXAN, ASIN, ACCS, ACCT, ACSC, and ASEC. 

7. Products of a tangent, cotangent, secant, or cosecant with 
another trig function of the same argumeit are simplified to 1 or to a 
single form where possible. For example, 

SEC(X)*COS(X) —> 1, 

and TAN(X)*COS(X) ~> SIN(X). 

For an expression such as SEC(X)''2*COS(X)''2 it is necessary to 
reevaluate with EXPBAS being a negative multiple of 2 in order to 
acheive the desired trig tranSformaticsi. 

8. When TPGEXPD is a negative multiple of 2, then negative powers 
of tangents, cotangents, secants, and cosecants are replaced by 
corresponding positive powers of the corresponding reciprocal trig 
functions. For example, when TK3EXPD * -6, lAAN{X+7)"3 —> CCT(X+7)*3. 
For techniced reasons, negative powers of sines and cosines 2 u:e treated 
in file TEGPOSJ^ 

9. When TRGEXPD is a negative multiple of 3, then sines and 
cosines of mulitple angles are expanded in terms of sines and cosines of 
non-multiple angles. For example, when TEGEXPD » -6, 

SIN (2*X) —> 2*SIN(X)*C0S(X) 

and COS {3*X) —> 4*C0S[X)^ - 3*C0S(X). 

•ntese transformations usually give the most attractive results if NOMNUM 
is a positive multiple of 6. 

10. When TBGEXPD is a negative multiple of 5, then sines and 
cosines of angle sums and differences are exp^ed in terms of sines and 
cosines of nonsums and nondifferences. For example, vtaox TPGEi!H>>-15, 
COS(X+y) —> COS(X)*COS(Y) - SIN(X)*SIN(Y). These transformations 
usually give the most attractive results if NUMNUM is a positive 
multiple of 6. 


11. When TBGExro is a POSITIVE multiple of 1, then sines and 
cosines are converted to complex exponentials. For example, when 
TPGEXPD - 14, then COS(X) —> (♦E"(«*X) + 1/#E"(#I*X)) / 2. The 
opposite transformation, provided in file ARITHjros, is requested when 
TRGEXPO is a negative multiple of 7. A worthwhile net trig 
simplification can sometimes be achieved by converting to complex 
expor»ntials, expanding or factoring judiciously, tt^. converting back 
to trig functions. 

12. In muMAIH-79 changing the value of an option variable does not 
affect the values of expressions which have already been evaluated. 
Thus, after changing the value of TSGSSSD and other relevant variables 
it may be necessary to use EVAL to get the desired effect. 

13. Function TBGEXPD reevaluates its first argument with IPGEXro 
temporarily set to the value of the second argument. Thus, it provides 
a coiveni^t way to accomplish a trigonometric transformation without 
the necessity of altering the global setting of the TBGEXH) control 
vauriable. 

14. File II^GPOSJ^ has other important trig transformaticns, many 
of which are the opposite of those provided in file TPGMEGJUJG. 
Generally, the positive settings yield a more canonical (but not 
necessarily more compact} representaticxv A net simplification is often 
achieved evaluating an expression with the relevant ^ion varieties 
set positive, then reevaluating with them set the other way. Thus, 
files TEGPC6JVDG and TBGNEG.AU3 emprise an important coroplenentary pair 
pair of files. Since together the files are relatively large, for some 
applications it may be desirable to extract and combine a few of the 
required features from both files, together perhaps with a few 
additional transformations modeled after them. 
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■tation 


PORPOSE; 

File DIF.AIi3 provides a function which returns the symbolic first 
partial derivative of its first argument with reject to its seowd 
argument. 

He^'ISITE FILE: ALGISRA.ARI 


OSAGE: 

DIF (expression, variable), 
DOM'fPLES: 

DIF (A*X"2, X) —> 2*A*X, 

DIF (m(X+A), X) —> 1/(X+A). 


REMARKS: 

1. When the differentiatiai rule for a function or operator is not 
known to the ^stas: 

a. The derivative is 0 if none of the arguments or operands 
contain the differentiation variable. For example, 

DIF (F(Y), X) —> 0. 

b. The derivative is not evaluated otherwise. For 
example, 

DIF {F(X), X) —> DIF (F(X), X). 


2. A careful study of file DIFJOLG reveals how additional 
differentiatioi rules can be inserted. 

3. The differentiation "variable" can actually be an arbitrary 
expression, which is then treated the same as a simple variable for 
differentiatioi purposes. (This is occasionally quite useful, such as 
when perfoming a square-free factorization or wh^ deriving idle Biler- 
Lagrange equations for a specific variational calculus problem.) 

4. Higher-order partial derivatives can be requested directly by 
nested use of DIF, such as DIF (DIF(SIN(X*Y),X), Y). However, beware 
that repeated differentiation can require dramatica^y increasing time 
and space, especially for products, quotients, and composite 
expressions. 


5. The useful utility function FREECF (exprl, expr2) is a 
predicate which returns TRUE iff exprl is free of (i.e. contains no 
occurrences of) e:^r2. 
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PURPOSE: 

Pile IMT.DIT provides facilities for indefimte syirbolic 
integration. 

PRERBQOISHE PILE; DIP J^ 

(mSE: 

INT (expression, variable}. 

EXAMPLE: 

INT {A*X + SIN(X) / X) —> A*X*2/2 - COS (X). 


REMARKS; 

1. When IWT is unable to determine a closed-form integral of 
portions of an expression, the returned expression will contain 
unevaluated integrals of ti^e porticns. Por example, 

INT (X + A*#E'VX, X) —> r2/2 + A*INr(#E'VX,X). 

2. BIT merely uses distributicn over sums, extractio: of factors 
which do not depend upon the integration variable, known integreds of 
the built-in functions, a few-reduction rules, and a "derivatives- 
divides" substitution rule. Ca::se^ently, integration succeeds cxily for 
a relatively modest class of integrands. Bowever: 

a. The class is large enough to be quite useful, 

b. Pile IMTMCRE.IKr contains additional rules, 

c. Integration of a truncated Taylor-series approximaticxi 
of an integrand can often yield a truncated series 
representation of otherwise intractable integrals. 

3. A careful study of files INTJ3IP and INTMORE.INT reveals how 
additional integrati^ rules can be inserted. 

4. The integration "variable” can be an arbitrary expression, 
which is then treated the same as a simple variable for integration 
purposes. 

5. Successful integration may depend upon the form of the 
integrand, after it is simplified according to the current flag 
settings. Generally sj^aking, it is best to employ conservative flag 
settings which do relatively little to adter the form of an expression. 
IKT will automatically expand, factor, employ trigonometric trams- 
formations, etc. as necessary. 
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SDRPOSE: 

File IKIIO&llir provides synitolic defuiite integration and extends 
the power of the indefinite integraticxi provided by file INIJOIF. 

PREREQUISITE PILE: lOT.DIF 

OSAGE: 

INI (expression, variable) , 

DEFINT (e:^cession, variable, Icwerlinit, bEperlisiit). 

EXA:4ELE: 

DEFINT (A*X"2, X, 0, 1) —> V3. 

RBARKS; 

1. OEFQir merely uses substitution into the indefinite integral, 
which is appropriate c^y for proper integrals. 

2. Whei DEFINT is unable to determine a closed-form integral, the 
unevaluated integral is returned. For exao^le, 

DEFINT (X+A^fE'VX, X, 0, 1) —•> DEFINT (X+A**E"X/X, X, 0, 1). 


3. Nested integration can be used to request directly an iterated 
integration, such as occurs for epprc^riate multiple-integrations. For 
example, to integrate the expression y*x*2 over the upper unit semi¬ 
disk, we could evaluate 

DEFINT (DEFINT(y*X*2,y,0,(l-X*2)“{V2)), X, -1, 1). 

However, beware that the class of expressions which is repeatedly 
integrable is dramatically smaller than the class which is once 
integrable. 

4. File INXJDOC contains other appropriate remarks. 
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The following is an index of all the important function^ 
variable, and constant names in both muSIMF and muMATE. Each name is 
followed by the module in which it occurs, a descriptor indicating the 
name's use, the page in the module's documentation on which it is 
explained, and finally the page in the module's muMATH source file on 
which it is defined. Function names are indicated by a set of 
parentheses following the name which contains the usual number of 
arguments given to the function. An asterisk {*) in the "Page 
Defined" column indicates that the item is incrementally defined in a 
number of places within the source. 


Name 


Desgxiptflc 

Page 

Documented 

Page 

Defined 

ABS (1) 

ARITH 

Numerical 

1 

3 

ADJOIN (2) 

muSIMP 

Constructor 

8 


AND (N) 

muSIMP 

Logical 

12 


APPEND (2) 

MATRIX 

Constructor 


2 

APPLY (2) 

muSIMP 

Evaluator 

. 28 


ABGEX (1) 

ARITH 

Selector 

1 

4 

ARGLIST (1) 

ARITH 

Selector 

1 

4 

ARRAY (1) 

ARRAY 

.Recognizer 


2 

ASSIGN (2) 

muSIMP 

Assignment 

13 


ATOM (1) 

muSIMP 

Recognizer 

10 


ATSOC (2) 

muSIMP 

Property 

14 


BASE (1) 

ARITH 

Selector 

2 


BASEXP 

ALGEBRA 

Control variable 

2 

4 

BASEXP (1) 

ALGEBRA 

Recognizer 


4 

BLOCK 

muSIMP 

Keyword 

31 


CODIV (1) 

ARITH 

Selector 

2 

3 

COEFF (1) 

ARITH 

Selector 

2 

4 

COL (1) 

ARRAY 

Recognizer 


2 

COMMA 

muSIMP 

Constant 

35 


COMPRESS (1) 

muSIMP 

Sub-atomic 

18 


CONCATEN (2) 

muSIMP 

Modifier 

9 


COND (N) 

muSIMP 

Evaluator 

29 


CONDENSE (2) 

muSIMP 

Storage 

33 


COS (1) 

TRGPOS 

Numerical 

1 

1 

COS (1) 

TRGNEG 

Numerical 

2 

1 

COT (1) 

TRGPOS 

Numerical 

1 

1 

CSC (1) 

TRGPOS 

Numerical 

1 

1 

DEFINT (4) 

INTMORE 

Numerical 

1 

1 

DELIMITER 

muSIMP 

Constant 

36 

22 

DELIMITER (1) 

muSIMP 

Recognizer 

36 

22 

DEN (1) 

ARITH 

Selector 

2 

2 

DENDEN 

ALGEBRA 

Control variable 

1 
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Name 

ttQdule 



DENDEN(l) 

ALGEBRA 

Recognizer 


3 

OENNUM 

ALGEBRA 

Control variable 

2 

iO- 

DENNUM (1) 

ALGEBRA 

Recognizer 


DENOM (1) 

ARITH 

Selector 

2 

4 

DIF (2) 

DIF 

Numerical 

1 

1 

DIFFERENCE (2) 

muSIMP 

Numerical 

19 


DIVIDE (2) 

muSIMP 

Numerical 

19 


DRIVER (0) 

muSIHP 

Evaluator 

31 


ECHO 

muSIMP 

Control variable 

22 


EMPTY (1) 

muSIMP 

Recognizer 

10 


ENDBLOCK 

muSIMP 

Delimiter 

31 


ENDFUN 

muSIMP 

Delimiter 

17 


ENDLOOP 

muSIMP 

Delimiter 

30 


ENDSUB 

muSIMP 

Delimiter 

17 


EQ (2) 

muSIMP 

Comparator 

11 


EVAL (1) 

muSIMP 

Evaluator 

27 


EVSDB (3) 

ARITH 

Constructor/Evaluator 

2 

1 

EXIT 

muSIMP 

Delimiter 

31 


EXPAND (1) 

ALGEBRA 

Evaluator 

4 

5 

EXPBAB 

ALGEBRA 

Control variable 

2 

4 

EXPBAS (1) 

ALGEBRA 

Recognizer 


4 

EXPD (1) 

ALGEBRA 

Evaluator 

4 

7 

EXPLODE (1) 

muSIMP 

Sub-atomic 

18 


EXPON (1) 

ARITH 

Selector 

2 

4 

FCTR (1) 

ALGEBRA 

Evaluator 

4 

8^ 

FIRST (1) 

muSIMP 

Selector 

7 


FLAGS 

ALGEBRA 

Global variable 



FLAGS (0) 

ALGEBRA 

Printer 

4 

8 

FOURONPI 

TRGNEG 

Constant 


1 

FREE (2) 

SOLVE 

Recognizer 


1 

FREE (2) 

DIF 

Recognizer 

1 

1 

FUNCTION 

muSIMP 

Keyword 

16 


GCD (2) 

ARITH 

Numerical 

2 

3 

GET (2) 

muSIMP 

Property 

14 


GETD (1) 

muSIMP 

Definition 

16 


HALF 

TRGPOS 

Constant 


1 

IDENTITY (1) 

ARITH 

Identity function 

2 

1 

IDMAT (1) 

MATRIX 

Constructor 

1 

2 

INFIX 

muSIMP 

Parse property 

35 


INT (2) 

INT 

Numerical 

1 

2 

INTEGER (1) 

muSIMP 

Recognizer 

10 


LBP 

muSIMP 

Parse property 

35 


LCM (2) 

ARITH 

Numerical 

2 

3 

LENGTH (1) 

muSIMP 

Sub-atomic 

18 


LINELENGTH (1) 

muSIMP 

Printer 

25 


LIST (N) 

muSIMP 

Constructor 

8 

lO 

LN (I) 

LOG 

Numerical 

1 

LOAD (3) 

muSIMP 

System function 

34 


2 


Name 


Descriptor 

D<?cuiBent„ed Dexin, 8 .d 

LOG (2) 

LOG 

Numerical 

♦ 

1 

1 

LOGARITHM (1) 

LOG 

Recognizer 


2 

LOGBAS 

LOG 

Control variable 


1 

LOGEXPD 

LOG 

Control variable 


1 

LOGEXPD (2) 

LOG 

Evaluator 


1 

LOOP (N) 

muSIMP 

Evaluator 

30 

30 

LOOP 

muSIMP 

Keyword 

30. 


LPAR 

muSIMP 

Constant 

35 


MAPPHN (2) 

EQN 

Mapping 


1 

MATCH (2) 

muSIMP 

Reader 

36 - 

22 

MATCHNOP (2) 

muSIMP 

Reader 

36 


MIN (2) 

ARITH 

Numerical 

2 

3 

MINOS (1) 

muSIMP 

Numerical 

18 


MKPROD (1) 

ARITH 

Constructor 


2 

MKRAT (1) 

ARITH 

Constructor 


5 

MKSOM (1) 

ARITH 

Constructor 


2 

MOD (2) 

muSIMP 

Numerical 

19 


MOVD (2) 

muSIMP 

Definition 

16 


MULTIPLE (2) 

ARITH 

Comparator 

2 

1 

NAt^E ( 1 ) 

muSIMP 

Recognizer 

10 


NEGATIVE (1) 

muSIMP 

Recognizer 

10 


NEGCOEFF (1) 

ARITH 

Recognizer 

2 

2 

NEGMULT (2) 

ARITH 

Comparator 

3 

1 

NEWLINE (0) 

muSIMP 

Printer 

24 


NOT (1) 

muSIMP 

Logical 

12 


NUM (1) 

ARITH 

'Selector 

,,3 . ... 

2 

NOMDEN 

ALGEBRA 

Control variable 


3 

NOMDEN (1) 

ALGEBRA 

Recognizer 


3 

NUMNUM 

ALGEBRA 

Control variable 

■■ 1 ' ■ .: 

1 

NUMNDM (1) 

ALGEBRA 

Recognizer 


1 

NUMBER (1) 

ARITH 

Recognizer 

3 

1 

OBLIST (0) 

muSIMP 

Constructor 

8 


OR (N) 

muSIMP 

Logical 

12 


ORDERP (2) 

muSIMP ■ 

Comparator 

11 


PARSE (2) 

muSIMP 

Reader 

36 

21 

PBRCH 

ARITH 

Control variable 

lr4 

10 

PI0N2 

ARITH 

Constant 


12 

PI0N4 

TRGNEG 

Constant 


1 

PLUS (2) 

muSIMP 

Numerical 

18 


POSITIVE (1) 

muSIMP 

Recognizer 

10 


POSMULT (2) 

ARITH 

Comparator 

3 

1 

POWER (1) 

ARITH 

Recognizer 

3 

1 

PREFIX 

muSIMP 

Parse property 

35 


PRIMES 

ARITH 

Global variable 

4 

11 

PRINT (1) 

muSIMP 

Printer 

24 


PRINTLINE (1) 

muSIMP 

Printer 

24 


PRODUCT (1) 

ARITH 

Recognizer 

3 

1 

PROPERTY 

muSIMP 

Keyword 

15 


PRTMATH (4) 

muSIMP 

Printer 

25 


POT (3) 

muSIMP 

Property 

14 
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PUTS (2) 
PWRpFD 


QUERY (2) 
QUOTIENT (2) '' 

RADfX (1) 

RBP 

RDS‘"(3) 

READ (0) : 

READCHAR (0 ) “ 
RECIP (1) 

RECLAIM (0) . 

REMPROP (2) ' 
REPLACE? (2)' 
REPLACER (2K 
REST (1) 
REVERSE (2) 
ROW'(l) 

RPAR 

RREST (1) 
RRREST U) 

SAVE (3) 

SCAN (0) 

SEC (1) 

SECOND (1) 

SIGN (1) 

SIMPU (2) 

SIN (1) 

SIN (1) 

SOLVE (2) 
SPACES U) 

SUB (3) 
SUBROUTINE 
SUM (1) 

SYNTAX (N) 

TAN (1) 

TERMINATOR (0) 
THIM) (1) 

TIMES (2) 

TRACE (N) 
TRGEXPD 
TRGEXPD 
TRGEXPD (2) 
TRGSQ 

UNION (2) 
UNTRACE (N) 

WHEN 
WRS (3) 


Muoule 

Descriptor 

Documented 


muSIMP . 

Definition 

16 


ALGEBRA;; 

Control variable 

2 


INT 

Reader/Printer 


1 

muSIMP 

Numerical 

19 


muSIMP 

Printer 

25 


muSIMP 

Parse property 

35 


muSIMP 

Reader 

22 


muSIMP 

Reader 

21 


muSIMP 

Reader 

21 


ARITH 

Recognizer 

3 

1 

muSIMP 

Storage 

33 


muSIMP 

Property 

14 


muSIMP 

Modifier 

9 


muSlMP 

Modifier 

9 


muSIMP 

Selector 

7 


muSIMP 

Constructor 

8 


ARRAY 

Recognizer 


1 

muSIMP 

Constant 

35 


muSIMP 

Selector 

7 


muSIMP 

Selctor 

7 


muSIMP 

System 

34 


muSIMP 

Reader 

21 


TRGPOS 

Numerical 

1 

1 

muSIMP 

Selector 

7 


INT 

Recognizer 


ARITH 

Evaluator 

3 

4^ 

TRGPOS 

Numerical 

1 

1 

TRGNEG 

Numerical 

2 

1 

SOLVE 

Numerical 

1 

3 

muSIMP 

Printer 

24 


ARITH 

Constructor 

3 

1 

muSIMP 

Keyword 

35 


ARITH 

Recognizer 

3 

1 

muSIMP 

Reader 

36 

22 

TRGPOS 

Numerical 

1 

1 

muSIMP 

Recognizer 


22 

muSIMP 

Selector 

7 


muSIMP 

Numerical 

19 


TRACE 

Debugger 

1 

1 

TRGPOS 

Control variable 

2 

1 

TRGNEG 

Control variable 

2 


TRGNEG 

Evaluator 

3 

3 

TRGPOS 

Control variable 

2 

1 

SOLVE 

Constructor 


1 

TRACE 

Debugger 

1 

3 

muSIMP 

Keyword 

31 

o 

muSIMP 

Printer 

25 


4 
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Name 


Descriptor 

Documented Defined 

ZERO (1) 

muSIMP 

Recognizer 


10 


ZEROBASE 

ARITH 

Control variablS^' 


1 

■ ■■ 

EERdEXPT 

ARITH ; 

Control variable^ 


1 


. 

muSIMP 

Comparator ^ ,, 


11 


> 

muSIMP 

Comparator 

.... H 

11 

. . -t. - - • 

> 

ARITH 

Comparator 



.. ..3./ 

< 

muSIMP 

Comparator 


11 


< 

ARITH 

Comparator 



.3.,r 

m 

• 

muSIMP 

Assignment 

^ >* '*'* 

13 

f , ■ r- * . r*' 


muSIMP 

Numerical 

>■ * 

20 . 

k} • . 


ARITH 

Numerical 





muSIMP 

Numerical I 


20 / 



ARITH 

Numerical 



' . : A. 


muSIMP 

Numerical • 


20 



ARITH 

Numerical 




/ 

muSIMP 

Numerical ' 


20 

- - ■■ 

/ 

ARITH 

Numerical 





ARITH 

Numerical 



■ ■ d■-* 

1 

ARITH 

Numerical 


5 

• la:--: 

*> 

• 

ARITH 

Error function . ' 



, „8-.: 

\ 

MATRIX 

Numerical : 


3 



MATRIX 

Numerical ’ ” 


2 


tARB 

SOLVE 

Constant 


1 

yyV'ii 

*E 

ARITH 

Constant j 


4 


#I ^ 

ARITH 

Constant 2 , 


4 

'J- 


ARITH 

Constant 

* ' 

4 







